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ABSTRACT 

Previous work and available theories concerning prestressed 
concrete members subjected to torsion, bending and shear are 
reviewed. The experimental study of this investigation consisted 
of tests on forty four eccentrically prestressed I section beams 
subjected to various proportions of torsion, bending and shear. Al] 
the beams had nominally the same lenath and cross sectional dimensions. 
The primary variables included two levels of average prestress, 
three eccentricities, three amounts of web reinforcement and 
various combinations of torsion, bending and shear. Only three 
beams did not have web reinforcement. Torque, transverse load, 
Strains in the web reinforcement and the strands, angle of twist 
and flexural deflections at selected locations were measured during 
the tests. 

Theoretical analyses were developed to predict the cracking 
and ultimate strength of prestressed I beams. Cracking strength 
was predicted by the following two methods. 
1). Torsional shear stress distribution was assumed to be in 
accordance with plastic theory. The stresses due to other forces 
were considered to be elastic. Principal stress failure criteria 
was used coupled with a tensile strength of concrete equal to 6vFr 
2) A general three dimensional finite element analysis using 
hexahedrons was developed. Cracking was assumed to commence when the 


maximum principal tensile stress equalled the modulus of rupture es- 
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timated at 7.5V7F¢. The out put of the program included maximum princi- 
pal tensile stress, initial crack inclination, angle of twist, flexural 
deflections, for any given loading. This analysis was applied to se- 
veral beams in the literature, in addition to the beams of the present 
investigation. 

Tests showed that the I beams fared in skew bending either in 
a bending mode or a torsion mode. Expressions were developed based on 
a skew bending mechanism of failure for each of the modes to predict 
the ultimate strength of I beams of this investigation and other in- 
vestigations. 

Comparisons were presented for theoretical and test results 
both at cracking and ultimate and the agreement was generally found 


to be good. 
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Dimensions and Section Properties 


A = cross sectional area of a beam 
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b, = width of the section at any height of beam 

e = eccentricity of prestress force with respect to 


centroid of cross-section 


h = total depth of a beam 

hy = height between center of longitudinal corner bars 

I = transformed-uncracked moment of inertia of a section 
Kx = depth of centroid of compressive force in concrete 


measured from the compression face of the beam 


Q = statical moment of the cross-section area above (or below) 
that level (where shear stress is computed) about the 
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Material Properties 


ES = modulus of elasticity of prestressing reinforcement 
E. = modulus of elasticity of stirrup steel 
EL = modulus of elasticity of concrete 
Ey = ultimate strain in concrete in pure bending 
ety = yield strain of stirrup 
in = compressive strength of concrete cylinder 
us = modulus of rupture of concrete 

fr, = tensile strength of concrete 
Sp = split cylinder tensile strength of concrete 
Fry = yield stress of stirrup 

my = yield stress of longitudinal mild steel 

py = yield stress of prestressing steel 


XX 


BRO! s2evoveneis mov? soatw2 wniticbeonr sited. 


phot dntog yp 08 nottose to wen st 


notensd Al vedi? emevdxs of noldoe2 to browtiaa wort aonedeth = 
«tod at vat? omeydxs of nordoar Fo bichitnad iat panera 


auyvit2 6 JO nOezrambb ont T-WiInedveprET = 


x 


A 


hh i ee a 


tnamsovntarey pniezsisae%q to Vitor iaefs to eufabom = 
[aad2 qinvidt? 20 Ytforsesla to evfubam = 

gteriono? to yitottzefs to euluhon = 

patbned awe nf stoyonog of} visite stemislu 
auwittie to ateste bisty 


il 


yebniiyo sisvsnes to dépnetwse ovizes gms 
sto1s00s to siwqus to zulebom = 

gteronos to ddpneyte sftanad = 

atsrano> to dsenetie slieney vwebarlyo drige = 
| quitite to 2vovt> blery = 

toate bfim fentbustpiof Yo ezsvse biaty + 
fasta pntezwes1y 0 eeorte biaiy < 


ly ’ as, 7 


Forces and Moments 


C = compressive force in concrete 

Bi = total effective prestress 

M = bending moment 

May = bending moment at cracking 

My = ultimate bending moment in combined loading 

Wee = ultimate shear-bending moment 

He = ultimate pure bending moment 

ber = transverse load at cracking 

ve = transverse load at ultimate 

if = twisting moment (torque) 

qh = bending component of applied torque 

ye = cracking torque based on elastic concept of concrete 
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in = cracking torque based on plastic concept of concrete 
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Ls = diagonal tension cracking torque in pure torsion 

Lee = diagonal tension cracking torque in combined loading 

Typ = ultimate torque of plain concrete members 
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Miscellaneous 


Oo = torsion factor in St. Venent theory 

B = inclination of the compression zone to the longitudinal 
axis of the beam 

By = empirical constant used by ACI 318-71 Code to define com- 
pressive stress block 

6 = loading ratio defined as T/V bY 

w = loading ratio defined as T/M 

) = inclination of the tensile crack to the longitudinal 
axis of the beam 

m = ratio of volume of longitudinal steel to volume of web 
reinforcement 

Py = web reinforcement percent based on total cross-sectional 
area 


XX111 


- 7 ' eeatee 
- . ~ - 
: . : 7 
" - 7 
_ : 7 


=noa'sntteb of pbod IT-SPE 10A' Ye 


fanotspe2-2e019 Tatod na bezed Steen trehaarra oe dow z 


- 


CHAPTER I 
INTRODUCTION 


1.1 General Remarks 


Extensive research with respect to torsion in reinforced 
concrete has resulted in the formulation of a design criteria for such 
members. The ACI building code (ACI 318-71) for the first time, now 
contains provisions for designing reinforced concrete members to 
resist torsion combined with bending and shear. Though research has 
been in progress since the late 1960's similar design specifications 
for prestressed concrete members have not been formulated mainly 
because of inadequate test data. The majority of torsion investigations 
in prestressed concrete beams have been carried out for rectangular 
cross sections; only a limited amount of work has been done on 
practical shapes such as I, T and box sections. 

Earlier studies of flanged sections of prestressed concrete 
(Gardner, 1960, Reeves, 1962) were of limited scope and dealt with 
I and T shaped sections without web reinforcement subjected to 
combined torsion and bending only. Zia (1961) conducted pure torsion 
tests of prestressed rectangular I and T beams with and without web 
reinforcement. The amount of web reinforcement was constant in his 
beams. Bishara (1969) described tests of prestressed rectangular, I, 
and T beams identical in cross-section to those of Zia (1961). These 


beams were tested in combined torsion bending and shear and no strength 
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predictions were proposed for the combined loading case. Gausel's (1970) 
investigation of prestressed I beams subject to torsion, bending and 
shear was limited to interaction studies. Wyss, Garland and 
Mattock (1969) made an experimental study of full scale I shaped 
bridge girders prestressed eccentrically with differing amounts of web 
reinforcement subject to pure torsion. More recently Wyss and Mattock 
(1971) conducted two series of tests on eccentrically stressed I 
girders for evaluating the interaction of torsion-bending and 
torsion-shear. 
1.2 Object and Scope 

Prestressed } beams used in bridges that are skew or 
curved in plan are subject to considerable amounts of torsion in 
addition to bending and shear. The test data and current knowledae 
about such beams under combined loading is very limited. The main 
objectives of this investigation were to obtain experimental data on 
the behavior and strength of prestressed I beams and to develop a 
Simple, approximate but reasonably accurate method for predicting 
both the cracking and ultimate strengths of such beams. A total of 
forty four prestressed I beams were tested to failure to evaluate 
the effect of such variables as level of prestress (two levels of 
average prestress), eccentricity (three eccentricities), amount of 
web reinforcement (three amounts of web reinforcement) and loading 
combination (five T/Vb., ratios). Measurements were made of torsional 
load, transverse load, strains in the web and prestressing 
reinforcements, angle of twist and flexual deflections at three 


locations in the test zone. 
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Cracking strength of these beams was predicted by a three 
dimensional finite element analysis and also by classical theories. 
Analytical expressions were developed based on skew bending mechanism 
of failure to predict the ultimate strength of prestressed I beams 
under combined loading. Comparisons were made of theoretical 
predictions and experimental results for beams of this investigation 


and several beams in the literature. 
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CHAPTER I] 
REVIEW OF PREVIOUS WORK AND AVAILABLE THEORIES 


Long neglected in comparison with other forces, torsion in 
concrete members assumed increasing importance in the 1950's; pre- 
viously torsional loads were considered to be secondary and were neg- 
lected as higher factors of safety of that time indirectly accounted 
for torsional effects. However, more recently, with reduced load 
factors and refined design methods it becomes mandatory to design for 
torsion explicitly. Further, in structural components such as edge 
beams, and girders curved in plan torsion is of primary consideration. 
This concern with torsion was reflected in the intensity of research 
work at various institutions resulting in significant contributions 
to the solution of the problem of torsion and its interaction with 
other forces. The accumulation of research findings to date is so 
vast that it is impossible to cover them all here and further, exten- 
Sive reviews of the literature with emphasis on rectangular beams 
have been conducted by various authors (Kemp et al 1961, Zia 1968, 
Mukherjee and Warwaruk 1970, Gangarao and Zia 1970, Chander, Kemp, and 
Wilhelm 1970, Henny and Zia 1971, Woodhead and McMullen 1972). 

In the sections that follow, first the cracking strength and 
the available solutions are discussed for concrete members presetressed 
or unprestressed under torsion alone or in combination with bending and 


shear. Later, research into the ultimate strength of prestressed rec- 
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tangular and flanged section beams in torsion, bending, and shear is 


covered. Finally, the space truss theory for torsion is briefly reviewed. 


2.1 Cracking Strength 
2.1.1 General 
Plain concrete beams under torsion alone or combined with 
other loadings fail abruptly at formation of the first crack. Providing 
reinforcement in such members in one direction alone does not yield 
any increase in strength beyond cracking. By providing reinforcement 
both longitudinally and transversely, the torsional strength can be 
increased substantially after cracking, however, the torsional stiff- 
ness reduces quite significantly once cracking occurs. In the case of 
prestressed concrete beams similar behavior is noted with respect to 
cracking. However, as the cracking strength is increased by the appli- 
cation of prestress for prestressed beams the level of cracking load is 
nearer the ultimate load. Hence, it is necessary to compute the strength 
of a beam at the onset of cracking, when torsional loads are applied 
alone, or in combination with other loads. 
Computation of the torsional cracking strength of a concrete 
beam includes the following three different phases. 
(a) Shear stress distribution across the section. 
(b) Failure criteria for concrete cracking. 


(c) Tensile strength of concrete. 


In the following sections these three phases are discussed 
and used to calculate the cracking strength of beams subjected to 


combined loadings. 
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2.1.2 Available Theories 
oslnerl  Grassicaliheor ies 

In early torsion research on concrete, two classical theories 
have been used to predict the strength of plain concrete beams sub- 
jected to torsion. These were the elastic theory (Saint-Venant, 1856) 
and the plastic theory (Nadai, 1911). In the elastic theory, concrete 
is assumed to be ideally elastic. In the plastic theory concrete is 
treated as a fully plastic material; the shear stress distribution 
across the section is uniform. 

Saint Venant's semi-inverse method of elastic solution of 
the torsion problem has limited application only to convex cross 
sections, without re-entrant angles as in the case of I, L, and T 
sections. However, for elastic shear stress distributionsof more 
complicated shapes the soap film analogy by Prandtl (1903) may be 
used. The governing differential equation of the deflected shape of 
a soap film stretched by a small differential pressure over an 
opening of the same shape as the cross-section of the member under 
torsion, has the same form as the differential equation of the tor- 
sional stress function. The differential equation of the soap film 
can be solved numerically to arrive at the shear stresses at discrete 
intervals on a cross section. 

A similar analogy, known as the Sand Heap analogy, was dis- 
covered by Nadai (1911) for the plastic torsion theory. If sand is 
piled on a plate having the same cross section of the member sub- 
jected to torsion the volume and slope of the sand heap are related 


to the applied torque and shear stress respectively. 
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Though the two theories adopt different assumptions regard- 
ing the material behavior, they use the same failure criteria, namely, 
principal stress failure theory in which the member is considered 
to have failed when the maximum principal tensile stress becomes 
equal to the tensile strength of concrete. Most previous researchers 
observed that plain rectangular beams failed suddenly coincident with 
the appearance of first crack and the failure surface intersected all 
four sides in the form of helical cracks. The observed mechanism of 
failure was in line with Saint Venant's theory. The strength equations 
for a rectangular concrete beam take the following form: 


Elastic Theory 
ee 
Tes Cbehnia ) (251) 
Plastic Theory 


Trazl-gpl] ony (2.2) 

Concrete however is known to have deformation properties 
intermediate between those corresponding to ideal elastic or ideal 
plastic materials. Hence neither of the two theories adequately re- 
flects the behavior of concrete beams under torsion. For a given 
limiting tensile stress, the elastic theory underestimates the strength 
by as much as 50 percent while the plastic theory accounts for the ex- 
cessive experimental strength, but fails to incorporate aspect ratio. 


The calculated elastic shear. stresses are susceptible to the shape 
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8 
of the cross section. Tests on concrete in torsion, of different cross 
sections as summarized by Marshall (1944) indicated that at failure 
plastic shear stresses were more consistent, less influenced by the cross- 
sectional shape, and, nearer to the so called tensile strength of 
concrete than the critical stresses computed on the basis of elastic 
theory. 

Reeves (1962) pointed out that in these theories maximum 
Stress was picked up and failure criteria applied without any regard 
for the stress gradient. He held that this method of approach was 
partly responsible for the general unsatisfactory results. Probably 
Hus (1968a) was the only investigator in torsion studies who took into 
account the strain gradient while deriving a relationship between the 
tensile strength and modulus of rupture, for his prediction equation 
of torsional strength. 

The inelastic behavior of concrete was believed to be re- 
sponsible for the increase in the strength over that obtained by the 
elastic theory. To account for this Miyamoto (1927), Turner and 
Davies (1934) and Marshall and Tembe (1941) using the non-linear stress- 
strain relationship of concrete in idealized forms, proposed torsion 
theories. These are called semiplastic theories as only partial 
plasticity of concrete was assumed. The strength computed by any of 
these theories should be intermediate between the predictions of 
elastic and plastic theories. However, the one advanced by Turner 
and Davies estimates torques higher than plastic torques whenever 
the height to breadth ratio exceeds 2.5 and thus becomes invalid. 


In any case reasonable correlation with test results was ob- 
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tained by various researchers using any of the above mentioned classi- 
cal theories coupled with appropriate values for the tensile strength 
of concrete. 
2.1.2.2 Skew Bending Theory 

Hsu (1968a) proposed a new skew bending theory for plain 
concrete beams, rejecting the classical theories on the ground that 
(i) the principal stress failure criterion might not be relevant to 
torsion problems and, (ii) the failure surface was not intersected on 
all the four faces of a rectangular cross section by a helical crack. 
Based on experimental observations, Hsu concluded that torsion failure 
is a bending phenomenon involving the development of a crack on three 
Sides and a compression hinge adjacent and parallel to the long side 
of the section but inclined to the axis of twist. Hsu's torque equation 


is expressed as 


bh 


{reas (0x85 fen) (2-32) 
An examination of equation (2.1) and equation (2.2) reveals that Hsu's 
bending theory differs with classical theories in two aspects. First, 
his bending theory uses a constant of 1/3 for the torsion shape factor, 
while for the other theories, it is a function of aspect ratio (h/b). 
The 1/3 corresponds to the highest value of the shape factor in the 
elastic theory. Second, in the bending theory and the classical 
theories, the principal tensile stress should respectively exceed 0.85 
times the modulus of rupture, and the tensile strength of concrete to 


cause failure. However, the difference between 0.85 Le and f, becomes 
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negligible in beams of practical sizes since the strain gradient for 
them is reduced from that of test specimens. Equation (2.3) was modi- 
fied by Hsu to involve the concrete compressive strength instead of 


modulus of rupture and is shown below. 
"A ule (2.4) 


In proposing his skewed bending theory, Hsu idealized the failure 
section as a plane, although it is actually a surface resembling a 
warped trapezium. Only the bending component of torque was considered 
in satisfying equilibrium, while the twisting component of torque on 
the failure surface was presumed only to cause a reduction in the modulus 
of rupture. He computed flexural tensile stress elastically at the 
ultimate stage. In spite of these defects, Hsu's skew bending theory, 
in using a torsional shape factor of 1/3 and modulus of rupture instead 
of other types of tensile stresses, improves on the classical elastic 
theory in that his strength predictions showed good correlation with 
experimental results. 
Zalecdso Statistical Approach 

All the preceding theories yield prediction equations of the 
same form shown below though in their specific assumptions they differ 


significantly. 
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where K = A dimensionless coefficient called shape factor which 
is a function of h/b in classical theories and a con- 


stant equal to 1/3 in bending theory 


Mukherjee and Kemp (1967) attempted to generalize the different shape 
factors of the theories to yield one expression. They carried out a 
statistical regression analysis of available test results and deter- 
mined the parameters of the shape function for concrete rectangular 


members under pure torsion. Their equation can be stated in the form: 


T= 0.4124 [1-0.23332] b°h £! (2.6) 
u h t 
where le = pure torsional strength of rectangular plain concrete beam 
Shape factor K = 0.4124 [1-0.2333 P43 


fincebvite 


Using equation (2.6), for slender beams the calculated torque capacity 
was closer to the elastic torque, but approached the plastic value at 
the height to width ratio equal to unity. 

Subsequently, based on statistical observation of test re- 
sults from different investigations Chander, Kemp, Wilhelm (1970) 
presented the following prediction equation. 
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where fi = by f. 
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This equation yields results much closer to those of the plastic theory 
rather than elastic and correlates well with test data. It shows that 
concrete though not capable of complete plasticity, need not be completely 
plastified as the contribution to torsional resistance is maximum from 
the outer fibres and only their plastification is necessary to yield 
results close to those obtained by the plastic theory. 
2.1.3 Flanged Concrete Beams Under Pure Torsion 

Comparatively little research work has been carried out on 
members of practical shape such as I, L and T sections subjected to 
pure torsion. Contributions to the knowledge of the behavior of the 
flanged concrete beams include those of Turner and Davies (1934), 
Marshall and Tembe (1941), Nylander (1945), Zia (1961), Arockiasamy 
(1964), and Ramakrishnan and Jayaraman (1968), Hsu (1968), Evans, Kemp, and 
Wilhelm (1970). Some of these investigators computed the strength 
of component rectangle sections and used Bach's (1911) approximation 
as detailed below to estimate the pure torque capacity of T and L 
sections. It may be recalled that the elastic shear stress distribu- 
tion and the torsional strength can be obtained by membrane analogy. 
Torsional strength is proportional to the volume under the membrane . 
Bach's approximate method computes the volume, as the sum of the volumes 
of the individual component rectangles, treating the membrane of each 
rectangle as a uniform cylindrical surface. Thus, it neglects the 
end effect and junction effect which results in gain and loss of volume 


at the ends of each rectangle and junctions of any two rectangles, 
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respectively. For a section composed of narrow rectangles, the approxi- 
mate method yields a volume close to that of the membrane analogy. But, 
in concrete sections the component rectangles are not narrow and Bach's 

method underestimates the volume and hence the torque capacity. Depend- 
ing on the relative dimensions of the cross section this can be Significant. 


Bach's equation for the maximum shear stress is 


Tmax snck (2.8) 
eee maximum torsional shear stress 
‘i = applied torque 
L = largest thickness of component rectangles 
b = width of each component rectangle 
h = depth of each component rectangle 


Plain concrete flanged beams (excluding the I sections) fail 
at the appearance of the first crack, thus the cracking and the ulti- 
mate torsional moments are the same. Recently Hsu (1968b) conducted 
tests on eight plain concrete T and L section beams subjected to pure 
torsion. He observed that non-rectangular beams under torsion fail 
by bending, and the calculation of failure torque based on his bend- 
ing mechanism is difficult. Extending his equation (2.4) derived 
for rectangular beams, and using Bach's approximation, he presented 


the following equation for the strength of flanged beams. 
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14 
He stated that the outstanding flange width was ineffective if it ex- 
ceeded about three times the flange thickness. 

The skewed bending theory was applied to predict the strength 
of symmetric I sections tested and reported in this investigation and 
was found to be inadequate. While Bach's method used by Hsu to compute 
cracking of flanged beams does not consider the relative position of 
rectangles, the bending theory takes into account the position of an 
element on the section in terms of moment of inertia. This results in 
a disproportionate increase in moment of inertia, and hence a similar 
increase in failure torque, if the cross sectional area farthest from 
the bending axis increases slightly, and it may be concluded that his 
theory is of limited use in flanged beams. 

In a recent study Evans, Kemp and Wilhelm (1970) investigated 
into the strength and behavior under pure torsion of 33 T-beams and 5 
L-beams either plain or reinforced. The initial crack always occurred 
in the web (web was 6 in. thick, while flange thickness was either 
3 in. or 4 in.) leading to immediate failure if the beams were plain 
or reinforced only in one direction. All the specimens having steel 
in both directions exhibited very large amounts of ductility and a 
concrete compression zone was observed only when a specimen was deformed 
well past its ultimate torsional strength. 

Evans et al disputed with the rationality of the ACI's pro- 
vision of associating flange effectiveness in resisting torsion to the 
ratio of the outstanding flange width to its thickness; and limiting 
such a ratio to three. Their test results showed that even at a ratio of 


five a limit had not been reached. However, they observed that if the 
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maximum tensile stress occurred in one of the component rectangles 
failure might occur before the other component rectangles developed 
their full capacity. It is generally believed that Bach's approximate 
method is adequate and yields conservative estimates of torsional 
strength for the flanged beams. Utilizing Bach's method, they com- 
puted the cracking strength of a flanged section on the premise that 
as the applied torque gradually increased each of the component 
rectangles undergoing the same rotation carried a portion of the total 
torque proportional to its torsional stiffness. The applied torque 
could increase until one of the component rectangles developed a 

crack while the others might not have been stressed to their capacity. 
The final cracking torque was then equal to the sum of the torques 
carried by each of the component rectangles at a rotation corresponding 
to the cracking rotation of that particular component in which the 
initial crack appeared. Additionally, junction effect had to be 

taken into account. This method of computation was similar to the 
moment distribution and considers the effectiveness of the component 
rectangles without restricting their size. They pointed out that the 
strength calculation by Bach's method for some large T beams in the 
literature was unsafe and the effective section concept yielded safe 


predictions. The equations derived by them are as follows: 
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effective ae (2.11) 
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However for design purposes, Evans et al stated that the ACI equation 
with its limitation of the flange width to depth ratio seemed adequate. 
2.1.4 Prestressed Concrete Beams Under Pure Torsion 

Nylander (1945) was the first to carry out an investigation 
into the torsion problem of prestressed concrete members. He tested 
a series of sixteen 8-inch square specimens under pure torsion. The 
prestress was applied in the form of axial compression. It was observed 
that the precompressed beams exhibited torsional strength about double 
that of their unprestressed counter parts at the first crack. The failure 
was similar to that observed in uncompressed beams. He concluded that 
the ‘cracking torsional strength of his beams was adequately predicted 
by the plastic theory in conjunction with principal stress criteria 
of failure. 

Cowan and Armstrong (1955) made tests on prestressed beams 
in pure torsion as well as in torsion and bending. Some beams had 
uniform prestress, while others had nominally triangular prestress 
distribution. They observed that the failure was sudden at the for- 
mation of the first crack in the case of beams subject to pure torsion 
or high torque/moment ratio. The addition of bending moment lessened 
the violence but increased the destructiveness of failure, and allowed 
an increase in load to be carried after the formation of the first 
crack. Their 6 x 8 in. eccentrically prestressed beams cracked on 
the top first in pure torsion. The equation proposed for the per- 


missible strength of a rectangular prestressed concrete beam, My is 
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MpaiMy Clee areuy/? (2.12) 


where Moo = maximum permissible twisting moment for a similar 
plain concrete section 


Fo = magnitude of prestress 


The experimental work reported by Humphreys (1957) related 
to the testing of 94 plain prestressed rectangular beams, four of which 
were eccentrically prestressed with zero stress at the top. Failure 
in the beams was so sudden that tensile cracking could not be detected. 
In the case of eccentrically prestressed 5 in-square and 10 in. x 5 in. 
rectangular beams the initial crack occurred on the top face. 

Humphreys stated that, though the crack occurred at a loca- 
tion of zero prestress on the top in rectangular beams, the strength 
was governed by the same conditions as that for a concentrically 
stressed member, because of the higher torsional factor associated 
with the short side of the section. The square and rectangular beams 
developed a torsional strength 13 and 43 percent respectively in excess 
of their unprestressed counterparts. He approximated the principal 
tensile stress theory of failure by a regression equation which yields 


the following ultimate torsional stress. 


med ays Ty + 0.307 Op (2.13) 


where i = ultimate torsional shear stress. 
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His equation for the torsional strength prediction of a plain pre- 


stressed member is of the following form. 


where 


iT] 


iT] 


iT] 


T = k(r#0.30,) beh (2.14) 


ultimate torsional strength 
a function of depth to breadth ratio 


ultimate shear strength of a plain concrete member 


Based on the tests in the literature, Hsu (1968c) extended 


his skew bending theory to concentrically prestressed rectangular 


beams without web reinforcement. The torsional strength of such a 


beam was obtained as the torsional strength of non-prestressed beams 


multiplied by a prestress factor which accounts for the effects of 


prestress. 


where 


The equation is of the following form. 


T,, = [e(b+10)h #2'/39p7 + Ley t/2 (2.15) 
Cc 


= ultimate torque of plain prestressed beams 


(1 + aoe le = prestress factor 


Cc 


In establishing the prestress factor, the modulus of rupture was taken 


equal to f./8.5. This prestress factor becomes identical to the factor 


that would be derived assuming a principal tensile stress failure 


criterion and a value of tensile strength equal to Fi/10. From the 


test results of Humphreys (1957), Hsu concluded that a compression 


failure would occur if average prestress exceeded 70 percent of the 
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19 
compressive strength. The applicability of equation (2.15) was limited 
to eccentrically prestressed sections which cracked on the wider face 
and in such cases the crack was assumed to occur at the centre of the 
wider face, neglecting the variation of prestress. This theory could 
not estimate the strength of nearly square beams which cracked first 
on the top face. 

Mukherjee and Kemp (1967) extended their statistically derived 
prediction equation (2.6) for rectangular concrete members, to include 
the effect of concentric prestress. The torsion factor of this ex- 
pression contains a term involving prestress which resulted in torques 
increasing slightly with the level of prestress. The enhancement of 
torsional shear resistance due to the presence of prestress was taken 
into account by using the maximum principal stress theory. The re- 


sulting expression is as follows. 


Dee ae pes ae 
i = [0.4124(1-0.2333 p) + 6 Fi b ht, (2816) 
where i = M1 + 0/f!3 ae 
pase nyt 


Equation (2.7) was modified in a similar manner to account for the in- 
creased strength of plain concrete members due to prestress. This 


equation for a plain prestressed beam may be written as 


2 
- b b 2 
i. = [0.4731(1-0.5924 p * 0.2763 2) b ine, aay) 
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where Testy (140,/F"4) 
fi=5/f! 


Chander et al observed from test results that the torsion factor was 
independent of degree of prestress. It was stated that tensile failure 
occurred if the average prestress did not exceed ten times the tensile 
strength of concrete. Equation (2.17) was checked by results of beam 
tests of several investigations in addition to their own and the corre- 
lation was found to be good. 

In the case of eccentrically prestressed rectangular beams 
under pure torsion the initial crack was observed by Humphrey (1957) 
and Zia (1960) on the top face whereas Mukjerjee and Warwaruk (1970) 
and Gangarao and Zia (1970) observed it at the top of the longer side, 
and Chander et al (1970) and Woodhead (1972) at the centre of longer 
Side. It was generally assumed that the principal tensile stress 
situation at the least compressed face determined the strength of 
such beams. Each of the above investigators observed experimental 
strengths in excess of that corresponding to a plain concrate beam. 
In some investigations, eccentrically prestressed beams exhibited the 
Same strength as concentrically prestressed beams having the same 
average prestress. 

With zero prestress at the top, Hsu's bending theory would 
predict a strength equal to that of a plain concrete beam if the first 
crack originates at the top of the vertical face. However, his theory 


cannot be applied if the top face cracks first. 
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Plastic theory cannot be applied in the usual form, as the 
Shear stress would depend on the level of prestress at each point and 
it varies nonlinearly on the cross section. 

For eccentrically prestressed beams only the elastic theory 
because of the higher torsion factor at the centre of the short face 
rather than on the long face predicts a higher strength than for a similar 
unprestressed beam. However, Chander et al (1970) observed that the 
elastic theory explained only a part of the increased strength and 
Showed that a skewed sand heap analogy resulted in better strength 
predictions of such members. They assumed that failure occurred only 
when the complete section was plastic. The shear strength varied 
in accordance with the prestress. With no tensile stress at the top, 
the slope of the sand heap at the middle of the longer side controlled 
the failure. But for practical use they recommended the use of inter- 
action diagrams as suggested by Mukherjee and Warwaruk (1970) for 
computing the torsional strength of eccentrically prestressed beams 
under torsion, treating the eccentricity of prestress as an external 
moment. 

Zia carried out an extensive test program consisting of pre- 
stressed and unprestressed rectangular, T and I sections subjected to 
pure torsion. He believed that the elastic theory was more justified 
than the plastic theory because for short time loadings high strength 
concrete such as used in prestressed concrete structures exhibited 
very limited plasticity. Elastic shear stress distributions for rec- 
tangular, T and I sections were obtained using finite difference techniques. 
The increase in torsional shear stress due to prestress was computed 


uding modified Cowan's failure criteria developed by Zia. 
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Failure criteria in general require a knowledge of both com- 
pressive and torsional strength of concrete. Generally, a premise is 
made that the torsional strength is equal to the tensile strength of 
concrete. As for tensile strength the following relationship based 
on the test results of Gonnerman and Shuman was adopted by Zia. 


3/4 (2.18) 


fy = 0.68 (0) 
Rectangular and T section beams failed immediately after cracking. 
The agreement between the measured and theoretical cracking results 
of these beams was good. Zia observed in his testing that the cracking 
and ultimate strengths of plain, concentrically and eccentrically 
prestressed I sections without web reinforcement were quite different. 
For these beams he reported only test results at ultimate but not at 
cracking. Comparison was made between theoretical crackings strengths 
and experimental ultimate strengths. Thus it is not possible to deter- 
mine how well the elastic theory compared with the first cracking 
load. But the fact was established that I beams of his investigations 
resisted substantial torques after cracking. This behavior of un- 
stirruped I beams, different from rectangular and T beams, depends 
strongly on the relative proportions of the flanges and the web, 
and if prestressed also on the eccentricity of prestress. Ina plain 
or concentrically prestressed beam, if the dimensions are such that 
first cracking occurs at the center of flanges then it spreads into the 
flanges and web and the beam collapses immediately. However, if the 


web is thicker than the flanges the more highly stressed web picks up 
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the first crack and it progresses towards the flanges at a slower rate 
because of the presence of less stressed flanges on either side of 
the web. Zia stated that at this stage a redistribution of stress 
takes place in the flanges transforming the torsional shear in them 
to lateral bending shear. The external torque is resisted by the internal 
couple of the lateral shear forces developed in the flanges. In this 
case failure occurs when the lateral shear force equals the lateral 
shear strength of the flanges. Zia's beams fall in this category 
having a web thickness of 3-1/2 in. and a flange depth of only 2 in. 
Such behavior for example, was non existent in the plain I girders 
tested by Wyss et al (1969) although the first crack occurred in the 
web: here the web and flange thicknesses were equal and failure 
occurred at the formation of the first crack. 

Irrespective of the thickness of the web relative to that 
of the flanges if the I beam is prestressed eccentrically Zia observed 
that the bottom flange became stronger than the top flange. When the 
top flange cracked first and the crack propagated towards the web the 
uncracked portion became a T section with the prestress more uniformly 
distributed on the uncracked portion. Such a section resisted torque 
in excess of the torque that caused top flange cracking. Unfortunately 
as pointed out earlier experimental cracking torques were not reported 
to assess the strength available between cracking and ultimate stages. 
If the experimental torques at the elastic limit in the I beams were 
assumed to be the test cracking torques, then the failure and cracking 


torques differ by about 11 to 40 percent depending on the eccentricity. 
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The cracking torques might have been greater than those at the elastic 
limit; then the difference between the strengths at cracking and ultimate 
would be less. 

|Ihough Zia did not draw any conclusions, his plain pre- 
stressed rectangular beams, and I beams showed essentially the same 
strength irrespective of the eccentricity. However in the case of 
plain prestressed T beams eccentricity of prestress had an adverse 
effect causing a decrease in strength. 

A major test series consisting of 18 large concrete I girders 
and six rectangular beams under pure torsion was conducted by Wyss, Gar- 
land and Mattock (1969). Twelve of the I girders were eccentrically 
prestressed. In an attempt to estimate the cracking strength of 
these beams they applied a two dimensional finite element method. 
Triangular elements were used with nodes located at the intersections 
of a 1/2 x 1/2 in. mesh. The output from this program consisted of 
a shear stress distribution which is essentially elastic. Cracking 
was assumed to occur when the principal tensile stress at a critical 
location exceeded the measured split cylinder strength. For most of 
the I girders the measured torques were considerably larger while for 
the rectangular beams they were consistently smaller than the computed 
values. 

In the case of web cracking, the calculated values compared 
well with the test results, if the high shear stress gradient in the 
vicinity of the web-flange juntion was neglected, and a smaller value 
of shear stress which existed at a distance 2 in. away from the 


top of the web and along the remainder of the web was adopted. With 
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regard to girders that cracked on the top, they stated that shallow 
cracks might have occurred at lower torques than those reported, 
but were not seen because of the rough trowelled top face. But no 
explanation was given for the opposite trend of lower test strengths 
than the computed ones in the case of rectangular beams. 
2.1.5 Reinforced Concrete Beams Under Combined Loading 

A beam under pure torsion develops equal tensile and com- 
pressive stresses on two perpendicular planes. The addition of a 
normal stress due to the application of any of the following loads 
prestress, bending, and direct tension or compression, further com- 
plicates the torsion problem. In such a case the maximum principal 
tensile stress is accompanied by normal stresses of varying degrees 
on orthogonal planes. Viewed differently, the shear stress distribu- 
tion is to be obtained under the simultaneous action of a known direct 
stress. This does not seem to be possible. Assuming the stress 
situation is accurately known it remains to be determined when failure 
occurs. The strength of a material is dependent on the state of 
stress. Though several mathematical failure criteria were advanced 
in the past for the combined stress problem experimental verification 
is very difficult in the case of concrete. 

The differences in the cracking analvses of various re- 
Searchers are as follows. 

First, elastic stresses due to bending moment, transverse 
shear, and prestress were separately obtained. The only difference 
related to the evaluation of the torsional shear stress distribu- 


tion. For this purpose, any one of the following methods was used: 
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(a) Elastic theory 

(b) Plastic theory 

(c) Semiplastic theory 

(d) Hsu's skew bending theory 


(e) Statistical approach 


Second, most of the investigators used a principal tensile 
stress criterion of failure. It was assumed that cracking would com- 
mence when the maximum principal stress due to the combined action 
of bending, shear, torsion, and prestress reached the tensile strength 
of concrete. However, very few used a failure criteria other than 
principal tensile stress. 

Third, differences are considerable, as also noted earlier, 
in assuming a value for the tensile strength of concrete. 

Most recently Barton and Kirk (1973) carried out a cracking 
analysis of their reinforced concrete T section beams subject to 
combined loadings assuming full plastic distribution of torsional 
shear stress and elastic distribution of bending and transverse 
shear stresses. They used the maximum principal stress failure 
criteria. The tensile strength of concrete was limited to 5 vie 
Two locations on the cross section, one at the neutral axis, the other 
either at the centre of top or bottom face, were examined for possible 
cracking. Ratios of test and analytical torques varied from 0.82 
to 1.58 with a mean value of 1.13. It is to be noted that the moment 
used in the cracking strength equations referred to the centre of the 


test span. 
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2.1.6 Prestressed Concrete Beams Under Combined Loadina 

As a first part of an extensive test program, Gangyarao and 
Zia (1970) reported the research findings of tests on prestressed rectan- 
gular beams (6"x12") under combined bending and torsion. Of the 42 
specimens 28 were prestressed concentrically and the rest eccentri- 
cally. All beams had transverse reinforcement and some of them were 
provided with longitudinal mild steel while others were not. They 
presented a cracking analysis using the skewed bending theory. Bend- 
ing was assumed to take place across a plane inclined at an angle to 
the longitudinal axis of the beam. The width of the failure plane 
was the length of the crack on the tension side of the beam. The 
initial crack was considered to occur when the internal stresses due 
to prestress, bending moment, and torque exceeded a reduced modulus 
of rupture. Stress due to flexural moment was computed elastically, 
and only the bending component of torque was considered to cause 
bending and contribute to the tensile stress. Two cracking modes, 
namely, bending mode and torsion mode were examined. In the former 
the tensile crack occurred in the bottom face while in the latter 
it appeared either at the top, bottom or centre of the wider face. 
Reduced modulus of rupture was approximated to f(/10. Finally, an 
equation was established involving torque and angle of cracking as 
unknowns. Minimization of torque with respect to crack angle pro- 
vided a second equation for solving the cracking torque and angle. 
Gangarao and Zia stated that elastic torsion theory did not work wel] 
and proved to be highly conservative. Beams subjected to pure torsion 


developed an initial crack at the centre of the wider face if con- 


sn st ns hl tt aie 
yout son vow evento af tw Tanda WFR Tantbutenot tty babtvayg | 
-baed .wioortt onticos bowoe ait ntan atavtans potions « batman, 


ot afons 6 s6 beaklont ons 6 eetine aaela ode? od bomueas 26m Qe? 


shalg avulist add Fo dfbtw oT mead Sag VO 2hxs tentbus tenet ott 


eit .magd ott YO abie noranat ails io" tomis ste 20 dapnet ant z6w 
sub eaeeattc [anratt? and netw wo90 0% Hevsbtanos zew does fetsint 
aufubom baoubey © bebesaks supros bas ,Inemta gatbied ,2eerteerq 03 


cy{feortesta bevsamas Saw Snead tesuxet? of sub agate .cxwdqus to 


ssus> of berebtemes 25w eupriot 2 Disnoqnes onbbnad ef? yline bas 
eeobon pridssts owl .22eette aftensd ody od studtydnos brs pntbned 
yam? ont nt ban tmexe tom eben notett bre esbom pntbaed -fomen 
wagneT od? nt Si tAw 938? motted off nF bownos0 doers sltenot aid 
gpa? dubtw adf Yo Sxdde2 Yo modtod .aod ard te teddts beyeoqge IF 
AL ath tent /OT\3Y od betemtxosaus 2ow syudqus to 2ufubon beoubsst 
26 entaasta te ipté bris supres pnivfovat boretidstee 26w not seups, 
OG ofpne Aosta Ot toaqgesy Atiw eupiot to notisstmiarh -2rwoniny 
-efpns bas syprad gntdsens ans patvfor yeh nofyeups bnosez 5 bebtv 

| [iow Seow ton brb yrostd wotzvet ofdesie tate batese ats bas osrspned 
notenos Seuiq of besosinua amead <jibipioane> idotd ot 08 enone be 
-09 Tt 9a8t ‘vabiw sid Yo etdnde afd I6 Apets Tetstnt ie 


Prey | . ; 


7 


28 


centrically prestressed, at the top if eccentrically prestressed. 
Beams with eccentricity cracked at slightly lower torques than the 
beams without eccentricity. 

Subsequently Henry and Zia (1971) investigated the behavior of 
32 rectangular prestressed beams (6"x12") loaded in combined torsion, 
bending and shear. All the beams were provided with longitudinal mild 
steel and the same amount of web reinforcement. One half of the total 
beams had concentric prestress and the remaining one half had eccentric 
prestress. The cracking analysis presented was identical to that deve- 
loped by Gangarao and Zia, but included the effect of transverse shear. 
In accordance with the elastic theory bending shear stress was assumed 
to be parabolic across the height of the section. They observed that 
the splitting tensile strength was an adequate measure of tensile stress 
as similar qualitative biaxial state of stress existed in the split 
test and the beam at initial cracking. Two cracking modes, bending 
and torsion, as in the earlier study, were identified. Two locations 
one at the centre and the other at the bottom of the side face on which 
Shear and torsion stresses were additive, were checked for cracking in 
the torsion mode. The possibility of an initial crack at the top of 
the long side was not considered. 

Flexural cracks were observed first near the bending load 
on the side at the bottom, even though the beam failed subsequently 
in a torsional mode. In some beams the initial crack formed at the 
centre of the side face. Henry and Zia observed that the cracking 
ave were very sensitive to loading combinations. The torque 
moment ratio used in the cracking analysis related to the failure 


section rather than the first crack section. This may then result 


to solvated ant bssaprteovet cisen) sf bne wna yisnsupsedu2 

-fotened bentdmo> st bsbsof isi . BLL 
bit fentbuttpnol mtw vettverq ane vsarz baw entbned 
fete sit Yo Vat snd jnansavotntey daw 0 tewoMs ame siz bas este 


sprinaone bed Pied ono gnintanen sd brie eagtdeeg atainsonos bat teed 


_sveb dada of (eottnabt caw besnsgetq, 2hzytens @rhlcen> oft .2zensestq 
sola sevevenent Yo dasits aid bebulont dud «ets bas omtepmed yd: beget 


bemees zow 22nvie vssd2 ontbnat wioens aidealo on? 3tw somabron36 1 . 


test bovuezdy ySAT = .nehioet ods. 10 sdptest ott 220726 2tlodevteq ad 03 
sseste offanst 19 ssueesm seupaie ns cow AdeNsTae gitensa entzstige sd3 
ttfga ote wt boaette epee! to etate Tetnatd aut pet Tsup vertmt2 26 
confined .zsbom yntdoers owl patanein fstatet 36 maod sd bas test 
enokiesaf owl .bettitnght sew. ybuse yetfags ont af 28 .notevos bas 
dntdw io sont Shite 249 to mostod of? 36 santo ofg bas sttnso ond 76 ano 
nt gatdss. st baioary avow ,sviatbbe svew eeteese notzi0s brs wesde 
to god Sd3 te dogo Isitiat as 40 yoiitdt2zog oT .sbom note703. 93 
-peiabtanoo tom eew sbfe paol oft 

beoT patbnad ont ween dz1tt bavysedo oie 242879 Terwxeld 
ultnsupezdue batrst mead edt fiouods nave .motsod edd jn. sbt2_ ant, no 
si? 26 bamvoy Az673 feldint o> emsed smaz at -sbom Tsnofewod & nt 
pittvasy edt tert bovvezde ofS bre yrs aos? sbte att to sxtna2 
suport sit anotsenidmos entbnol ot ev}ztzns2 vray ae 2beo! 
sviha? 96 ot batefor etavisns patiser six nt beew often, sneman 
dTueey neat vem etdT noises das Sevit ond nsrit toddey noldase 


\4n° aibnl ; 


29 
in computing a cracking strength at a section different from where 


the crack actually had formed. 

More recently Johnston and Zia (1971) conducted a series 
of tests on 37, 12-in. square hollow eccentrically prestressed beams 
in various combinations of torsion, beriding and shear. All beams ex- 
cept five contained web reinforcement. Stirrups were spaced either 
at 3 in. or 6 in. Two shear span to depth ratios were adopted. 

Theoretical analyses were presented to predict cracking 
and ultimate strengths. They compared their test results at crack- 
ing with two methods of analyses. The first one based on skew bend- 
ing theory, was presented by Gangarao and Zia; suitable modifications 
were applied to account for the hollow portion, presence of shear, and, 
the addition of the top cracking mode. The second one was elastic 
theory developed by themselves. A finite difference solution of the 
governing differential equation yielded torsional stresses which were 
combined with other stresses computed elastically. The initial crack 
was assumed to occur when the principal tensile stress exceeded the 
splitting tensile strength of concrete. The justification offered 
for the use of the splitting tensile strength was the insignificant 
strain gradient for a thin walled hollow beam as is the situation 
in the splitting test. Three cracking modes were analysed as the 
beams of this investigation cracked either on the bottom (Mode 1), 
or on top (Mode 3), or on the side (Mode 2). 

Johnston and Zia showed that the bending theory failed to 
forecast the cracking strength of hollow beams adequately, yielding 


unconservative results for Mode 2 and Mode 3. But it gave good corre- 
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lation with test results for bending (Mode 1). This is not unexpected 
as cracking in Mode 1 is mainly due to bending. The elastic theory 
agreed well with test results though its predictions proved to be 
slightly unsafe for Mode 2. 

To estimate the cracking strength of rectangular prestressed 
beams under combined loading Woodhead and McMullen (1972) used the 
elastic theory to compute the torsional, flexural shear, and bending 
stresses. Cracking was presumed to commence when the principal 
tensile stress equalled the tensile strength of concrete which was 
assigned a value of 7 VFo. Three critical locations, namely centres 
of the top, bottom and vertical faces, were examined. Several rec- 
tangular prestressed beams from different investigations were analysed 
and the correlation with the test results was good and ranged from 
1.035 to 1.33. For flanged sections they suggested the same method 
but did not analyse any beams even though the strengths at first 
cracking in the case of prestressed I girders by Wyss et al (1969) 
were available. 7 

Gardner (1960) reported tests on 16 plain I-beams nearly 
concentrically prestressed (eccentricity was only 0.2 in.) loaded 
under combined bending and torsion. He reported only the moments at 
the ultimate stage. Further, substantial bending moment was applied 
first, in some beams causing cracking, before any torque was applied. 
No beams were tested under pure torsion. He found that for the beams 
which did not develop tensile stresses on the application of bending 
moment, the limit of elastic behavior (defined as the limit of the 


linear part of the torque-twist curve) could be estimated using the 
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elastic theory. A finite difference method.was used to determine the 
Shear stress distribution. 

The prestress used related to the centre of the web in spite 
of the eccentricity of the prestress. A value of 0.707 times the 
measured modulus of rupture was taken as the tensile strength of 
concrete. Another calculation assuming full plasticity at the ulti- 
mate stage, using a uniform prestress and shear stress distribution, 
principal tensile stress failure criteria, and tensile strength equal 
to that used in predicting the limit of elastic behavior yielded 
strengths very close to the ultimate torques. 

Reeves (1962) described the testing of 42 plain prestressed 
concrete T section beams under combined bending and torsion. The 
prestress alone resulted in nominal stresses of 0 psi at the top and 
2000 psi in the bottom. The load application was similar to that 
adopted by Gardner (1960). In some beams there were cracks due to the 
bending moment before the final torque was applied. Strengths were 
reported at the limit of elastic behavior and at the ultimate stage. 
He attempted to interpret the ultimate strength using the classical 
theories and concluded that the elastic and plastic theories could 
not adequately predict the torsional strength of his T sections at 
failure. If the section is cracked it is not possible to determine 
the torsion factor, the prestress distribution, and the critical shear 
strength. 

Bishara reported tests on 24 prestressed rectangular, I 


and T beams subjected to combined torsion, bending and shear in various 
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combinations. Half of the beams were eccentrically prestressed. 
Rectangular beams developed initial cracks on the longer side where 
the torsional and flexural shear stresses were additive. These beams 
failed in a skew bending mode with the compression hinge on the side 
opposite to that where cracking started. In flanged beams cracking 
started randomly either on the web, or on the top face or at both 

the locations. It was observed that the flange(s) and the web appeared 
to behave as individual component rectangular members with regard to 
the extension of cracks to the adjacent sides implying discontinuity 
of cracks at the web flange junctions. However, failure finally took 
place by crushing of the concrete along a part of the web and the 
flange. 

Bishara used the elastic theory and observed that it failed 
to predict the cracking strength of prestressed rectangular and flanged 
beams but did not suggest any alternate method. He obtained the stress 
distribution due to all the forces including the torque by the elastic 
theory. Principal tensile stress failure criteria was used. The 
principal tensile stresses at the critical locations varied from 400 
psi to 1260 psi for rectangular beams and from 580 to 1600 psi for T 
beams. In the case of I beams the variation was from 900 to 1600 psi. 
These high stresses did not relate to any of the computed strengths 
such as 5 vi for the tensile strength and 7.5 ae for the modulus 
of rupture corresponding to his test specimens. He attributed this 
discrepancy to the inelastic behavior of concrete prior to cracking. 

In accordance with the design philosophy of the ACI, 


Wyss and Mattock (1971) investigated the torsion-shear interaction 
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on diagonal tension cracking in a test series consisting of twenty full 
scale I section girders stressed eccentrically. The dimensions of 

the test section were identical to those adopted in their earlier 

Study of pure torsion (1969). In the study of shear-torsion inter- 
action, two cases, web cracking, and shear-flexure cracking were 
examined. Depending on the cross sectional dimensions and the eccen- 
tricity of prestress, initial cracking may occur either on the top 

or in the web when torque alone acts. Their study concerns the former 
case, while the latter case would prove to be more severe in combined 
loading, as the torsional and flexural shear stresses add on one verti- 
cal face. To establish the torque-shear interaction, stress. distribu- 
tion due to prestress, flexure and transverse shear was evaluated by 
the elastic theory and combined with the torsional shear stress com- 
puted by either elastic or plastic theory. Assuming the tensile 
Strength as 6 Agee tensile stress failure criteria was applied. 

Wyss et al noted that the computation of torsional stresses by the 
elastic theory was unrealistic and favoured the plastic theory. Three 
possible modes of diagonal tension cracking were recognized. In the 
case of cracking at the top, its center was examined. In the web the 
neutral axis was considered to be associated with the maximum principal 
tensile stress. However at a cross-section higher tensile stresses may 
occur at the bottom junction with the flange if the rate of increase of 
tensile stresses due to flexural moment is more than that of prestress 
assuming both torsional and transverse shear stresses do not vary in 
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In the development of their shear-flexure cracking strength 


Wyss and Mattock started with the ACI shear equation (1963). 


Wao ee 0b b dy fet ot + Vy (2.19) 
jae 3 
where Lee = diagonal tension cracking shear in a region of a 
beam cracked in flexure 
b' = minimum width of a flanged member 
d = effective depth of a beam 


M/V = moment shear ratio due to applied loads 


Vq = dead load shear 


= 
I 


net flexural cracking moment 


The tensile stresses due to the second and third terms of the fore- 
mentioned equation were combined with the plastic torsional shear 
stress (only in this case tensile strength was assumed equal to 
modulus of rupture estimated at 9 VRE The first term of equation 2.19 
was assumed to be affected by the presence of torque in the same way 
as the shear, torque interaction in the web cracking case. Their 
strength equations were not derived in general terms but certain 
known values were substituted and they are as follows: 

Cracking on the top face 


2 


(2520) 


od V 
tS ea Ont (L2G Ft 39-5) 


« Yo notpsr «nt setile pnd hahaa Tanoonth= 5. 
suxal? AE badosto mssd | 

ee ee 

mead gs Yo Nee ave < bo 7! 

ebaot bal iqqn of sub otter ‘sotz Swsmom = ~=V\M il 
aeatiz-bsol bash = 4V 

2 pasta pot asaya LetukstY gen = |W 


_sto? add Yo awrey Uivtad brs biege2 ef} oF sub zazearte, sttenst ont 
xgoke (edotesot attaslq end mtiw bentdmos sisw folteupa benotinem 
ot Tsups bamuees 2aw diaterte eftenes e260 enna nt ylno) 2egts2 | 

or.S noriaves Fo wread t2ext9 aft Ps & te badamites siudoy" Yo ai Coben 

Yi dme2 ond Nt supio}d to sonazang aid yd batoa?ts od 03 bamuees e6w 

sfad? e269 patigeys dew eft of noteosiatar supios , 16arie ont al 
ntedv92 dud ame? (etanen nt baviasbh son suaw enotteups toned aan 
-ewor lo? 26 ove veld bre eetleitedue stew 2euTay mot : 

ont aps edt no pitdoew | 
1. ia 


ae 


ah 


| git f 
(OS.3) | . Wg ag +°85, 1 TSE a 


35 


Web Cracking 


fits, Se4025ne08-95 View (2s21,) 


Shear flexure cracking 


rd 
e it ii 
Ves = 5.88 (1 = 79>) + (29.52 = 20,600) am | Oe, (2«22) 
where i = diagonal tension cracking torque in combined loading 
T = torque applied to a section 
V = shear applied to a section 
vee = diagonal tension cracking shear in combined bending 


EE Ove diagonal tension cracking shear in a region of a 


beam cracked in flexure 


The predicted results were generally in good accord with the test re- 
Sults except in the case of some girders which had initial cracks in 
the web. They believed that the highly conservative predictions in 
the case of web cracking were due to either; 

(i) the less stressed side of the web reinforcing the highly 
stressed side of the web delaying cracks, or 

(ii) redistribution of stresses occurring which could not be 

taken into account by the elastic theory. 
Probably the high stress gradient also might be responsible for this 


conservatism. 
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2.2 Ultimate Strength 
2.2.1 General 

The general methods developed in the literature for predicting 
the ultimate strength of reinforced concrete beams under pure torsion 
and combined torsion, bending and shear are indicated in the following 
sections. Reinforced concrete flanged beams analyzed by the skew bend- 
ing theory are discussed. Usually the analyses available for reinforced 
concrete beams have been extended with certain modifications to evaluate 
the ultimate strength of prestressed members. The most recent investi- 
gations of prestressed rectangular beams under combined loading and 
prestressed flanged section beams are reviewed closely. 

To achieve post cracking strength and ductility concrete 
members are provided with longitudinal and transverse reinforcement to 
resist torsion occurring alone or with other forces. These reinforced 
beams can be broadly divided into three categories, namely under-rein- 
forced, over-reinforced and partially over-reinforced. In an under- 
reinforced beam both longitudinal steel and stirrups yield before the 
concrete crushes at ultimate. Strain compatibility need not be considered 
in such a beam. Most of the available analyses cover only this category 
of beams. A over-reinforced beam is one in which neither the longitudinal 
steel nor stirrups yield before the concrete fails. Over-reinforced 
beams are normally not permitted to be designed. A partially over- 
reinforced beam has an unbalanced ratio of longitudinal steel to stirrups. 
Only one reinforcement, either the longitudinal steel or the stirrups, 
yields while the other does not as the concrete crushes. For strength 


evaluation of such beams the strain in the steel which is not stressed 
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to yield at failure, must be known. Most of the prestressed concrete 
beams in the literature belong to this category. 

To obtain an under-reinforced or a balanced failure, not only 
must the total amount of longitudinal and web steel not exceed a certain 
amount but also the distribution between the longitudinal steel and stirrups 
must be within certain limits. Though several investigators observed 
that the distribution ratio (ratio volume of longitudinal steel to web 
steel) can vary widely to produce balanced failures, the most often 
used value is unity. Knowledge about the variation of the distribution 
ratio is limited. 

2.2.2 General Theories in Reinforced Concrete 

Early investigators' and the ACI's approach for predicting the 
ultimate pure torsional strength of reinforced rectangular concrete beams 
can be related to an expression of the following form, based on the pre- 
mise that the torsional strength is composed of the sum of strengths 
due to concrete and reinforcement 


A,f 


= weaty. 
Tyerabligweanigys tas (2523) 


Ho 
i} 


where 0 pure torsional strength of a plain concrete beam 
Qeie anfuncttonroft xy and Vy: and ratio of volumes of 


longitudinal steel and stirrups 


ana 
il 


a non dimensional factor 


The first term accounts for the torsional resistance of the concrete and 


second term of the reinforcement. There were considerable differences 
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between estimates of these two torsional resistances cited in the 
literature. 

Rausch (1929), proposing a truss analogy for the study of 
torsion of reinforced concrete beams, obtained a value of 2 for &, as- 
suming uniform yield stress in the stirrups. Concrete strength contri- 
bution was considered to be zero. Anderson (1937) and Cowan (1950) 
disputed the uniform stress in the stirrups and based the stress dis- 
tribution according to the elastic theory. Cowan's study resulted 
in a value of 1.6 for & and the concrete contribution equal to the 
torsional strength of a corresponding plain concrete member. 

The first rational analysis for the ultimate torsional 
strength of reinforced concrete beams under combined loading was that 
of Lessig (1958,1959), who proposed that rectangular beams failed 
in one of two skew bending modes, namely bending mode and torsion mode. 
Using two of the six equilibrium conditions, one for forces and the 
other for moments, together with several assumptions she derived ulti- 
mate strength equations for each of the modes of failure. In her 
analysis the failure surface was idealized by a spiral crack at a 
constant inclination on the three sides of a beam, joined by a com- 
pression crack on the fourth. Inclination of this crack was determined 
by minimizing the torsional resistance of the member. Lessig's theory 
considers the combined resistance of concrete and reinforcement. 
However, it neglects the tensile strength of concrete. 

Those who used Lessig's skew bending theory for reinforced 


concrete members with certain modifications and additions resulting 


{oaer) nswod bas (Veer) noanstinn <b brah a ; 
-ath aasit2 oft boend, bas equiatt2 sit int epeide motion: oft bosuqath 
basfvesy vbute 2'nswod .wragila ets IE 
aig of [supe sottud! stnos atoronoo ad bis. 10% 3.1 Yo aufev 6 ab 
visdnem adevoto> tela ontbneqesen 6 to dipnerte Isnotever’ 
renatevod stembify ady 10t eteyfens: Isaattes suit ont “se 
sett ae odhbeo! bentdnes xobmi enasd etstoqop Sadrdtnter to dignense 
hatts® amped telkdnstoen dend beeoqorq: Onw , (ABeT, BACT) pkezel Yo. 
aban werenad ae ebont Sakbned yYlemsm , zohan’ gnttnied wade ow Yo-sro nt 
afty bas 2zaviot Tot ano anotd tones mu rd ftps ate of? Yo ows pot2u 
_itty boviveb aflz enotighdees: fevevge Wdtw ‘adsgpod 2tnemom vot vento" 
sot aT .switst to 2eohom edd Yo nose to? enofssups dipastte etem 
6 6 dvews fentge & “i bestfeabt caw soetrwe owltet odd 2teyisons: — 
-map 6 vd bantol ,msed 6 to-esbte savla ett no aotdentfont tnstentos 
bontnietet cow aoeto ett to nottentfony .diswe? afd no Joss not2z87q 
vioont 2'gtzes) .sodmem ait Yo sonstetzsy Isnoterod ent entstatnim yd 
tnemestornter brs 9d699n00 ‘Yo aonstztess bentdmo sit evsbtenoa 
aterioios Yo rtenesae sfbenst add ztoofpen Jt .vevewoH! 
besvotntey vo? yrneds pnibasd woda e'ptzzed been orlw pzodT aie 
onttfuzss enotithbs bas weninelvione nitsiie Adtw 2vedmem Siete 


39 


in refinements included Yudin (1962), Gesund et al (1964), McMullen and 
Warwaruk (1967), Collins et al (1968), Goode and Helmy (1968) Brad- 
burn (1968), and Kirk and Lash (1971). Gangarao and Zia (1970) Henry 
and Zia (1971), Johnston and Zia (1971) and Woodhead and McMullen (1972) 
applied her theory in prestressed concrete. Yudin (1962) believed that 
it was accurate enough to assume the failure surface at an angle of 45 
degrees with the longitudinal axis. He claimed it was more rational to 
consider moment equations about two axes, one about the axis parallel 
to the longitudinal axis and the other perpendicular to it; and one force 
equation. Gesund et al (1964a,1964b) considered the effect of dowel 
forces in their study of ultimate torsional strength. McMullen and 
Warwaruk (1967) observed that a beam having more longitudinal steel 
near the bottom face than the top failed developing a compression zone 
in the bottom face, and referred to this as a Mode three failure. 
Collins et al (1968) simplified Lessig's equations and also observed 
the Mode three failure. Goode and Helmy (1968) used a third equilibrium 
equation which furnished a condition to solve another unknown. They 
utilized this condition for estimating the steel stress in a partially 
over-reinforced case. Kirk and Lash (1971) applied the skew bending 
theory to T section beams, idealizing them as rectangular ones. 
Observing that Lessig's theory did not account for certain 
experimental observations and did not numerically fit the test results, 
Hsu (1968d) proposed a theory for pure torsional strength of symmetri- 
cally reinforced rectangular beams based on his test program. Dif- 


fering from that of Lessig he idealized the failure surface as a plane 
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inclined at 45 degrees to the axis of twist but perpendicular to the 
wider face. He stated that the inclination of the failure surface 

was constant irrespective of the cross-sectional dimensions and the 
ratio of longitudinal to web reinforcement. He believed that the 
failure surface should not cut the horizontal legs of the stirrups. 
The dowel forces in the longitudinal steel and the shear resistance 
of the shear-compression zone were taken into account. Establishing 
the equilibrium of moments about the longitudinal axis of the beam, 
the following expression was derived to predict the torsional strength 


of rectangular reinforced concrete beams: 


A Oe By 
2 Lk Lésthhy 
is To + : (2.24) 
where To = 2. Abeh VF é/vb 
2 = 0.66m + 0.33 ¥4/%) 
m = ratio of volumes of longitudinal steel and stirrups 


Hsu evaluated the term 2 experimentally. Thus the equation heavily 
depends on his test results. The first term of this equation is about 
40 percent of pure torsional strength of a plain concrete beam which 
is also equal to the intercept on the Tyo 2X78 of the experimental 
curve between T,,versus reinforcement factor. However, Hsu obtained 
this term as the resistance of the shear compression zone. 

An extensive continuing investigation of reinforced T and 
L section beams under combined loading has been reported by Ferguson 
and Farmer, Ersoy, Victor, Liao, Behara, and Rajagopalan (1967, 1967 
and 1968, 1968a and 1968b, 1969, 1970a, 1970b, 1971 and 1972). The 
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4) 
equilibrium approach as proposed by Lessig (1959) for rectangular 
beams was considered complex for a flanged beam and they analyzed their 
beams by interaction curves and surfaces. 

Another intensive investigation by Mattock et al (1967, 1969) 
concerning rectangular and flanged beams was based on the interaction 
of torsion with other forces. 

In a review paper based on the tests in the literature Zia 
(1970) indicated that I, T, and L sections also fail by skew bending 
and because of the complexity of the problem, one has to depend on 
the approximate approach using the summation procedure of component 
rectangles for computing the strength. Evans, Kemp and Wilhelm (1970) 
stated that the behavior and ultimate strength of such members depended 
on the degree to which the component rectangles act together. Based 
on their investigation referred to earlier, they proposed the following 
equation for predicting the ultimate strength of reinforced concrete 
T and L beams in which there were equal volumes of longitudinal and 
stirrup steel, were component rectangles reinforced transversely 
by closed stirrups and longitudinally at least at the corners, and the 
reinforcement from the component rectangles sufficiently overlapped 


to enable the flanged section to act as a unit. 


Atty 
Tage 0.45 Se 2.02 X13 : (2.25) 
where Us = cracking torque of a plain concrete beam. 


The flanged sections tested by Hsu (1968b) at the PCA did 


not have any web reinforcement. 
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2.2.3 Flanged Reinforced Concrete Beams - Skew Bending Theory 
Gvozdev, Lessig, and Rulle (1968) reported the test results 
of reinforced concrete I section beams under combined loading. The I 
beams had a flange width of 8.66 in., an overall depth from 15.75 in. 
to 18.11 in. and web thickness from 1.97 in. to 3.15 in. and were 
11 feet long. The variables included the arrangement and amount of 
longitudinal and web reinforcements, cross-sectional dimensions and 
T/M ratios. The failure of under-reinforced I beams tested under 
combined torsion and bencing was observed to be similar to that of 
rectangular beams. Two modes of failure namely bending mode and torsion 
mode were identified, which are shown in Fig. 2.1. The bending mode 
of failure was characterized by the development of a spiral crack on 
the bottom and the side faces, with the concrete failing in the com- 
pression flange along the line joining the ends of the spiral crack. 
In the torsion mode the failure crack was formed on the side, top and 
bottom faces. Splitting of concrete on the side faces of the top and 
bottom flanges took place, followed by the compression distress on the 
side of the web opposite to that on which the tensile cracks opened. 
The compression distress in the web was along a line connecting the 
ends of concrete split in the flanges. Gvozdev et al found from their 
tests that the crack inclination of the failure surface to the longi- 
tudinal axis was independent of torque/moment ratio and it followed 
the direction of principal tensile stresses due to torsion only, unlike 
in the case of rectangular beams. Hence they assumed a 45 degree 
inclined crack on all the three faces in the development of strength 


formulas irrespective of the relative proportions of the torque and 
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bending moment. Since the depth of the compression zone was small, it 
was taken equal to zero. The strength contribution of only those stirrups 
crossed by the tensile crack in the web portion was taken into account. 
With the foregoing assumptions, the same skew bending analysis (Lessig 
1958, 1959) developed for the rectangular beams was applied to I section 
beams. The analytical strengths correlated well with the test results. 
However the strength prediction was slightly unconservative for beams 
tested in pure torsion. This is due mainly to the fact that Gvozdev et al 
chose to locate the skew bending axis in the torsion mode, close to 
one of the side faces of the flanges increasing the lever arm of the 
steel forces; though they observed compression distress in the web 
also. Gvozdev et al stated that the post cracking strength in pure 
torsion did not exceed 30 percent even for the highly reinforced I 
beams they tested and the torsional stiffness reduced sharply after 
cracking. They observed that the ultimate capacity of an open cross- 
section member such as an I beam, under torsion was limited by its 
deformation rather than by its strength. 

Thirty seven reinforced concrete T beams were tested by 
Kirk and Lash (1971). Under combined bending and torsion all the 
beams were provided with a symmetric arrangement of top and bottom 
longitudinal reinforcement and most beams had web reinforcement. Al] 
beams were under-reinforced both longitudinally and transversely when 
torsion alone acted. The beams failed in a skew bending mode, bending 
taking place about an inclined axis on the flange connecting the spiral 
tensile crack on the remaining sides. The inclination of the failure 


surface varied from 45 degrees in the case of pure torsion to 90 degrees 
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in the pure bending. Considerable ductility was displayed by the rein- 
forced beams at the ultimate stage. The compression zone in the flange 
was assumed to extend only in the web portion, neglecting the outstanding 
flanges. In essence it was treated as a rectangular beam for analysis. 
The ultimate strengths were computed by two methods, one proposed by 
Lessig (1959), and the other proposed by Collins (1968) for Mode 1 
failure in the rectangular beams. Both methods gave good correlation 
with the test results, but Lessig's method correlated better. 

In an extension of this investigation, Kirk and Loveland (1972) 
carried out tests on 18 T beams with an unsymmetrical arrangement of 
the longitudinal steel in the top and bottom of the T beams and all the 
beams had web reinforcement. At ultimate either the top or the bottom 
longitudinal steel yielded. The stirrups either yielded or approached 
yielding. The location of the compression zone depended on the re- 
lative amounts of the steel in the top and bottom and T/M ratio. The 
compression zone always formed in the flange for the beams with steel 
more at the top than the bottom, this is a Mode 1 failure. In the 
case of beams having more bottom steel, the compression zone formed in 
the flange for low torque to moment ratios, and in the bottom across 
the web for high torque to moment ratios, which is usually called a 
Mode 3 failure. The T section was assumed as a rectangular section 
and a Mode 3 Analysis similar to those of Lessig and Collins was pro- 
posed. For this series also, Lessig's method correlated better. 

The formation of a compression zone on the side did not 
occur in T-beams reported in their investigation; as the transverse 


stiffness with respect to an axis perpendicular to the flange was much 
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larger than that parallel to the flange. In Mode 1, the over hanging 
portion of the flange was neglected in calculating the depth of com- 
pression zone. This assumption is not necessary as the failure loads 
are practically the same, whether the outstanding flanges are excluded 
or included for T beams tested in their investigation (Rao and Warwaruk, 
V973)u 
2.2.4 Prestressed Rectangular Beams Under Combined Loading 

As a first part of a continuing research program dealing with 
prestressed concrete members under combined loading at the University 
of Alberta, Mukherjee and Warwaruk (1970) reported the tests on fifty 
four rectangular (6 x 12 in) prestressed beams. One half of the beams 
were prestressed eccentrically and the other one half concentrically. 
Barring two, each beam had the same amounts of non-prestressed longitudinal 
and web reinforcement. Twenty-eight of the beams were tested under 
differing combinations of bending and torsion, and the remainder 
additionally had shear. 

The inclination and the location of origination of the cracks 
depended on the level and eccentricity of prestress as well as T/M 
ratio. Initial cracking varied from the bottom side for the beams with 
high bending to the centre of side face for beams with high torque. In 
beams loaded under a medium range of T/M ratio cracks started simul- 
taneously at the bottom and side faces. Eccentrically stressed beams 
cracked first at the top of vertical side. Presence of shear caused 
crowding of the flexural cracks in the vicinity of transverse load and 
the crowding of torsional cracks away from the bending load. 


The torsional strength of the beams tested in their investigation 
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increased when bending was present in moderate amounts. This increase was 
greater in the case of eccentrically prestressed beams. Similarly the 
bending capacity was not increased with small amounts of ois but 
was not reduced significantly until the coexisting torque was about 40 
percent of the pure torsional strength. The addition of flexural shear 
reduced the torsional strength regardless of the combination of torque 
and bending moment. They reported that the initial torsional stiffness 
was unaffected by the level of prestress, T/M ratio and the presence 
of shear. The rotational capacity reduced as the level of prestress 
increased and T/M ratio decreased. The reduction in torsional stiff- 
ness was more gradual when bending was predominant. 

Mukherjee and Warwaruk indicated that most of the beams failed 
either in Mode 1 or Mode 2. However eccentrically stressed beams loaded 
under high T/M ratios showed upward deflections at failure. These 
beams may have failed in Mode 3 with the compression hinge in the bottom 
face. When transverse shear was present, as the T/M ratio increased, 
the failure zone shifted away from the bending load. The location of 
the failure zone could not be satisfactorily related, to any known 
parameters for a given T/M ratio. The failure surface is generally 
near the bending load in the bending mode (Mode 1), in the minimum 
moment region for Mode 3, and at any location in the pure torsion 
case. For Mode 2 it may partly depend on T/M ratio. However, it is 
to be noted that in this mode the influence of moment is not consider- 
able. For correct computation of flexural moment, when it varies in 
the test zone, the knowledge of the location of the failure zone is 


necessary. No satisfactory solution is yet available to this end. 
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Mukherjee and Warwaruk observed that the interaction be- 
tween bending and torsion was similar at the cracking and ultimate 
stages, and proposed a parabolic interaction curve involving the 
level of prestress and the eccentricity. This curve reflects the 
beneficial effect of increase in the torsional strength in the pre- 
sence of moderate bending moments. The proposed interaction curve 
showed good correlation with their test results and those of Gangarao 
and Zia (1970) whose beams contained various amounts of web reinforcement. 
Further, Wyss and Mattock (1971) showed that this torque-moment inter- 
action, though derived for rectangular sections, exhibited good agree- 
ment with their prestressed I girders tested in combined bending and 
torsion. 

The effect of simultaneous action of torsion, bending and 
shear can be expressed in terms of an interaction surface. Mukherjee 
and Warwaruk proposed such a surface for the prestressed rectangular 
beams. On the torque-moment plane, the trace of the interaction sur- 
face was the fore-mentioned relationship proposed by them. For torsion- 
shear interaction a straight line relationship was utilized. In the 
past for reinforced concrete beams circular interaction was often 
proposed at cracking based on test results. As the web reinforcement 
amount increases this circular curve flattens and tends to a straight 
line as evidenced by the tests of Klus (1968). Thus a straight line 
is a lower bound in reinforced concrete beams. Because of the paucity 
of test results it is safe to assume straight line interaction relation- 
ship between torque and shear. Assuming that full flexural capacity 


would be reached without a premature flexural shear failure due to 
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enough web reinforcement a rectangular trace was used on the bending-shear 
plan2. This interaction surface predicted very well the ultimate strengths 
of their beams. Subsequently it has also been confirmed in tests of 
prestressed rectangular beams under combined loading by Woodhead and 
McMullen (1972) who reported very good correlation. 

The study of prestressed beams under combined loading by 
Gangarao and Zia (1970), Henry and Zia (1971), and Johnston and Zia 
(1971) was introduced while discussing the cracking strength. Here 
a brief review of their ultimate strength analyses follows. 

Gangarao and Zia observed from their tests that rectangular 
prestressed beams (6"x12") failed in either one of the two skew bending 
modes proposed by Lessig (1959), for reinforced concrete beams and de- 
rived prediction equations for ultimate strength. A bending mode was 
identified by the formation of an inclined compressive zone adjacent 
to the top face of the beam joining the two ends of the spiral crack 
on the remaining sides and by the yielding of either the prestressing 
strands or the horizontal legs of the stirrups near the tension face. 

In the torsion mode the compression hinge formed on the vertical 

side of a beam. This mode was considered to occur when the transverse 
reinforcement near the tension face just yielded. In both the modes 
the strain across the failure surface was assumed to be proportional 
to the distance from the neutral axis situated in a plane having the 
compression zone. This is an approximation only. The spiral crack 
inclination at failure, which is related to the inclination of the 
compression hinge, was based on cracking analysis. At ultimate, if 


either the stirrups or the longitudinal steel yielded and the other 
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did not, the strain in the steel that did not yield was obtained by 
Gangarao and Zia using the following assumption. The strain components 
parallel to the failure surface of both longitudinal and stirrup steel 
were equated assuming that the relative displacement along the direction 
of the failure surface was negligibly small for small rotations. How- 
ever, the rotations at ultimate in general are not small. Strains in 
other steels in the beam were then assumed to be compatible with the 
yield strain of the steel which yielded. However, yield strain was 
assumed in the longitudinal mild steel on the tension face regardless 
of their assumption of compatibility. 

Eccentrically stressed members developed slightly higher 
ultimate pure torsional strength than the concentrically stressed ones. 
Eccentricity of prestress increased the ultimate strength of a member 
under combined loading with a decrease in the T/M ratio. The addition 
of longitudinal mild steel did not result in a significant increase 
in the pure torsional strength, but it increased the ultimate strength 
of a beam under combined loading as T/M decreased. This is partly be- 
cause torsional strength is usually limited by the transverse steel, 
while the flexural strength is limited by the longitudinal steels. Square 
and circular non-dimensional interaction relationships between torsion 
and bending were proposed conservatively, for beams with and without 
longitudinal mild steel respectively. Their theoretical predictions 
were consistent and in good agreement with the test results for high 
and low T/M ratios, but less reliable in the transition from the 
torsion mode to the bending mode which was characterized by the 


yielding of horizontal legs of stirrups in the bottom face. 
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For lack of sufficient data, limits for the amounts of steel 
were suggested only in the case of pure torsion to identify a balanced 
failure. The total volume of steel was not to exceed 3000 VF / fry: 
Further, to ensure that both the longitudinal and the web steels yield 
the total steel needs to be distributed properly between the former and 
the latter. This ratio m, given by the following equation, was suggested 


fo De Ie 0 to) 2s 


m = (A refiny*Ans* py)S/Ag fey (x44¥ )eot@e (2.26) 
For beams under pure torsion, the torsional stiffness was 
unaffected by the spacing of the stirrups, but the ductility was much 
lower for beams with stirrup spacing at 4 in. than at 3 in. Eccentricity 
considerably increased the ductility. Presence of longitudinal mild 
steel resulted in little effect in the deformation prior to cracking, 
in great increase in ductility, negligible increase in pure torsional 
strength but considerable strength increase as T/M ratio decreased. 
The bending stiffness prior to cracking decreased with a decrease in 
T/M ratio, but the torsional stiffness changed only insignificantly. 
It was observed that for Mode 2 stirrups closer to the compression 
zone nearly yielded which indicated that the depth of the compression 
zone was negligibly small. 
The study of Henry and Zia included the effect of shear. 
They observed that all the beams (6"x12") failed in either bending mode 
or torsion mode. The first cracks to form were generally flexural near 


the bending load, though a beam subsequently failed in a torsion mode. 
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Stirrups yielded when intersected by the failure surface at the ulti- 
mate stage. In the torsion mode the strands did not yield while in 
the bending mode they did. Yielding of horizontal legs of stirrups 
observed by Gangarao and Zia did not occur in these beams in bending 
mode failure. Contrary to the finding of Gangarao and Zia, Henry and 
Zia observed that the prestressing strands registered strain prior to 
cracking when T/M ratio was low. The strains recorded in the longitudinal 
mild steel were compressive, but they were assumed to have yielded in 
the analysis. The failure surface established by the cracks was less 
distinct in the bending mode than in the torsion mode. The location 
of the failure surface which is essential for computing the bending 
moment was found to be somewhat arbitrary. 

Analytical expressions based on skew bending theory similar 
to those of Gangarao and Zia were presented to predict the ultimate 
strengths. The agreement between the test and theoretical results was 
good. The following conservative interaction relationships were pro- 


posed between torsion and bending and torsion and shear. 
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It was observed that the ultimate torsional strength of a beam could 
be increased up to about 30 percent of pure torsional strength when- 
ever the amount of moment and shear was less than 20 percent. 


The initial torsional stiffness was found to be nearly 
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independent of the loading ratio, prestressing force and its eccentri- 
city. and span of beam. Post cracking behavior depended mainly on the 
loading ratio. There was little reduction in torsional stiffness be- 
tween first flexural cracking and diagonal cracking. Increases in 
transverse shear decreased the torsional ductility of a beam under 
combined loads. 

Johnston and Zia (1971) observed from their test results 
that prestressed hollow beams developed the skew bending mechanism 
of ultimate failure similar to solid reinforced and prestressed 
rectangular sections. Three modes of failure corresponding to a com- 
pression zone on either the top, vertical side or bottom face occurred. 
When transverse shear was high another mode having a compression zone 
in the top corner was identified, but was analyzed as mode 2. Theo- 
retical analyses based on skew bending theory similar to those of Henry and 
Zia were presented. In mode 1, the failure was considered to occur 
when either the vertical legs of the stirrups yielded or, the longi- 
tudinal prestressing strands yielded. For mode 2 it was assumed that 
failure would occur when the horizontal stirrup legs first reached the 
yield stress. Generally but not always the cracking mode and the 
ultimate mode were the same. Theoretical predictions compared wel] 
with the test results. 

The torsional stiffness was found to be independent of trans- 
verse shear, moment and web reinforcement, but increased with the 
concrete strength in the precracking stage. Post cracking torsional 
and flexural stiffnesses depended on the loading ratios and the amount 


of reinforcement. The necessity of more tests on beams with low torsion 
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and high shear was suggested. 

Woodhead and McMullen (1972) carried out an investigation into 
the strength of thirty rectangular (6"x12") prestressed beams loaded in 
torsion, bending and shear. Two beams had concentric prestress and 
the rest eccentric prestress. Except for nine beams all had web rein- 
forcement. Three of the web reinforced beams had longitudinal mild 
steel reinforcement. They presented analyses for cracking and ulti- 
mate strengths; the former analysis based on the elastic theory and 
the latter on the skew bending theory. The cracking analysis was dis- 
cussed earlier. Ultimate strength analysis was proposed for each of 
the two modes, namely the bending mode and torsion mode observed in their 
tests and also for tests reported in the literature. Analysis of 
each mode covered prestressed beams with or without web reinforcement. 

(a) Bending Mode: A pure flexural failure occurs with the for- 
mation of a rectangular compression zone adjacent to the compression 
side of the beam. With the addition of increasing torsion, this zone 
may become trapezoidal, triangular and may even extend partly on to 
the side before the torsion increases sufficiently to transform this 
mixed bending into a torsion mode. Though this was recognized by 
earlier investigators a rectangular compression zone was used in 
analyses, presumably for simplicity. Woodhead and McMullen took into 
account the changing shape of the compression zone by assuming that 
the neutral axis met one of the vertical faces at an angle which is 


given by the following expression: 


M 
a = tan”! = (2.29) 
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They stated that this refinement resulted in little difference in the 
ultimate strength. Their equilibrium analysis was carried out similar 
to that for over reinforced beams in pure bending using the following 
empirical ultimate strain in concrete both for pure bending and combined 


loading as suggested by Evans and Khalil (1970). 


My 2 
alec Py Ce (2.30) 
uo 
where See longitudinal compressive strain under combined loading 
aah 0.0034 
ste = ultimate bending capacity of beam in pure bending 


If the beams were web reinforced the resistance of web steel was taken 
into account in establishing equilibrium. The strain in the longitudinal 
steel was known and the strain in the stirrups was obtained assuming 
the strain compatibility suggested by Gangarao and Zia (1970). Wood- 
head and McMullen showed good agreement of their bending mode analysis 
with several rectangular and flanged beams reported in the literature. 
(b) Torsion Mode: Woodhead and McMullen applied their cracking 
analysis to compute the ultimate strength of prestressed beams with- 
out web reinforcement failing in torsion. For such beams cracking was 
regarded as constituting failure. The critical location for cracking 
was assumed at the centre of the vertical face for rectangular and I 
beams and centre of the top face for T beams regardless of eccentricity 


and T/M ratio. For web cracking in I beams their theory grossly under 
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estimated the test results as in the case of I beams of Wyss and Mattock 


(1971) Li = 1.62). 


test)/"u(theory) 
For web reinforced beams Woodhead and McMullen considered two 

cases of torsion failure. The first case related to the beams longitudinally 

and transversely under reinforced and the analysis was similar to that pro- 

posed for reinforced concrete beams (Lessig 1959, McMullen and Warwaruk 


1967, Collins et al 1968). Failure was assumed to occur with the yielding 


of longitudinal and transverse reinforcement away from the inclined com- 


pression face. The induced strains due to the loading were neglected 
in the reinforcement near the compression zone. The second case con- 
cerned beams transversely under-reinforced but longitudinally over- 
reinforced and failure was considered to take place when the vertical 
legs of the stirrups attained yield strain. The induced strains in 
the other reinforcements were related to the stirrup yield strain by 

a certain strain compatibility assumption. In both cases the two un- 
knowns, the depth of compression zone and the torque, were obtained by 
solving the two equations established considering the equilibrium of 
forces and moments in a skewed plane. Prestressed beams of several 
other investigators were analyzed by both cases of the torsion mode 
and it was found that only a few beams were under-reinforced; many 
were under-reinforced only transversely but not longitudinally, some 
were over-reinforced. Empirical limits were proposed for the amounts 
of reinforcements to define what constituted under and over reinforced 
beams. Despite a variety of beams from different investigations, their 


strength predictions were in good agreement with test results. 
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Woodhead and McMullen extended the analysis of torsion mode 
to I beams with simplifying assumptions. They computed the ultimate 
strength assuming the compression zone adjacent to the side faces of 
the flanges. The stirrups on the tension side were assumed to yield. 
The induced strain in the strands was arbitrarily taken equal to the 
yield strain of the web reinforcement. It was implied that the force 
in tne leg of the stirrups adjacent to the compression zone was neglected, 
thus the strength was partly underestimated. This may be justified if 
the compression zone is very near to this vertical leg, but for I 
beams it is not. It was not indicated how the inclination of the com- 
pression hinge was arrived at. They tested the validity of the theory 
by applying it only to the test results of Gausel (1970). 

In addition to several other conclusions, it was stated that 
the shear force reduced the ultimate strength in the torsion mode of 
failure. The theory was conservative near the torsion axis of the 
interaction diagrams as small amounts of bending moment increased the 
torsional strength above the pure torsional strength. The bending strength 
was equal to the ultimate flexural capacity until the coexisting torque 
was of the order of 30 percent of pure torsional strength. Circular 
interaction was conservative in both the torque-moment and torque- 
shear planes. The three dimensional interaction surface proposed by 
Mukherjee and Warwaruk (1970) predicted well the ultimate strengths 
of their beams. 

2.2.5 Prestressed Flanged Beams 

Gardner (1960) described the testing of 16 identical plain 

I beams, prestressed at a small eccentricity of 0.2 in., under combined 


bending and torsion. The beams were loaded initially with a predeter- 
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mined bending moment which varied from 10 to 80 percent of the pure bend- 
ing capacity of the section, and then twisted to failure. None were 
tested under pure torsion. Within the range of flexural moments ap- 
plied he observed nearly the same ultimate torsional strength for 

all the beams. The combined effect of prestress and the normal bend- 

ing stress resulted in a moment with eccentric prestress before torsion 


was applied. 


Gardner's beams seem to have failed in Mode 2, triggered by 
cracking in the web. Some I beams because of eccentric prestress might 
have developed initial cracking on the top face, but failed in Mode 2 
as the stiffness about an axis parallel to the web was less than that 
about an axis perpendicular to the web. Mode 2 failure most probably 
is the reason for the uniform failure loads in all the beams. No new 
theory emerged from study of I beams but he computed the ultimate torque 
based on the plastic theory which was close to the ultimate torsional 
strengths. 

Reeves (1962) described 42 tests carried out on three series 
of plain prestressed concrete T-beams subjected to combined bending 
and torsion. Beams within each series had identical properties and 
between the series only the breadth of the flange was varied. Pre- 
stress distribution in all the beams was triangular. Load application 
was similar to that used by Gardner (1960). The application of a pre- 
determined bending moment caused in some beams initial cracking before 
the torsion was added. 


The failure was abrupt and occurred by the development of a 
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tensile crack on one side of the web which passed through the section 
resulting in crushing on the opposite side of the web if the coexisting 
bending moment was between 0 to 70 percent of the pure bending capa- 
city of the section. For larger moments a combination of bending and 
twisting failure occurred. 

His beams exhibited considerable increase in torsional strength 
up to a certain increase in the bending moment. The optimum increase 
in the torsional strength, as much as 1.65 times the pure torque capa- 
city depending on the geometry of the section, occurred when the 
flexural moment was of such magnitude as to cause a direct tensile 
stress in the bottom face equal to about the split cylinder strenath. 
The torsional strength did not fall below the pure torque capacity 
until a bending moment equal to 80 percent of the pure flexural 
capacity was applied. Any increase in the moment above this value 
caused a rapid reduction in the torsional strength. 

Reeves proposed a cubic polynomial interaction equation 
for each of his series based on a regression analysis of his test re- 
sults. Other typical interaction curves available suggest that the 
torsion-bending interaction depends on geometry of section, distri- 
bution of prestress, and the reinforcement details. Reeve's inter- 
action considers only geometry of the section. For estimating the 
ultimate strength of plain prestressed T beams under moment and tor- 
sion use of interaction curves was suggested. Reeve's compared three 
methods to obtain pure torsional strength. 

(1) Plastic theory, with the direct stress at the critical lo- 


cation, using the failure criteria due to Brester and Pister (1958). 
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(2) Plastic theory, with the direct stress at the critical lo- 
cation using principal tensile stress theory and the tensile strength 
equal to split cylinder strength. 

(3) Elastic theory with the direct stress at the critical lo- 
cation, using principal tensile stress theory and the modulus of 
rupture. He recommended the third method as he found it statistically 
best. 

Bishara's experimental investigation (1969) consisted of 
tests on prestressed concrete rectangular, I and T section beams with 
web and longitudinal mild steel reinforcement, under combined torsion 
bending and shear. For each type of cross section eight beams were 
tested. One half were concentrically stressed and the other half 
eccentrically stressed. These beams were similar to Zia's (1961). 

The load deformation curves were observed to be linear up 
to about 90 percent of the cracking loads. Tors‘onal stiffness as 
estimated by the elastic theory was close to the average measured 
values in the case of I and T beams but for the rectangular beams the 
correlation was poor. It was found that the initial torsional stiff- 
ness was independent of the T/M ratio and the eccentricity. His 
description of failure modes implied that the flanged beams generally 
failed in skew bending. On the average the ultimate load was higher 
than the cracking load only by about 15 percent. 

Bishara found from his test results that the torsional 
strength of prestressed beams under combined loading could be in- 
creased from 1.3 to about 2.5 times the pure torque capacity, irre- 


spective of the eccentricity, by the addition of bending moment. The 
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higher figure corresponded to flanged, members. It was further stated 
by Bishara, that the torsional strength was below the pure torsional 
capacity only if the moment and shear exceeded 80 percent of the strength 


of the beam when bending and shear acted alone. 


Based on his limited test results Bishara presented the 
same form of parabolic interaction diagrams both for the torque- 
moment, and the torque-shear planes for each shape of the cross section. 
In establishing the interaction relationships there were no experimental 
values for pure torque and bending and shear capacities. The small 
span of 30 in. between the loading brackets might have (Pandit, 1970) 
influenced the test results. Woodhead (1969) showed that the torsional 
strength of a specimen significantly increased if the test length be- 
came small. Further the interaction relation will be influenced by 
the eccentricity of prestress; however the proposed curves were in- 
dependent of eccentricity. Also the magnitude of the failure bend- 
ing moment was taken as that relating to a section adjacent to the 
loading bracket in the test zone. Bishara used the following empirical 
expression to compute the pure torsional strength to establish his 


interaction curves 
Tio = 108, rat + 120,) + 1.6 x1y, Ay fey/S (231) 


where ie = ultimate pure torsional strength of a prestressed 


concrete rectangular, I or T section member 
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Vo = volume of plastic membrane on section having unit 


slope 


In the above expression the contribution of stirrups was added to the 
torsional resistance of plain prestressed concrete beam estimated by 
the plastic theory. The correlation between Zia's (1961) pure torsion 
test results and those computed by this equation was shown to be good. 
However, Henry and Zia (1971) showed that it grossly overestimated 


the pure torque capacity of Gangarao's (1970) beams. 


Zia (1961) found that the ultimate strength of reinforced 
and prestressed concrete rectangular and T sections without web rein- 
forcerent under torsion was the same as the cracking strength, com- 
puted using the elastic theory of torsion, and the modified Cowan 
failure theory. However, for I section beams the ultimate strength 
was substantially different from the cracking strength; this was 
attributed to the redistribution of stress and transformation of the 
load resisting mechanism, discussed already. The ultimate strength 
of reinforced and prestressed members with web reinforcement was taken 
as the strength of that member without web reinforcement, to which 
the strength contribution of web reinforcement was added. The tor- 
sional strength of web reinforcement was estimated using the proposal 
of Cowan, referred to earlier in this section. In addition Zia tested 
three types of sections, which were non-prestressed. For such beams, 
the presence of web steel increased only the ductility but not the 


strength. He reasoned that this was due to the location of longi- 
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tudinal bars at the corners of narrow stirrups where stresses were low. 
Thus there was no longitudinal steel at highly stressed regions to re- 
sist the horizontal component of diagonal tension. His argument may 
not be completely valid as cracking terminates the elastic distribution 
of stress. However he reported that in prestressed beams, except for 
one series of T beams, the stirrups yielded resulting in their full 
contribution to the strength. He explained that though in prestressed 
beams the longitudinal bars were at the corners of the stirrups as in 
the non-prestressed specimens, the horizontal component of diagonal 
tension was resisted by the strands. In Zia's computation of ultimate 
strength of web reinforced specimens conrete contribution is equal to 
the cracking strength of the plain specimen, however, in the ACI method 
based on a similar hypothesis this contribution is only 40 percent of 
the cracking strength. His beams had web reinforcement equal to about 
1.5 percent of the volume of care enclosed within the stirrup. It was 
observed that the web reinforcement did not influence the elastic 
torsional behavior of the member, but enabled the specimen to fail 
gradually developing a good amount of ductility. Spacing of stirrups 
was recommended at not more than 0.4 times the depth of the member. 
Gausel (1970) described tests on 16 eccentrically prestressed 
I section beams. One half of each beam was provided with web rein- 
forcement and the other half was not. Thus 32 tests were performed. 
The beams were subjected to differing transverse loads and then gradually 
increasing torque. Eight beams were tested under torsion and bend- 
ing and the remainder had additionally shear. Only one beam was tested 


in pure torsion. In several of these beams, the amount of web rein- 
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forcement used was very small. In some beams the one half portions of 
the beam having web steel failed at a lower torque than the other half 
without web steel for which Gausel could not give any explanation. 
However, he observed that the web reinforcement used had no effect on 

the torsional strength of beams tested under torsion, bending and 

Shear. Even a small amount of web reinforcement resulted in increases in 
ductility and prevented an explosive failure. It was noted that the 
Sequence of loading did not have any effect on the ultimate strength, 

an observation made earlier by Pandit and Warwaruk (1965). 

Beams were divided into three groups for the purpose of 
analysis. The first group related to beams which developed little or 
no cracking under pure torsion or torsion and bending. The ultimate 
strength of these beams was calculated as the sum of the cracking 
strength and the contribution from stirrups; the former being esti- 


mated by the plastic theory and the latter by Hsu's equation. 


ar I 4 T Stirrup (2.32) 
where Pee ae eo a 
Tee erin = (0.66m + 0.33 y_/%y) A, XYq ays 
m =] 
VY = volume of sand heap with unit slope 
fi, = tensile strength taken equal to split tensile strength 


it 


Though the prestress was eccentric, Op was taken as the average pre- 
stress. In I beams Zia (1961) and Gausel considered the full contri- 


bution of concrete strength, while Wyss et al (1969) took into account 
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only a fraction of it (about 40 percent) and increased the lever arm 
of the stirrups, in estimating the ultimate strength. Gausel assumed 
that stirrups contributed only to the torsional strength but: not to 
other strengths. 

The second group consisted of beams under torsion, bending, 
and shear with minor cracking. Gausel found that the principal stress 
criteria gave poor correlation with test values and the following 


straight line interaction was suggested. 


EE A yg ree Eley 


The third group consisted of beams subjected to torsion 
and bending or torsion, bending and shear with large flexure-shear 
cracking; the analysis was based on the shear failure theory of Walther, 
assuming uniform distribution of normal, and flexural stresses over the 
compression zone, and torsional shear stress over the whole cross- 
section. For torsion-bending a nearly rectangular interaction diagram 
and for torsion, bending and shear a circular interaction diagram 
were proposed. The rectangular interaction relationship was checked 
by Gardner's (1960) unstirruped beams. 

In a study of torsional strength of eccentrically prestressed 
concrete bridge girders, Wyss, Garland and Mattock (1969) tested 
eighteen full scale I section girders and six rectangular section 
members with web reinforcement. The principal variables were level 
of prestress (average prestress 0, 530, 1150 psi) and amount of web 


reinforcement. I section girders were divided into three aroups. 
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One group did not have any prestress and the other two had two levels 


of prestress. Within each group the level of prestress was constant 

and the girders were provided with varying amounts of web reinforcements 
ranging from 2 in spacing of stirrups to no stirrups. Rectangular 
members had the same height, width of web, and web reinforcement as 

the I section girders and served the purpose of comparison. For a 

given amount of web steel, longitudinal reinforcement in the form of 
straight deformed bars was designed by the following equation suggested 


by Hsu (1967). 


(x, ty, ) 
Ae ete Neal ba eg ce 
Ay = 2 5 A, Z (2.34) 
ty 
where F = effective prestress 
A, = total area of longitudinal reinforcement to resist 


torsion 


This equation implies equal volumes of longitudinal and transverse 
reinforcements and that the initial prestress in the tendons pro- 
vided the force that would otherwise have been supplied by the longi- 
tudinal torsional reinforcement. By this scheme some girders did not 
require any longitudinal torsional reinforcement. For I sections 
the web reinforcement was made up of two separate vertical legs and 
two horizontal bars across the full width of the flanges. 

They observed that the initial torsional stiffness up to 
torques less than 80 percent of cracking torques was linear. It was 
taken: as the average slope of the linear part of the torque-twist 


curve and was found to be in reasonable agreement with that calcu- 
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lated using elastic theory. They further observed that the initial 
torsional stiffness increased with the level of prestress but it was 
unaffected by the amount of web reinforcement. The eccentricity of 
prestress resulted in progressive cracking and because of this the 
reduction in torsional stiffness at torques slightly above the crack- 
ing torque was not as abrupt as in the case of reinforced girders. 
However further slight increases in torque caused large reductions 

in stiffness. As the jevel of prestress increased, the torsional 
ductility reduced. 

The basic mode of failure for all the I girders was observed 
to be skew bending. They stated that the skew bending axis was pro- 
bably inside the web, close to one of its faces. But the Russian 
investigators Gvozdev, Lessig, and Rulle (1968) proposed that the 
skew bending axis was close to one of the side faces of the flanqes. 
Wyss et al ruled out the possibility of location of the neutral axis 
as suggested by Gvozdev et al. They asserted: 

(a) The critical crack in the web did not penetrate the full 
thickness of the web and was not visible on the compression side 
of the girder as shown in Fig. 2.2. 

(b) Shear compression forces developed near the inside faces 
of the flanges but the compression distress in flanges was never 
oriented as postulated by Gvozdev et al. 

(c) If the neutral axis were located as hypothesized by Gvozdev 
et al, the internal lever arm of the stirrup forces increased very 
much resulting in gross over estimates of strengths. 


Hsu (1967) proposed the following equation for calculating 
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ultimate torsional strength of concentrically prestressed rectangular 


section beams. 


Ai Xa Vat 
ims Tiere pele (2.35) 
uO C S 
ye 100 
where T= 7 eeeav es [2.5 v(1+ —P) - 1.5] 
VX Cc 
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where X width of a rectangular section 


=< 
HI 


depth of a rectangular section 


Hsu's equation (2.35) was found consistently to underestimate 
the strength of eccentrically prestressed I section members. Wyss et 
al modified this equation by increasing the lever arm of stirrup forces 
from x, to (x, +b, )/2 and arrived at the following equation. 

(x4 +b,,) Y) 
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x = width of a rectangular section 
y = depth of a rectangular section 


The ultimate capacity was limited to 


2 
=147ft 54% 
Cc 


To(max) 3 


a © vit 
wees o@ 
ieee | - : ; 


se : s 7 iy 


ge sayy cage #47 


1 a ~avthts 


[é.f- Ce 


aettoee veluntadosy s to dtbtw = x ayoriw 
notiose knee tne & Yo dtgsh = YO 


waste of vitnetatenan bawe? enw (26.8) notseupe 2 ‘wel 
$9 2avW .2radmem notdoee | baeeaytests yf fsotities39 to Azeneite a << 
asovet qurty2 20 me vavel 9A? gataservont yd motvaups etd? bot ribo Is 
_wotteues pnfwotlo? siv ts bavtwre bra S\C dt p>) ot 4x ort 


(98.8) vs £ — + a Am 


at ae 
ot ~ Glog 4118 2.8) 8 a ¢ ent oe 


notise2 setupnedaey 6-70 AibIw = 9x) ok Ee 
aottszae “eTueistoen & to Mqeb= Yo = = 
od bettmel cow wirosqss stant fu oat i 

oll wn 


nh - 
WET ANA taaapuT ni 


o@ 


69 


In deriving the equation (2.36), the prestress was assumed to affect 
only the contribution of concrete to torque, hes As in reinforced 
concrete, the value of te was less than the cracking torque. But 
Chander, Kemp, and Wilhelm (1970) observed that the level of prestress 
influenced the reinforcement factor. However the equation was very 
well substantiated by test results of their girders. 

The minimum reinforcement necessary to develop strengths 
greater than cracking torque, for the girders was obtained using the 


equation shown below derived for rectangular concrete beams. 


A 
Bee eye 
es F tty) x; fal) 
where x = width of a rectangular section 


depth of a rectangular section 


<< 
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It was stated that the presence of prestress or the flanges in these 
beams did not appear to change the minimum reinforcement requirement 
The maximum stirrup spacing was not found to be controlled by the 
web thickness. To limit crack widths it was recommended that stirrup 
spacing should be less than one-half of the flange width or one and 
one-half times the web width, whichever was greater. The form of web 
reinforcement used in their I section girders was found to be adequate 
in lieu of closed rectangular stirrups. 

| A beam of given dimensions can be provided with flexural 
reinforcement as specified by ACI to resist a given moment. A beam 


of identical dimensions can be designed in an under-reinforced way 
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with certain transverse reinforcement and corresponding longitudinal 
torsional reinforcement to resist a particular pure torque. If the 
beam is provided with a total reinforcement corresponding to pure 
moment and pure torque, it is reasonable to expect the beam to fail 
in an under-reinforced way developing simultaneously a bending moment 
equal to pure moment capacity and a torque equal to pure torsional 
capacity. This was possible in reinforced concrete beams (Osburn et 
al 1969). ACI permits this approach and the longitudinal reinforce- 
ment can be designed independently for bending and torsion and added 
in a combined loading case resulting in a nearly rectangular interaction 
relationship. 

In the second part of the investigation consisting of eleven 
girders divided into three series Wyss and Mattock (1971) examined 
the possibility of a rectangular interaction relationship between 
torsion and bending, of eccentrically prestressed I girders. The 
Six specimens of the first series were provided with #4 stirrups at 
9 in. centers and the four specimens of the second series at 6 jin. 
centers. The eleventh girder together with two others constituted 
the third series and contained only longitudinal reinforcement with- 
out transverse and unprestressed reinforcement. Girders within each 
series differed only in the unprestressed longitudinal steel. One 
in each series was tested to obtain the pure torsional strength. One 
girder furnished the pure flexural strength for all the series. The 
remaining girders of the three series were subjected to combined 
torsion and bending in a ratio of approximately that of pure torsional 


strength to pure flexural strength. 
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The girders under combined loading failed in a over-reinforced 
manner without developing either the pure torsional strength or the 
pure flexural strength. Thus the rectangular interaction could not 
be obtained. They observed that the torsionally induced cracking tended 
to split off the overhanging parts of the flange restricting the ef- 
fective compression zone to web only and sufficient compression could 
not be mobilized to produce the yield force in the reinforcement. How- 
ever it may be pointed out here that the prestressed beams of Gangarao 
and Zia when augmented with unprestressed longitudinal reinforcement 
developed rectangular interaction between bending and torsion. The test 
specimens of Wyss and Mattock were concentrated in one region of an 
interaction diagram and were very closely represented by the inter- 
action relationship proposed by Mukherjee and Warwaruk (1970). 
2.2.6 Space Truss Theory 

The behavior of reinforced concrete beams under torsion differs 
considerably before and after cracking. Hsu (1968d, 1968e) and Lampert 
and Thurlimann (1968) showed experimentally that hollow beams developed 
initial cracks at a smaller torque than the corresponding solid beams 
of the same size and concrete strength. This can be easily explained 
by the contribution of the concrete in the core to the cracking torque. 
Their tests also showed that the ultimate torques were the same for 
hollow and solid beams. This is probably due to the fact that at 
failure the concrete core is penetrated by cracks. Lampert and Thurli- 
mann found that in the advanced post cracking stage solid and hollow 
beams had the same deformations and the same stresses. Hence Lampert 


et al concluded that due to the elongation of the stirrups an outer 
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comparatively thin concrete supporting shell was formed which alone 
carried the load i.e. the strength depended only on the rein*orcing 
cage. The remaining core did not contribute materially to the ultimate 
strength. Because of this, in their theoretical study of the space 
truss, equally reinforced solid rectangular sections were considered 
equivalent to fictitious thin walled box sections . 

Based on their tests carried out on beams with square cross- 
section of 20 x 20 in., Lampert and Thurlimann confirmed the general 
validity of simple truss theory of Rausch and improved upon it. Their 
Space truss model consists of longitudinal reinforcement which is assumed 
to be concentrated into stringers at the corners and intermediate shear 
walls. Stirrups act as posts in these shear walls and the concrete 
between the inclined cracks provides the compression diagonals. All 
the reinforcement is assumed to yield at ultimete. The concrete com- 
pression diagonals are necessary for equilibrium only; concrete does 
not provide any additional contribution to the torsional resistance. 
Lampert et al found from tests that the crack angle, 6, depended on 
the ratio of transverse to longitudinal reinforcement amounts. The 
inclination of the inclined cracks and hence the inclination of the 
concrete compression diagonals changed in spite of existing cracks to 
accommodate a redistribution of the forces between the longitudinal 
and transverse reinforcements especially after yielding of one part of 
reinforcement unlike in simple truss theory. It is interesting to note 
here that Hsu's surface of failure does not depend on the amounts of 
reinforcement but is a plane perpendicular to the wider face and in- 


clined at 45° to the axis of the beam. In reinforced concrete investi- 
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gators using the Lessig-type skew bending theory obtained the crack 
angle by minimizing the expression for torque which was a function of 
amount of stirrups and of longitudinal steel. Others working with 
prestressed concrete arrived at this angle from the stress situation at 
cracking. 

Lampert et al used the location of longitudinal bars to deter- 
mine the cross-sectional dimensions of the failure model because they 
observed that the diagonal forces in the shear walls were deflected 
into the adjacent wall by the longitudinal corner bars. The lengths 
of the moment lever arms corresponded to the distances between the 
longitudinal bars in their theory. Other theories assumed the center- 
lines of the stirrups, Xy> Yy> as the lengths of the moment lever arms. 
In fact these lever arms depend on the depth of the compression zone. 

The former assumption is more conservative. However Elfgren (1972), 

and Kuyt (1971) working with the truss model used the latter assumption. 
This difference is significant for the small cross-sections commonly 
used in testing. Rational determination of the cross-sectional dimen- 
sions of the failure model is difficult because of the complex situation 
existing at the corners of the beam. 

For solid sections Lampert et al suggested a wall thickness, t, 
for the corresponding hollow beam to be h/6 or h,/s whichever was smaller. 
In space truss theory the strength of under-reinforced beams is indepen- 
dent of the wall thickness and a nominal wall thickness is needed only 
to check the possibility of an over-reinforced failure. Some investiga- 
tors and codes such as the ACI 318-71 have assumed that the concrete 


contributes to the torsional strength due to the fact that torsion is 
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considered similar to shear. However Leonhardt (1970) has argued that 
in torsion there are no favourable effects such as the shear carrying 
capacity of the compression zone or the flat inclination of the com- 
pression struts. The aggregate interlock does not come into play as 
the torsional cracks opan perpendicular to the crack inclination. 

In this respect torsion is more analogous to flexure rather than shear. 
For flexure no theory is proposed for the independent contribution of 
concrete to the strength of the beam. 

The function of the longitudinal reinforcement is mainly to 
resist the horizontal component of the diagonal tension i.e. to satisfy 
the equilibrium. For this purpose the longitudinal steel need not be 
placed at the corners. The longitudinal steel located in the center of 
cross-section such as prestressing steel is suggested by Lampert et 
al to be effective. However to anchor the concrete diagonals and to 
prevent a push-out of the diagonals (Mitchell et al, 1971) longitudinal 
reinforcement in the corners is stated necessary. But based on certain 
test observations concerning ductility and brittleness of failure McGee 
and Zia (1973) have recommended that the prestressing tendons may not 
be considered effective as longitudinal torsion reinforcement. 

The space truss model for rectangular hollow cross-section 1s 
shown in Fig. 2.3. The following expression was obtained for the pure 


torsional strength of a symmetrically reinforced beam. 
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FIG. 2.3 SPACE TRUSS MODEL FOR PURE TORSION 
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where A = area of one longitudinal mild steel bar 


tay = yield force of longitudinal mild steel bar 


The angle between the concrete diagonal and the beam axis, a, was given 


by 


WIP ON MSS eas nee ys (2539) 


An unsymmetrically reinforced section could be replaced by a symmetri- 
cally reinforced section having the longitudinal steel with the lowest 
yield force, on all sides since failure was governed by yielding of the 
weakest longitudinal bar. 

The space truss model for pure torsion was extended to beams 
subjected to torsion and bending taking into account the differing 
forces in the top and bottom longitudinal reinforcement. Lampert and 
Thurlimann (1969) used the same torsion lever arm, ho» for the pure 
bending capacity of a section instead of the flexural lever arm. How- 
ever recently Collins and Lampert (1971) suggested that the skew bend- 
ing theory developed by Collins (1969) be used for computing the pure 
bending capacity of a section since it contains flexural lever arm, dd. 

Leonhardt (1970) observed that if the longitudinal reinforce- 
ments required individually for torsion and bending were superimposed, 
the resulting interaction diagram between torsion and bending was 
rectangular apart from a small section in the outer corner. 

: Tests aimed at web compression failure showed that the con- 


crete compression struts failed at a lower nominal torsional shear stress 
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than in flexural shear. Lampert and Thurlimann attributed this to the 
warping of the initially plane walls of the hollow box beams into 
hyperbolic paraboloids. This caused eccentric stressing of compression 
struts forcing failure sooner than in flexural shear where the web 
remained plane. Hence in torsion and torsion and shear lower values 
are to be permitted for the upper limit of shear stress than in flexural 
shear. 

Kuyt (1971) showed that the space truss theory and Lessig's 
skew bending theory (Ultimate Equilibrium Method) were not each others 
antipoles but on the contrary would yield the same results if certain 
assumptions were made. He pointed out that the skew bending theory con- 
Sidered the equilibrium in a possible failure surface, provided no infor- 
mation with respect to failure surface or steel stresses in the several 
parts of the reinforcement crossing the failure surface. Kuyt preferred 
the space truss theory because it was on a more definite basis with 
regard to the internal equilibrium of the beam as a whole. He dif- 
ferentiated between beams with longitudinal bars concentrated in the 
corners and beams with these bars distributed along the four sides 
and stated that the assumption that all the steel crossing the failure 
crack was valid only in the latter kind of beams. It was also pointed 
out by him that for the beams with the longitudinal bars at the corners 
only, the assumptions that the inclination of the failure crack was the 
same on each side of the beam and that the stresses in the stirrups 
crossing the failure surface were of the same magnitude, were invalid. 

Based on a truss model for combined torsion, bending and 


shear Kuyt (1972) established equilibrium for each one of the sides of 
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the beam; but this results in the violation of compatibility of the 
failure mechanism. 

Lampert, Luchinger and Thurlimann (1971) extended the space 
truss theory to cover prestressed beams. They carried out tests on 
two 20 in. square, pretensioned beams with stirrups. Of the two, one 
beam was solid, concentrically stressed and tested in pure torsion 
while the other was hollow eccentrically stressed and tested in com- 
bined torsion and bending. They found that there was no basic difference 
between reinforced and prestressed beams with regard to ultimate strength 
and concluded that the space truss theory was applicable to prestressed 
beams if the initial prestress in the tendons was at the usual level 
and the strands yielded after a load-induced stress nearly equal to the 
yield stress of the mild steel. It was shown by them that a prestressed 
and an unprestressed beam of about equal longitudinal yield force sub- 
jected to torsion and bending exhibited the same ultimate moment and 
angle of twist at yield. The strands are to be bonded to the concrete 
for the space truss theory to apply. 

Woodhead and McMullen (1972) found from their test results 
that the truss theory was far more conservative than a design procedure 
based on the ACI code (318-71) requirements for reinforced conrete 
Since the truss theory neglected the concrete contribution to the tor- 
sional strength. 

Lampert (1971) stated that the space truss theory was valid 
for solid and hollow reinforced and prestressed concrete beams of general 
Shape subjected to torsion, bending and axial force, provided the longi- 


tudinal and transverse reinforcement yielded at failure and the St. Venant 
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torsion was predominant. From his study of a general cross-section he 
found that in the collapse mechanism of such a section except the two 
longitudinal bars defining the compression hinge at failure, all] 
transverse and longitudinal reinforcement should yield. 

In a recent paper Lampert and Collins (1972) have contended 
that the problem of torsion and bending is basically solved. They 
suggest that space truss theory furnishes a rational accurate and wel] 
verified model to predict the pure torsional strength. The second step 
necessary for the design of torsion and bending is the interaction re- 
lationship between torsion and bending. It has been shown by them 
that both skew bending and space truss theories are in close agreement 
and predict the same parabolic interaction between the torsional and 
flexural strengths. If torsion is predominant there is a truss-like 
behavior up to failure as in pure torsion. If bending dominates the be- 
havior is similar to that in pure bending except that the compression 
zone is inclined. When space truss theory is extended to combined tor- 
Sion and bending the flexural capacity predicted by the space truss theory 
deviates from the conventional flexural strength of the beam as this 
theory uses the torsional lever arm, bo: in the whole range of torsion- 
bending interaction. Similarly skew bending theory uses the flexural 
lever arm, though in space truss theory it should change to torsional 
lever arm for strength predictions near the torsion axis of the inter- 
action diagram. Lampert and Collins for design suggested equal volume 
of longitudinal and stirrup steel as has been suggested and used by 
others. The longitudinal reinforcements required individually for 


torsion and bending are to be superimposed. Lampert and Collins allowed 
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for a reduction in the torsional longitudinal reinforcement in the 
flexural compression zone due to the presence of flexural compression. 
In the case of continuous beams reduction of longitudinal reinforcement 
as a result of presence of flexural compression may be only slight since 
reinforcement for anchorage in the vicinity of inflection points must 


still be provided and would exist in the compression regions. 
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CHAPTER ITI 
ANALYSIS OF CRACKING STRENGTH 


3.1 Introduction 

This chapter deals with the development of equations for the 
cracking strength of prestressed Concrete I beams subjected to com- 
bined torsion bending and shear. Two methods of cracking analysis are 
presented here. The first is based on the classical theories and the 
second is concerned with the finite element method using three dimen- 
sional hexahedron elements. The inclination of the initial crack ob- 
tained from the cracking strength analysis by Method 1 is used in the 
ultimate strength analysis in defining the symmetric skewed failure 
surface. The experimental cracking strength is defined to correspond 
to that load which causes initial visible cracks on the surface of the 


beam. 


3.2 Method 1 
3.2.1 Assumptions of Analysis 
The following assumptions are made in developing the equations 
for predicting the cracking torque of prestressed concrete I beams. 
1. The torsional shear stresses are distributed on the cross 
section according to plastic theory. 
2. The flexural shear stresses and the flexural normal stresses 


are distributed according to elastic theory. 
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3. Cracking initiates when the maximum principal tensile stress 
equals the tensile strength of concrete. 
4. The transformed area of steel is taken into account in the 


calculation of section properties. 


5. The stresses in the steel due to the applied loads are 
negligible prior to the formation of the initial crack. 

Nylander (1945) was probably the first to recognize that 
the maximum principal tensile stresses computed using the elastic 
theory in the case of beams subjected to torsion, were much dif- 
ferent from the usual tensile strength of the concrete. He proposed 
that the torsional shear stresses be assumed in accordance with the 
plastic theory. This approach has been preferred by Kemp, Wyss and 
Mattock. 

According to assumption 1, the torsional shear stress is 
computed considering full plastification of the section. On the cross 
section of a bar under torsion, as the shear stress reaches the yield 
strength of the material at any point, there the two orthogonal com- 
ponents of the shear stress t, and Ty have to satisfy the condition 


X 
of plasticity and equilibrium (Nadai, 1950). 


oh + ce = xe (3.1) 
where K is a constant 
OT OT 
see Sere Reis (32) 
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83 
Expressing the shear stress components in terms of a stress function 


F(x,y) 
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the condition of plasticity can be given as 


2 2 
OF 
s (ay) 


(2h) 2c 


ay = K (3.3) 
The plastic stress function F(x,y) represents the sloping plastic 
stress area and K is the steepest slope of this area. From Nadai's 
Sand heap analogy, the plastic stress function at any point is pro- 
portional to the height of the sand heap at the corresponding point 
and the uniform plastic shear stress is proportional to the slope of 
the sand heap. The applied plastic torque is proportional to the 
volume of the sand heap. The ideal sand heap covering the plastified 
section of the I beams tested in the present investigation is shown 
in Fig. 3.1. Using unit slope the volume of this sand heap varied 
FEHOM S40t0P S925 in? for various nominally identical beams. 

The stress distribution due to all other forces acting on 
the beam is determined using elastic theory. 
3.2.2 Theories of Failure for Concrete 

Cracking of concrete under combined loadings will occur 
under a complex state of stress. This of course is the case with 


crushing also under such loading. Quite often the stress condition 
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FIG. 3.1 IDEAL SAND HEAP COVERING PLASTIFIED SECTION 
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in a structural member is idealized to a biaxial state of stress 
assuming the stress in the third direction equal to zero. Various 
strength theories have been proposed in the past to determine the 
failure of structural materials. Of these, Mohr's theory o* failure 

is the most appropriate for a brittle material like concrete. Various 
stress conditions in a concrete differential element can be represented 
by Mohr's circles. A failure envelope can then be obtained by joining 
all those Mohr's circles that represent failure conditions. Such a 
failure envelope has to be established for given material based on 
tests under combined stress conditions but this is not easy. Because 
of this, different theoretical shapes of this failure envelope have 
been proposed using simplifying assumptions. 

| The commonly used failure hypothesis in torsion studies are 
Shown in Fig. 3.2. Rankine's principal stress theory assumes that 
failure is attained as soon as the magnitude of the principal stress 
reaches the strength of concrete in uniaxial compression or tension. 
The failure envelope on the Mohr's diagram consists of a pair of 
straight lines parallel to the shear stress axis. Coulomb's internal 
friction theory implies that failure occurs by sliding along a certain 
plane and the resistance to failure is composed of a constant shearing 
strength and an internal friction proportional to the normal stress on 
the sliding plane. The failure envelope consists of a pair of straight 
lines tangent to the pure shear and compression circles. Mohr generalized 
Coulomb's internal friction theory and assumed that the resistance 


to sliding is a function of the normal stress on the sliding plane. 
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The resulting failure envelope is curvilinear, the shape of which is 
to be determined by tests. Cowan's theory is based on a dual criterion 
of failure. He combined Rankine's maximum stress theory and Coulomb's 
internal friction theory; the former for cleavage failure and the later 
for shear fracture. Zia's modification to Cowan's failure theory is 
aimed at a closer approximation to the Mohr's generalized internal 
Triccione theory=, Fine modification predicts the failure of concrete transi- 
tional from tension failure to shear compression failure. Failure 
theories for concrete subjected to combined stresses have been dis- 
cussed in the literature (Bresler and Pister 1955, Zia 1961 and 1968). 
It is to be noted that the difference between these theories is greater 
when the compression (prestress) is high in which case failure is likely 
to be governed by principal compressive stress rather than the principal 
tensile stress. But when the prestress is of the order of 0.5 be 
which is more practical and about the maximum permitted by the ACI 
the theories due to Rankine and Cowan are identical and modified Cowan 
theory differs by a maximum of only about 12 percent. 

None of the above mentioned theories consider the effect 
of the intermediate principal stress and the strain gradient. It has 
been suggested by several investigators that the effect of the inter- 
mediate principal stress is significant in failure theories for 
concrete subjected to complex three dimensional states of stress. 
Bresler and Pister (1958), Seth (1958) McHenry and Karni (1958) 
have proposed failure theories taking into account all the principal 


stresses. In many cases such as in relatively slender beams, it is 
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considered accurate enough to use a failure criteria based on biaxial 
Stress condition. It is interesting to note that in the three dimen- 
Sional finite element analysis of cracking strength presented later in 
this dissertation, the intermediate principal stress under combined 
loading was negligibly small. Hence it can be concluded that the in- 
termediate principal stress can be safely neglected in such studies. 

Further, the above theories do not consider the effect of 
strain gradient. For example, in the thin web of an I section the 
Strain gradient is very high. However, it is small in a hollow con- 
crete beam. 

Recently Kupfer, Hilsdorf, and Rusch (1969) investigated 
experimentally the biaxial strength concrete over the entire range of 
biaxial stress combinations. They used square concrete plates 
(7.9 x 7.9 x 2 in.) subjected to an in-plane loading. Historically 
in tests on concrete under biaxial tension and compression hollow 
cylinders subjected to either torsion and axial compression or to 
internal hydraulic pressure and axial compression were used. For such 
stress states the test data from various investigations do not deviate 
from each other considerably. However, in biaxial compression tests 
the friction between the bearing platens and the specimens was suffi- 
ciently large to yield strengths of unusal magnitudes; for the case 
of equal compressive stresses in both principal directions the result- 
ing strengths reported were in the range of 80 to 350 percent of the 
uniaxial compression strength for identical test specimens. Kupfer et 


al replaced the usual solid bearing platens by "brush bearing platens" 
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which were flexible enough to follow the concrete deformations without 
CauSing appreciable restraint of the test specimens. They concluded 
contrary to the previous investigations that the maximum strength in- 
crease under biaxial compression was only 27 percent in excess of the 
uniaxial compression strength. Their tests in the region of combined 
compression and tension corroborate the findings of other investigators; 
i.e. stress at failure decreases as the simultaneously acting tensile 
stress is increased. The strength of concrete under biaxial tension 
iS approximately equal to its uniaxial tensile strength. 
3.2.3 Tensile Strength of Concrete 

It is evident from the literature review that the tensile 
Strengths used in torsion research by various investigators deviate 
from each other significantly for the following reasons. 

The tensile strength of concrete depends mainly on the state 
of stress at a point and also on the strain gradient. Beams acted 
on by combined loadings are subject to widely varying combined stress 
conditions; and different strain gradients. Further, concrete is not 
a perfectly elastic material. 

Each method of test available for evaluating the tensile 
strength of concrete has a certain stress situation associated 
with it. For example in the direct tensile test, there is a uniform 
stress along only one axis without any strain gradient. This is the 


only test which does not use any assumptions of elasticity or plasti- 
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city to calculate tensile stress. In the flexure test, tensile stresses 
are along one axis only with no stresses on the orthogonal planes; but 
there is a strain gradient. The splitting tensile strength test is 
characterized by a stress system which varies from biaxial compression 
immediately beneath the packing material to compression-tension over 
the most part of the diameter. The magnitude of the compressive stress 
varies and is about three times the magnitude of the tensile stress. 

The magnitude and sign of the third principal stress along the longi- 
tudinal axis is unknown. The measured strength may correspond to a 
complex function of the strength under the multiaxial stress system. 

Each investigator presumably chose a tensile test having a 
stress distribution which was nearer to the stress condition in his 
test beams under combined loading though it may not necessarily be 
satisfactory for all combinations of loadings. 

To predict cracking strength by method 1, the tensile strength 
of concrete taken in the present investigation, is equal to splitting 
tensile strength. It is believed that this tensile strength represents 
qualitatively the biaxial stress condition of compression-tension en- 
countered in beams subject to combined loading. The relationship be- 
tween splitting tensile and compressive strengths for each beam is shown 
in Fig. 3.3. For each beam the concrete compressive strength is obtained 
as the average of three cylinder tests and the tensile strength as the 
average of two split cylinder tests. The visually best fit line through 
the data yields a value of 6 v ce for the splitting tensile strength of 


concrete; this value is adopted in the cracking analysis of method 1. 
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3.2.4 Analysis 


location can be calculated as follows. 


sidered positive. 


Fig. 3.4. 
he 


The stress components causing the initial crack at any 


Shear stress due to torque 
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Shear stress due to flexural shear 
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Normal stress due to prestress 


At top face, 


fo = - z+ _+ 
At neutral axis, 
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At bottom face, 
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Tensile stresses are con- 
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The distributions of these stresses are shown in 


(3.4) 


(376) 


(320) 


(2.8) at “vr 
reoiszerq ot sub zeaidz TommoM 66 
soe? qot 3A 
(-€) | ipej-=4 
vetxe Gardusn JA 7 
(8.£) | 7 


(8.6) 


93 


SS3a1S 8Vans 
IWNOISAYOL 


SJOYOI SNOIYVA OL JING ONINIVYD OL YOTYd 
NOILIAS 3LFYINOD GISSIYLSIYd NI NOILNGIYLSIG SSsulS ay ee la 


LNIWOW $$3yLS3ud 
S$S3Y1s ONIONI@ OL JING Ol 3NG WvV38 3H1 JO 
YV3IHS 1VYNXII4 SS3¥LS IWWYON S$S3dlS IVWYON NOILD3IS SSOdD 


worT332 3TRAING? AFZ2IATLING HI -MOLTUAIATZTO-2e3ATe #.£ 917 
239H02 2UOTSAY OT 3G OMINDARD OT AOTAA 


94 


4. Normal stress due to flexural moment 


At top face, 


MY, 
104 Gr (359) 
At neutral axis, 
an = 0 (3.10) 
At bottom face, 
MY. 


At diagonal tension cracking, the stress components are related as 


(f+f 
Ty er hy Of ge se (3.12) 


where fy is the tensile strength of concrete equal to 6 Y ee Ex- 


follows: 


pressing equation (3.12) in terms of forces and ratios W(= T/M), and 
6(= T/Vb,.) yields 
es Fry 
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Solution of this quadratic equation yields the cracking 
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95 
torque. The equation can be specialized for each of the following three 
critical locations. 

le Gracking at the top face, 
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2. Cracking at neutral axis 
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3. Cracking at bottom face 
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Pal 2 ibe Selcrg ele aS ihe Ue Oxge 
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Any other critical locations can be examined by substituting appropriate 
quantities in the equation (3.13). 
The inclination of the initial tensile crack 6, with respect 


to the longitudinal axis can be obtained by the following equation 


Ketan aie 


tan(180-26) = FF 72 rales 


solution of equations (3.14), (3.15) and (3.16) furnish three dif- 
ferent values of torques and the cracking torque is chosen as the 


smallest of these three values. 
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3.3 Method 2 
3.3.1 General Remarks 

The finite element method has proved to be a convenient and 
versatile tool to aid in understanding the behavior of structural 
systems. It has been demonstrated in the past that this method can 
be applied with considerable success to various static and dynamic 
problems. Recently the method has been extended to include both ma- 
terial and geometric nonlinearities. 
3.3.2 Finite Element Method 

An attempt is made herein to predict by the finite element 
method, the cracking strength of reinforced and prestressed concrete 
beams under combined loading. Torsional cracking is a three dimensional 
phenomenon; accordingly the finite element model for such a study should 
be three dimensional. Clough (1969) compared various three dimensional 
elements in terms of their structural performance and computational 
efficiency, and recommended an eight node isoparametric hexahedron with 
24 degrees of freedom for the analysis of elastic solids. These iso- 
parametric elements are completely isotropic and can be fitted to 
quadrilaterally curved boundaries and were shown to be very effective 
finite elements. If the finite element system is analyzed by a dis- 
placement method of analysis to yield nodal point displacements and 
internal stresses, the models are stiffer in general. It was shown 
by him for this element, that any possible modification of the assumed 
interpolation functions would change the flexibility characteristics 


only to a small extent. Clough further showed that for structural 
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systems in which shear stresses are predominant the 8 node element 
is preferable to a 20 node refined element, but it is not very effective 
in representing plate bending behavior where shear stresses are low. 

The three dimensional eight node solid element used in the 
program in this dissertation is based on the element developed by 
Irons and Zienkiewicz. In analyzing beams of usual proportions, the 
bending characteristics of this element are observed to be poor even 
at aspect ratios of the order of 2.5. Hence the element has been 
modified here by adding incompatible deformation modes as suggested 
by Wilson (1966); this improved the bending behavior of the element 
Significantly. These deformation modes are eliminated at the element 
level by static condensation. The necessary integration for the 
element stiffness formulation cannot be done explicitly for this 
element and hence numerical integration is generally used. Two point 
Gauss quadrature is used in the present investigation. This gives exact 
values for rectangular prisms. Element stiffness formulation is skipped 
for identical elements. The structure stiffness matrix is generated 
in blocks and solved in blocks by direct Gaussian elimination. For 
solving the system of simultaneous equations an equation solver deve- 
loped by Gopalakrishnayya (1973) is used in this program. The pro- 
gram developed herein can solve various load conditions with the 
stiffness matrix inversion made only once; which results in con- 
siderable economy of computation time. The facility of automatic 
nodal and element data generation has considerably reduced the input 


data. 
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The theoretical background, essential features and listing 
of the finite element program developed herein are presented in 
Appendix B. The analysis developed is an elastic analysis, and the 
basic assumption made is that concrete is homogeneous and perfectly 
elastic until the formation of an initial crack. Before crack forma- 
tion the stresses in the steel induced by the applied loads are neglected. 
Cracking is considered to commence when the principal tensile stress 
equals the tensile strength of the concrete. 

3.3.3 Elastic Constants and Tensile Strength of Concrete 

For forming the modulus matrix, knowledge of the basic 
elastic constants of modulus of elasticity E. and Poissons ratio wv 
of concrete is necessary. The modulus of elasticity of concrete is 
taken as 57,000 vY us in accordance with ACI 318-71. Poissons ratio 
is assumed as 0.15 or 0.16. It is to be noted that an error in the 
assumption for the values of E. and v affect only the deformations, 
not the stresses; in statically determinate structures such as the 
beams analyzed herein. 

The tensile strength of concrete is taken as 7.5 v ie for 
the reasons listed below. 

1. The finite element analysis at best is an elastic analysis. 
However the torsional shear stresses are distributed at least semi- 
plastically. Hence a higher value of tensile strength than the one 
assumed in method 1 is appropriate. 

2. ACI 318-71 specifies modulus of rupture to be equal to 
yoy le : 
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3. CEB suggests modulus of rupture as 9.5 7 fe 

4. Woodhead and McMullen (1972) in their elastic analysis of 
prestressed rectangular beams for cracking used a value of 7 / Te 
This value consistutes a lower bound to the values that would result 
from the equation (2.17) of Chander et al (1970) if the elastic theory 
and their equation were to give the same torque. 
3.3.4 Finite Element Layout 

Unfortunately three dimensional analyses demands very large 

computational effort, even for problems of very small size. Hence it 
is not surprising that a compromise must be achieved between the de- 
Sirable accuracy and the computational time. The selection of the 
present finite element layout for rectangular beams in the literature 
and author's I beams was based on the experience gained by testing the 
program on similar smaller problems. It is believed that the mesh 
Size in terms of bandwidth and number of equations used produces an 
adequate degree of accuracy without further increasing the already ex- 
cessive computational effort. For the rectangular beams analyzed the 
maximum aspect ratio of the element varied from about 3 to 6, while 
in the case of I section beams it was 4. Outside the test zone slightly 
higher sepeee ratios were used. The number of equations did not exceed 
2500. With the usual solution techniques the computation effort is 
roughly proportional to the number of equations and to the square of the 
band width. The computations are performed in single precision for 
reasons of economy of memory and computational time. With large systems 
of equations round off starts to play an important role. This can be 


remedied significantly by doing the computations in double precision. 
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But this roughly doubles the core storage requirement. Another al- 
ternative is to keep the number of equations to the minimum possible 
which is adopted in the present case. 

The general layout and discretization of the beams and the 
boundary conditions at either support are shown in Fig. B.2 and Fig. 
B.3. The beams were simply supported in the tests and the analysis 
was carried out accordingly. The right hand side support is torsionally 
restrained, allowing flexural rotation, and longitudinal translation. 

To acnieve this all the nodes in the bottom face have displacements 

in all except the longitudinal direction equal to zero. All other 

nodes at this support have sideway displacements equal to zero because 

the cross section is torsionally restrained. The torsional load is 
applied at the left hand support and the nodes on this cross section 

are allowed to displace except the node at the centre of the cross section 
which is fixed in all directions. Both end sections of the test I 

beams of this investigation were made rectangular to facilitate easy 
gripping in the torsional heads. In the finite element analysis these 
were idealized as I section beams throughout the length for simplicity. 

The boundary forces must be distributed exactly in the same 
fashion as the calculated stresses if the exact stress distribution is 
to be obtained throughout the member. If the distribution of these 
forces does not conform to the exact solution but still has the same 
statically equivalent system of forces, then the resulting stresses are 
significantly different only in the vicinity of the applied forces. 


Inspite of these local perturbations, the exact stress distribution 
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develops a short distance from the areas where the forces are applied. 
To have the correct stress distribution as near as possible to the 
loaded ends, the transverse, prestress, and torsional loads are ap- 
plied at the nodes to reflect the exact stress distribution reasonably. 
Prestress was applied as external load. In the present case the stresses 
are in good agreement with the exact stresses at distances about one 
and one half times the depth of a beam or beyond the third cross-section 
from the ends of the beam. However, in each beam the end regions where 
the perturbations distort the stress distribution correspond to non- 
test zones which are of little interest in analysis. 

3.3.5 Procedure of Analysis 

In this method of analysis the input consists of a set of 
forces and the output in stress and strain distribution at every node 
and element face. As the cracking loads are known from tests before 
hand these loads are applied to obtain the stresses. This is not 
necessary in that any arbitrary set of loads in the same torque-shear 
proportion used in the test could be applied. The output consists of 
nodal deformations, the six components of stress and strain at the centers 
of element faces and average values at the nodes. It also includes 
principal stresses, strains and their directions. From this the re- 
sulting stress and strain distributions associated with the applied 
cracking loads are obtained. 

The loading combination for each beam test corresponds to a 
load vector consisting of transverse, prestress and torsional loads. 


If the beam geometry, material properties, and the boundary conditions 
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do not change and only the load vector changes, solutions for 10 to 
12 load vectors or load combinations are obtained by inverting the 
stiffness matrix just once. It took 2132 seconds to get solution for 
the first problem, while for the remaining nine problems, the time re- 
quired was 3946 seconds when ten problems with a band width of 150 
and number of equations of 2340 were solved at a time. 
3.3.6 Stresses 

Once the stress distribution at the cracking loads is known, 
the maximum principal tensile stress is located in the test zone either 
at a node or centre of an element. Cracking is assumed to initiate 
at that point. The maximum principal tensile stress occurred in the 
bottom face if the bending moment was predominant, and in the web or in 
the top face if the amount of torsion was high. The location of this 
principal tensile stress agreed with the test location in most of the 
beams. The cracking starts when the maximum principal tensile stress 
occurring in the test zone equals the tensile strength of concrete. 
These principal stresses are presented for each beam in Table 7.1 and 7.3. 
From the nodal deformations, the flexural deflections, and the angle 
of twist in a certain span length are obtained. 

The beams analyzed were subjected to various possible 
combinations of moments and forces. Of the resulting three principal 
stresses, the third was the tensile stress and related to the 
tensile strength of concrete. The first principal stress was 
usually compressive having varying dearees of magnitude. The 
intermediate principal stress was either compressive or tensile 


and had a magnitude of the order of 10 psi or less in all the beams. 
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This study has thus reinforced the possibility of neglecting the 


intermediate principal stress in determining the failure criteria 


of concrete for slender beams; because its magnitude is very small. 


Cracking is assumed to occur when the maximum principal 


tensile stress equals the assumed tensile strength of the concrete. 


However in the theoretical computations the test cracking load 

may produce a tensile stress different from the assumed tensile 
strength. In such a case a different load will produce a tensile 
stress equal to the assumed tensile strength, which is taken to be 
the theoretical cracking load. This theoretical cracking load is 
computed eliminating the effect of prestress from the stress 
distribution obtained by the finite element analysis associated 
with the test load at cracking. In making this calculation it is 
assumed that the principle of superposition holds good and the 
location of the maximum principal tensile stress does not change. 


The equations used are given in Appendix B. 


103 


, 
7 - 


Foy 


= 


“ia nF 
° Pstghontag: 
‘ainou (ao area fend 
bao! ontviwhd deod 5} bnohehagnos. Atengot va 
ofianss hamuezs ont mor? snovethb 2zavt2 off ened pecieterl: bee 
atiansd 6 souborg (ftw beat snatottth «© 9209 6 rove at «Atenas: . 
ad’ of tates 8t Watdw (idonerae ot renad bomuzes fa ot ‘a, ae a 
21 beat prtdsésa (sofderosdd 2tnT */baot @itoe9 fsotderosny od : 
seats: aid mont zeastéorg to Jos¥is sad anttentatte behenoe 
bedatooeee’ steyfenk tnomets ottntt and yd bentstdo nord: | 
27 3 motyatuoles 2tAd patdem nt -pntaosys 46 bsol ‘test oid AotW 
ailt ine boop Zor nattteoqiaque To sfatoning SMd Joris ‘\enlan 
Spneds ton esub 22ast2 altenst Taqtantq aes of) to iokseet 
a sree al ni navig sv6 be2u snutdeues emt 


oy 


Pr. 


CHAPTER IV 
ANALYSIS OF ULTIMATE STRENGTH 


4.1 General Remarks 

In this chapter equations are developed based on skew bending 
theory to predict the ultimate strength of prestressed concrete I section 
members subjected to torsion, bending and shear. The I section members 
tested by the author in the experimental phase of the program, failed 
generally in skew bending about an axis inclined to the longitudinal 
axis, either in Mode 1 or Mode 2. In mode 1 the compression zone formed 
on the top face while in Mode 2 it was observed inside the web close 
to one of its faces. Previously in Chapter II Mode 1 and Mode 2 were 
referred to as bending mode and torsion mode respectively. It is evi- 
dent from the literature review that flanged beams failed mainly in 
the skew bending mode. To date no adeauate theory is generally avail- 
able to compute the ultimate strength of prestressed I beams under 
combined loadings. The strength of such beams depends considerably on 
the degree with which the component rectangles act together, in other 
words, it depends on the reinforcement detail. In the sections that 
follow, the skew bending theory is extended and applied to I section 


members. 
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4.2 Analysis of Bending Mode of Failure 
4.2.1 Assumptions of Analysis 
The following assumptions are made in the analysis: 

1. The compression zone is rectangular; it is adjacent and parallel 
to one of the flange widths of the beam. 

2. The failure surface is defined by continuous tensile cracks 
at a constant inclination 6 with respect to the longitudinal axis, on 
the three sides of beam and joined on the fourth side by a compression 
zone. The crack inclination 6 at failure is same as the initial crack 
angle. 

3. The strain across the failure crack is proportional to its 
projected distance from the neutral axis measured in the vertical plane, 
referred as the 8 plane, containing the compression zone. 

4. The concrete compressive Strain, €,at the extreme fibre in 
the direction of the longitudinal axis at ultimate is given by a limiting 
value 0.0035. 


5. The stress-strain relationship of concrete in compression is 


given by 
£ Ee 
Oe i ee (4.1) 
re) 

where o = stress 

© -="ctyain 

E =e Omax! “0 

k. = (1.2-0.05 ie) < 1 and f, in ksi 
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6. The shear force in the compression zone is neglected. 

7. A premature shear failure does not occur. 

8. The dowel forces in both the longitudinal reinforcement and 
the stirrups are neglected. 

9. The mild steel reinforcement has a well defined yield plateau. 

10. The tensile force resisted by the concrete is neglected. 

11. The failure zone contains no external loads. 
4.2.2 Discussion of Assumptions 

Assumption 1] states that the compression zone is rectangular. 
This is true only when flexural moment and transverse shear act. 
Presence of torsion makes it either trapezoidal, triangular, or part 
of the compression zone may extend onto one of the vertical sides de- 
pending on the relative proportions of the torque. Rectangular com- 
pression zone was used by Lessig (1958) and others for reinforced con- 
crete beams and by Henry and Zia for prestressed concrete beams. Wood- 
head and McMullen (1972) used an inclined neutral axis, the inclination 
depending on the relative magnitudes of torque and moment and found that 
this refinement resulted in little difference to the theoretical ulti- 
mate strength. It is therefore considered adequate to assume a rectangular 
compression zone. 

Assumption 2 defines the failure surface. The three lines 
composing the spatial crack are assumed for simplicity to be straight 
and at the same inclination. The crack inclinations on the south and 
north faces differed mainly as the torsional and flexural shear stresses 
were additive on the former and subtractive on the latter. Crack in- 


clinations computed based on the uncracked section, are assumed to hold 
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good at ultimate. This is an approximation, however, a slight change 
in the inclination has not significantly affected the calculated ulti- 
mate strength. 

Assumption 3 is approximate and has been used by Henry and 
Zia (1971) and Woodhead and McMullen (1972). 

Assumption 4 relates to the useful limit of strain for concrete 
in compression under combined loading. This magnitude under bending 
only also depends on several factors. A limiting strain value of 
0.0038 was assumed by Hognestad, 0.004 by Warwaruk, and 0.0034 by Lin. 
As a lower bound the ACI suggests 0.003. Evans and Khalil (1970) ob- 
served that the longitudinal strains in the compressed concrete were 
essentially the same under pure bending and under combined bending 
and torsion for equal bending moments. Hence the same ultimate strain 
in concrete was used by them both in pure bending and combined loading. 
However they did not specify which value they adopted. Woodhead and 
McMullen (1972) used a value of 0.0034. In the analysis developed here- 
in a value of 0.0035 is used as an average between 0.004 and 0.003. 

In assumption 5 equation (4.1) expresses the stress-strain 
relationship of concrete in compression; this is due to Desai and 
Krishnan (1964). 

It is to be noted that in this analysis of the six equations 
of equilibrium, two equations are satisfied. This has been the usual 
approach followed by the previous investigators. Of these two,one is 
the equilibrium of forces perpendicular to the 8 plane, and the other 
is the equilibrium of moments about an axis in the 8 plane which is 


parallel to the top face. As the 8 plane is vertical, the transverse 
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shear acting vertically on the plane does not enter in either of the Mode 
1 equations. The transverse shear is resisted in part by the compression 
zone and the dowel action of the reinforcement. The sixth and the 
eighth assumptions concerning the neglect of this shear resistance are 
introduced to simplify the analysis. Since the failure section is 
assumed symmetric the forces in the two legs of the stirrups on either 
side of the web must be equal and opposite. It can be seen that torque 
is partly carried by the two legs of a stirrup developing equal but 
opposite forces. Also transverse shear is in part carried by the same 
two legs of a stirrup with the development of equal forces. Thus the 
transverse shear does not enter into the theoretical expression as the 
equilibrium of the vertical forces is identically satisfied. In essence 
no attempt is made to satisfy the vertical equilibrium of forces. 

Assumption 7 is introduced to exclude the possibility of a 
premature shear failure. If a beam acted on by only bending and shear 
fails in shear, such a beam probably develops a shear mode of failure 
when acted on by torsion, bending and shear. To date the question of 
shear strength alone is far from settled. Consequently the problem 
of combined torsion and shear is even further from solution as this is 
a much more complex problem. The effect of transverse shear is dis- 
cussed later in Chapter VII. 

Assumptions 9 through 11 are used to simplify the analysis. 
4.2.3 Scope and Steps of Analysis 

The analysis can be applied to under-reinforced as well as over- 

reinforced beams. From the assumptions it is clear that failure is 


considered to occur when the concrete crushes on a plane inclined at 
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an angle 8 to the longitudinal axis due to the attainment of a limiting 


strain along this axis. The flanged beams in which the torsional 
cracking splits off the outstanding parts of the compression flange 
are analyzed assuming that the width of the beam in the compression 
zone is equal to the web width. For others full flange width may be 
used. The analysis details are similar to those used for computing 
the ultimate flexural strength of an over-reinforced prestressed beam. 
The equilibrium equations are established with reference to the 8-plane 
which is vertical and inclined to the longitudinal axis at an angle 8, 
containing the compression zone and the neutral axis. The strain dis- 
tribution in the 6 plane is shown in Fig. 4.1. Components of internal 
and external forces perpendicular to the 8-plane are equated to zero. 
Moments of the forces are taken in 8-plane about an axis parallel to 
the top face and passing through the point of action of concrete com- 
pressive force. With the maximum compressive strain obtained from 
equation (4.3), a depth of neutral axis is assumed in the 8-plane. 

The induced strains in all the reinforcements are then compatible 

with the maximum concrete compressive strain and the assumed depth 

of the neutral axis. fhe compressive forces in the concrete and in 
any steel present in the compression zone, and the tensile forces are 
determined perpendicular to the 8-plane and equilibrium checked. Dif- 
ferent values of the depth of neutral axis are tried until the forces 
equilibrate. After finding the neutral axis, moments of all the forces 
internal and external are taken about the axis mentioned above and 
equated to zero. Solution of the moment equation yields the ultimate 
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4.2.4 Inclination of Failure Surface 

At failure the beam sections adjacent to the failure crack 
rotate about a neutral axis inclined at an angle 8 to the longitudinal 
axis. The developed failure surface for the bending mode is shown 
in Fig. 4.3 from which the inclinations of the tensile crack and the 
compression hinge can be related by equation (4.2). The tensile crack 
inclination 6 at ultimate is same as that at the instant of initial 


cracking and is obtained from the cracking analysis using equation (3.17). 
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4.2.5 Compressive Force in the Concrete 

It is recognized that the equations of equilibrium can be 
established with reference to any plane. In the present analysis the 
8 plane is used instead of a plane perpendicular to the longitudinal 
axis as the torque cannot be involved in either of the equations of 
equilibrium if the latter plane is used. 

In Mode 1, the concrete on the top face of the flange is sub- 
jected to longitudinal compressive stresses and shear stresses due to 
flexural and torsional moments, respectively, and this state may be 
resolved into that of principal compression and tension. At failure 
the beam rotated about an axis in the B plane with crushing taking 
place across the 8 plane accompanied by an orthogonal tensile strain. 
Because there is tensile strain on the orthogonal plane the ultimate 
compressive strain perpendicular to the B plane is considerably less 


than the strain at which concrete fails under pure compression. Kupfer, 


V5 ak iw 


doev> owitst edt ot 92 ne 
fentbuttpnot st of 8 efynts ns $e b 


ot bas dasta eftenet ond | 

doers sttenag off .(S.4) dotisups bla aatoantd ba wevgicnal 

letdint Yo Snstent sit te Jord 26 omee et stemtstu 38 6 fans ar 

(Sf.£) notisups ontew efeyfons pntdsevs edt mov? bantssdo shor Ser it 
: 7 

(#S#,d) | i: ‘ a 

($.6) anes a * Bt @ ¢ > : 7 


-* 


7 at912n0d aft nf sowed evt2zestqmod 2.8.8 
ad ne> ywind?liups to enok¥sups of? tedg boxtngoasy et 21 

old 2teyfens tnoeoq eit nl .onstq yas of sonsystey dziw bedetideses 
fentbuttonol ad? oF xatyatbnaqreq enafq 5 to bestent beau 2t ansig & 
to enokisups at to eddie nt bevlovat ed sonns> sup70s ond 26 2txe 
| | bow at ansiq redtet otf Vt mutvdt Tropes 

-dua et spnef? afd to aoet got ad mo stevonao edz ,F sboM al 
of aub zezecrd2 ~e9H2 bos zaaserde avtzzerqmon Fentbustpnos od betzat 
ad yom otetz atdy bus .ylovidosqesy ,2insmom Tsnotz102 bis Teywxslt 
ewifs? 2A .neienad is AObezarqmod Isqroatyq Yo Jedd aint bev (enbe 
ontdss pnid2us rdfw onelg A edt nt etxe m6 Juodé bedstor masd ond 
nfavte sfizned Isnqpori as yd betrsqmsos snaia’s oft 220726 s96lq 
atawtity ait onwiq Tenopodate oft no ntaite sfizned 2t avodt 925998 
zest yldsobienos et Snetg 4 adt of wefustbnsayeq ntert2 svizssrqma 
Listqu% -nofe2orqnos oyug yhny effet stsronoo dotdw te ntprie odd neds a 


| 2 eure 
: | es 


Hilsdrof and Rusch (1969) showed that in biaxial compression-tension 
of concrete when the tensile stress was about 20 percent of the com- 
pressive stress failure occurred at a compressive strain less than 
0.0005. This is almost an order of magnitude less than the uniaxial 
compressive crushing strain which is of the order of 0.0035. It can 
also be noted from their tests that even when the orthogonal tensile 
stress was only 5 percent the failure compressive strain was less than 
O20 TS: 

Experimentally Evans and Khalil (1970) measured the longitu- 


dinal concrete strains on the vertical face along the depth of the 


Tz 


beam, for the beams tested under pure bending and for other beams tested 


in combined bending and torsion. They stated that the presence of 
torsion had no significant effect on the distribution of the longi- 
tudinal strains, and on the depth of the compressed concrete for equal 
flexural moments. Having recognized that the concrete was crushing in 
a combined stress state at strains far less than the pure compression 


Strain, Evans and Khalil introduced an empirical relationship given 


: : t 2 
by equation (2.30) i.e. Seer eee (M/M ip) 
where . Se compressive strain in concrete in combined stress, 
aoe; longitudinal compressive strain in pure bending, 
My = ultimate bending moment in combined loading, 


uns = ultimate pure bending capacity. 


For example if Mi = 0.7 Me and on is assumed as 0.0035, then the 


failure strain ee will be 0.0017 from equation (2.30). 
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Woodhead and McMullen also noted a drastic reduction of the 
compressive failure strain across the B plane. They obtained this 
strain by starting with the longitudinal compressive strain and re- 


solving it across the plane. i.e. e€, = strain perpendicular to the 


B 
plane = Eo sing. The longitudinal compressive Strain, €., was Ob- 
tained by the empirical relationship suggested by Evans and Khalil 
i.e. equation (2.30). 
In light of the previous research and the appearance of the 
tensile cracks on the top flange just before failure in the beams of 


this investigation failing in Mode 1, the ultimate compressive strain 


perpendicular to the 8 plane is obtained by the following equation: 
z ee 
Sih ee (A) 


For the assumed depth of neutral axis, this compressive force can be 
computed knowing the stress-strain relationship of concrete in com- 
pression given by equation (4.1). By numerical integration the com- 
pression force and its centre of action can be determined. The steps 
are shown in Fig. 4.4. 
4.2.6 Forces in the Longitudinal Reinforcement 

The prestress results in a compressive strain in the 


concrete. As the load is applied the strain in the reinforcement 
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MODE 2 


Note: 


North side - Side on which shear stresses due to torsion 
and shear are subtractive 


South side - Side on which shear stresses due to torsion 
and shear are additive 


FIG. 4.3 DEVELOPED FAILURE SURFACES 
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increases first an amount equal e_. and later a small quantity of 


ce 
strain sufficient to cause cracking in the unstressed concrete. How- 
ever the strain necessary to induce cracking in the unstressed concrete 
is of the order of 0.0001 and is neglected. Further increase in the 
strain in the reinforcement is defined as induced strain es and is 
related to the useful limit of strain in the concrete by the strain 
compatibility. The total strain thus is equal to E. z 28 + Ey. 
The induced strain can be determined by knowing the depth of neutral 
axis and the concrete strain at the extreme fibre. The maximum force 
in the strands is assumed not to exceed the force corresponding to a 
strain of one percent. 
4.2.7 Forces in the Transverse Reinforcement 

The number of stirrups intersected by the tensile crack on 
the bottom face can be obtained from Fig. 4.5. As the crack inclina- 
tion near the compression zone is approximately equal to 8, in estima- 
ting the number of horizontal legs intersected by the crack, 6 is re- 
placed by B. The strain in the horizontal leas can be obtained as 


/cosB where e is the strain perpendicular to the 8 plane at the 


“stp stg 
level of the stirrups. The forces in the vertical legs of the stirrups 
do not enter the force equation. However, their moment contribution 

is to be taken into account. Hence, the force in the vertical legs 

of stirrups is to be estimated. The maximum strain in these branches 
of the stirrups is obtained by equating the components of the strains 
in the longitudinal direction and in the vertical legs of the stirrup 


in the direction of the failure surface, assuming no relative dis- 


placement along the failure surface for small rotations (Gangarao and 
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Zia, 1970). This is shown in Fig. 4.8. The number of vertical legs 
of stirrups intersected by the tensile crack is obtained from 
(h-a,-a,)coto/s. The force in the vertical legs is assumed to vary 
triangularly with the maximum force at the bottom and the minimum at 
the neutral axis. The force is considered not to exceed the yield 
force in any part of the stirrup. Hence, if the strain exceeds the 
yield strain, the stirrup carries a force equal to yield force. 


4.2.8 Analysis 


Figure 4.7 shows the freebody diagram with the forces on it 
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in mode 1 failure. Components of all the internal forces perpendicular 


to the 8-plane are as follows: 


Compressive force in the concrete 


C 


* sing ta 
Tensile force in the strands near the tension face 
= y any. E. been te ated sing (4.5) 
Tensile force in the strand at the second level 
= Aus a legte, te. | sink (4.6) 
Tensile force in the strands near the compression face 
=a) Ang ae le.gteete;! sing (4.7) 
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Tensile force in the horizontal legs of the stirrups near the tension 


face 


cote |*cv N-X-ay 


6) S cosB X 


= A. E,(b, -2a cose (4.8) 


Tensile force in the horizontal legs of the stirrups near the compression 


face 


6) eaiecaes cosp (4.9) 
There are no components of external forces perpendicular to the 6- 
plane. 
Summation of forces of equations (4.4) through (4.9) furnishes 
an equation which can be solved by an iteration procedure to evaluate 
the depth of the compression zone. 
The moments of all the forces are taken about an axis in the 
B-plane which is parallel to the top face and passes through the point 
of application of the concrete compressive force. External moments include: 


Component due to torsional moment 
= T x cosp (4.10) 
Component due to bending moment 


= Mx sing (4eqa) 
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Internal moments include: 


Moment contributions due to prestressing strands on the tension face 


ce [ 


ps te, |sing [h-a,-k,x] (A424) 


Ese ce 
Moment contribution due to force in prestressing strand at the second 
level 

=A CE Leuetes ten ]sine [h-a,-kox] (4.13) 


Moment contribution due to force in prestressing strands near the 


compression face 


=hE nee E. le. te te; Ising [a3-k,x] (4.14) 


Moment contribution due to force in horizontal legs of the stirrups 
near the tension face: 


(h-x-a,) 


= a\eeE (b,,-2a¢) oe = re oe (h-ag-k,x) (4.15) 


Moment contribution due to force in horizontal legs of the stirrups 


near the compression face: 


(a,-x) 
= 


b, -2a¢) 3 cots (ap-k,x) (4.16) 


Le EY ( CV Sx 


Moment contribution due to force in vertical legs of stirrups 


ee cote h-x-ay,  (h-aq-a) 2 


Piet Sistine or x S 
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cote 3(h-x-ay) cotesing (4277) 
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In equation (4.17) (h-ay-a,) cote/s gives the number of 
vertical legs of stirrups intersected by the crack. As the strain 
and hence force varies linearly with a maximum value at the level of 
the horizontal legs of the stirrups in the bottom face and a zero 
value at the neutral axis, the strain in the stirrups intersected at 
different levels along the height is taken as the mean of the fore- 
mentioned maximum and zero strains. The factor 1/2 accounts for this. 
The moment contribution due to the vertical legs on both sides of the 
web is same and is taken care of by the factor 2. The stirrup forces 
are distributed triangularly and the lever arm of the resultant of 
these forces from the neutral axis is to be accordingly calculated. 
This explains the presence of the factor 2/3. 

The external moments of equations (4.10) and equation (4.11) 
are equated to the equation (4.12) through (4.17) to establish the 
equation of moment equilibrium. The flexural moment can be expressed 
in terms of the torque and the only unknown, the ultimate torque can 


be solved. 


4.3 Analysis of Torsion Mode of Failure 
4.3.1 General Remarks 

In rectangular beams failing in torsion mode, the inclined 
compression zone forms close to one of the vertical sides and the two 
parts of the beam at failure rotate about the compression hinge. 
But for I beams, the literature review reveals that each investigator 
considered such a compression zone was either closest to the vertical 


sides of the flanges or in the web close to one of its sides. The 
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author noted the following observations during his I beam tests with 
reference to the location of the skew bending axis. The tensile 

cracks on one side of web did not penetrate completely to the other 
side i.e. tensile cracks on either side of the web were oriented in 
opposite directions. Shear-compression distress was observed clearly 
on the sloping portion of the flanges and faintly in the web. The 
vertical legs of the stirrups away from the tension face of the web 
recorded tensile strains. The compression louise On the vertical 
Sides of the flanges was not evident as suggested by Gvozdev et al 
(1968) and Woodhead and McMullen (1972) though it was clearly visible 
in the flanges of most of the beams tested. From these observations 

it was concluded that the shear-compression zone is along the sloping 
portion of the flanges and near the vertical face of the web, the magnitude 
and the centre of action of corresponding forces is very difficult to 
establish precisely. It is equally difficult to establish the axis 
about which the rotation at failure takes place. However to develop a 
simple, approximate, reasonably satisfactory method for estimating the 
ultimate strength in the torsion mode, it is assumed in the following 
analysis that the inclined skew bending axis is inside the web parallel 
and close to one of its faces. The presence of the outstanding flange 
on the compression side is neglected. Failure is considered to occur 
when the strain in the vertical legs of the stirrups on the tension 
side reaches yield strain. The strain in other reinforcements is 
assumed compatible with this yield strain of the stirrup leg. In addi- 


tion to assumptions 3 and 6 through 11 of section 4.2.1 the following 
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assumptions are made. 
4.3.2 Assumptions of Analysis 

1. The rectangular compression zone is parallel and adjacent 
to one of the vertical sides of the web extending throughout the height 
of the beam. 

2. The depth of the compression zone B, x computed according to 
the ACI equivalent rectangular stress block is very small. 

3. The failure surface is symmetric skewed and is defined by 
continuous tensile cracks at a constant inclination 9 with respect 
to the longitudinal axis on the three sides of beam and joined on 
the fourth side by a compression zone. The compression zone makes an 
angle B with the longitudinal axis. The crack inclination 6 at failure 
is same as the initial crack angle assuming that the initial crack occurs 
in the web. 

4. The transverse reinforcement is uniformly distributed along 
the length of the beam. 
4.3.3 Analysis 

The free body diagram for the torsion mode showing the failure 
surface and the forces acting on it, is shown in Fig. 4.8. The 8 plane 
is perpendicular to the sides and is inclined at an angle B with respect 
to the longitudinal axis. It contains the compression zone and the 
neutral axis and serves as a reference plane for establishing force and 
moment equilibrium equations. 
The assumed developed surface for the torsion mode is shown 


in Fig. 4.3. The relationship between the inclination of the tensile 
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crack and the compression hinge obtained from Fig. 4.3 can be written 


as 


tane = ( r ) tang (4°,18) 


As the compression zone is assumed rectangular and the bend- 
ing moment vector is orthogonal to the neutral axis, the bending moment 
does not enter the equations developed here for the torsion mode. Since 
the moment equilibrium about one axis is considered, only the bending 
component of the torque is taken into account. 

According to assumption 2, a concrete stress of 0.85 Us 
is assumed uniformly distributed over an equivalent rectangular com- 
pression zone. The depth of compression is 8] x as defined by ACI 
318-71. 

All the beams of this investigation failing in the torsion mode 
were under reinforced transversely and over-reinforced longitudinally. 
As noted before, failure is defined to occur when the transverse steel 
first reaches its yield strain and the strains in other reinforcements 
at that instant were assumed to be related to the yield strain of the 
transverse steel by the same strain compatibility noted earlier for 
Mode 1. The strain distribution in the 8 plane along the web is 
given in Fig. 4.2. The strain in the longitudinal direction at the 
level of the stirrups can be related to the strain in the vertical 
legs of the stirrups by the following expression (Gangarao and Zia, 


1970) 
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This relationship is obtained from Fig. 4.6, similarly the strain in 
the horizontal legs of the stirrups and the longitudinal strain are 


related 
Sr th tang (Az 21) 


For the torsion mode in the following analysis two equations 
of equilibrium, equilibrium of forces perpendicular to the 8 plane and 
equilibrium of moments about an axis of the 8 plane, are satisfied. 
The two equations allow the two unknowns namely the depth of neutral 
axis and the ultimate torque to be determined. Summing the forces 


perpendicular to the 8 plane on the free body: 
compressive force in the concrete = 0.85 i, B,xh/sing (4.21) 


Tensile force in the strands near the tension face 


b -a-x 


" W 
=) Bish Le eenery tané (5 mag sing CRD) 


If the longitudinal reinforcement is outside the web near the tension 
face then a4 is taken equal to zero. Tensile force in the strands 
near the compression face 


a+-X 
= y A_E bere iceet ty tang ( Jiesine (4.23) 
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If the longitudinal reinforcement is outside the web near the com- 
pression face then az is taken equal to zero. If nonprestressed 

steel is provided, its contribution can be obtained by substituting 
eee and oes equal to zero in the equations (4.22) and (4.23). Tensile 
force in the vertical legs of the stirrups near the tension face 


h-a -a3 


ENCE -—L-2) cote] cosp (4.24) 


tFtety Ll 
Tensile force in the vertical legs of the stirrups near the compression 
face 


h-a “a. 


—1) cote] cosB (4.25) 
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ALE. Lety 5 magax) IU 
Of the external forces only the flexural shear has a component per- 


pendicular to the 8 plane 
= V cos B (4.26) 


Using the relation between torque and shear & = NDS and assuming 
approximately one-half of stirrup steel cover to be equal to the depth 
of neutral axis in the equations (4.22) through equations (4.25), the 
expression for the depth of neutral axis may be written as follows: 
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a -a¢/2 
+ apse leseteceteeytanbly =a 72) 
(h-a,-a,) 


tFLety cote] cots 


ag/2 (h-a,-a ) 
Equilibrium of moments is established by taking moments of 
all the forces about an axis in the B plane at the level of the vertical 
legs of the stirrups close to the tension face of the web. 
Internal moments include: 
Moment contribution due toforce on concrete 
By xh 


= 0.85 ve Sing 


[b,-a¢- : 8.x] (4.28) 


Moment contribution due to force in prestressing strands on the tension 


face and compression face of the web 
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: W : 
= pst lesetecetpytanelp a ox J sing (a, ag) (4.29) 
a5-X 
= - Bipses' seco. sty becaeen (b -ae-as) (4.30) 


Moment contribution due to force in the vertical legs of stirrups 


near the compression face of the web 
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Moment contribution due to force in the horizontal legs of stirrups 


at top and bottom 


b -2a 


] 6 ] ' 
+ 2(5 ALE Let) CE et) (= [b,-a,-x]cots sing) (4.32) 


The term (b, -2a¢)coté/s gives the number of stirrups intercepted by 
the crack. The factor 2 accounts for the top and bottom legs of the 
stirrups, while the factor 1/2 averages the strain as it varies from 
a maximum value of ety on the tension side to zero on the compression 
side. The term 1/3[b, -a¢-x ]cotés ing gives the lever arm of the tri- 
angularly varying force in the stirrups. The external moments are 


shown in Fig. 4.8 and include: 


Due to the bending component of torsional moment 


= T cos B 


Due to the flexural shear assuming that it acts at the centre of the 


web 


Using the relation between the torque and shear. 6 = Way and equating 


the internal and external moments an expression for the ultimate torque 
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Equations (4.27) and (4.33) contain both x and i and hence have to be 
solved either simultaneously or by trial and error to determine the 


ultimate torque. 
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CHAPTER V 
EXPERIMENTAL INVESTIGATION 


5.1 Test Specimens 

The experimental program involved tests on forty four pre- 
stressed I section beams, which were subjected to various combinations 
of torsion,bending and shear. The beams were divided into four series. 
All the beams in each series nominally had identical properties and dif- 
fered only in the amount of transverse reinforcement. Among the series 
either the prestress, or its eccentricity, or both, were varied. The 
beams are designated by a system consisting of one letter followed 
by two numbers. The letter refers to the series, while the first 
and second numbers indicate the spacing of #2 stirrups, and an index 
of loading combination as shown in Table 7.1, respectively. For 
example, Beam A-3-2 related to A-series, had a stirrup spacing of 3 
in, and was tested at loading combination 2. For loading combination 
2 the torque-shear ratio (T/V bus equal to 5.2. 

The nominal details of the beam cross sections are shown 
in Fig.5.1. The arrangement of #2 stirrups used for web reinforcement 
is shown in Fig. 5.2. Actual beam dimensions were measured before 
test at three locations in the test zone, along the length of the beam. 
The positions of strands were similarly measured before casting and 
after testing. For calculation purposes dimensions at the measured 


location nearest to the failure section were chosen and are furnished 
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No. 2 STIRRUPS 


Nowe 3-@2.2' No. 2 No. 3 
STIRRUPS No. 2 STIRRUPS STIRRUPS | STIRRUPS 
6-@ 3" @ 3" OR 4" OR 5"OR | 13-@2" 6-@ 3" 
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FIG. 5.2 ARRANGEMENT OF STIRRUPS 
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in Table 5.1. The properties of all beams are summarized in Table 


Saas 


5.2 Materials 

Three sizes of prestressing strand were used in the test 
beams. Series A and C beams had two 1/2-in. and two 3/8-in. dia- 
meter strands while Series B had six 3/8-in. diameter strands. Series 
D beams were provided with four 5/16-in. and one 3/8-in. diameter 
strands. All were seven wire stress-relieved strands meeting the 
Specifications of ASTM A-416 with a minimum guaranteed ultimate 
strength of 250 ksi. The physical properties of typical strands are 
SNOW) IMBrNd. 5, 3. 

In practice such as used by the Highway Department of the 
State of Washington each two legged stirrup for prestressed I section 
beams subjected to torsion, bending and shear consists of four pieces 
involving two vertical legs and two horizontal pieces, for ease of 
fabrication of torsion stirrups around the tensioned strands. However, 
this may lead to difficulties in anchorage. Torsion specifications 
normally require a continuous stirrup as tensile stresses are deve- 
loped on all faces of the beam. In I girders because of the narrow 
web, the aspect ratio of the stirrup is such that the vertical legs 
would develop most of the torsional strength and short legs would be 
stressed less. Wyss et al (1969) used a 4 piece stirrup similar to 
that noted above and showed that it was not necessary to use closed 
rectangular stirrups for torsion web reinforcement in I section beams 
as the main contribution to the strength in torsion came from the 


vertical legs. Nevertheless in the present investigation, the two 
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a fe 5" 


ULTIMATE TENSILE STRENGTH (kips){ 39.800] 22.050] 16.050] 
LOAD AT 1% ELONGATION  (kips)| 35.200] 19.900| 14.225] 
MODULUS OF ELASTICITY (ksi)| 27500 | 28900 
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FIG. 5.3 LOAD-STRAIN CURVES AND PROPERTIES 
OF PRESTRESSING STRAND 
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TABLE 5.1 137 
DIMENSIONS OF TEST SPECIMENS 


ee 8 i a i ee ee ee eee 


Top Flange Bottom Flange Web Strand Cover 

B@amcNOjceme OV Css a a al oss) a ae 
Depth Depth Depth at Depth Depth at 
at end junction Width at end junction Width Depth Thickness Top Bottom 
with web with web 

A-0-4 W2eili2 2.19 2.94 7.06 2.03 2.78 eae 6.40 2706 1.66 0.97 
A-5-4 l2et2 ea22 ZA 7.03 2.06 2.81 Te5 6.34 2.56 ahs 1.06 
A-4-4 ete 2.24 2.99 Vresie 2.02 od 7 28 6.36 2.60 1.69 1.08 
A-3-4 12.12 rel 3.02 7.0 2.00 Om 7.43 6.35 2.59 1.81 1.06 
A-5-3 lem (asia 3.00 7.06 1.94 2.69 "7.40 6.43 2258 1.81 legit 
A-4-3 ler 2 (Bole 2.97 7.00 202 ah Woot 6.38 (ait) 1.81 1.03 
A-3-3 enazZ 226 SO) 7.00 1.93 2.68 7.40 6.43 a5) 1.80 1.00 
A-5-2 zine 2.20 ego 7.00 2.00 AS Wes) 6.42 2.56 eid 1.00 
A-4-2 Wa NZ eae? 2097 #05 Ze 2.77 asi! 6.38 2.56 1.81 1.00 
A-3-2 W2n 12 Zee 2.97 7.05 2.02 (Ae UT) aes 6.38 2n08) 1.81 1.00 
A-0-1 hee k2 2.238 2.98 7.00 1.98 ole oH 6.41 258 1.8) 1203 
A-5-1 2eii2 2.19 2.94 7.06 2.00 oils Pos 6.43 Zbo 1.69 1.00 
A-4-] Taiz Pee ZSDih 7.06 1.97 Bile Tieey 6.4] 2,03 1.81] 1.00 
A-3-1 Vee 2aN9 2.94 7.00 2.05 2.80 Hele 6.38 2.56 Wao 1.07 
A-3-0 Wee li2 Aorass 3.00 7.06 1.98 Ces 123% 6.39 53 1.65 1.10 
B-3-4 Wee 2 once 2.97 7.03 2107 2,82 oe 5538 2.58 1.70 1.10 
B-5-3 Weal 2ae5 3.00 7.00 1.96 Zi HP SH 6.4] AS qW75 Tees 
B- 4-3 lol Ao ie) 2,94 Wee Aekoiy 2.82 Poll 6.36 2.60 1.74 ewe 
B-3-3 Wei? Zolli enge 7.00 2.03 2.78 ool! 6.42 2.44 ets 1.03 
B-5-2 Welle 2ne0 3.00 Hale, 2.02 2 Hit 7,28 6.30 2.56 1.69 1.10 
B-4-2 ee Zee Zoi 7.00 2.00 (othe Ti 6.40 2.50 175 1.03 
5-3-2 leral2 2.24 2.99 7.00 i 2s! 2.73 dsc! 6.40 oe) 1.78 1.03 
B-5-1 In NZ oC 2.96 Tes 0, 2.08 2.83 Tachi 6.33 2.60 1.8] 1.00 
B-4-] lizaiiz 2.24 29S 7.00 25105 2.80 Pehl 6.33 Zn53 ils Ae 1.03 
B-3-] We We Bias 2597 7.00 2.08 2.83 i Sil 6.32 2008 1.8] 1.03 
B-3-0 W2NZ 2.20 3.00 7.03 2.01 2.76 Unies) 6.36 Cog 1.70 eal@ 
c-0-4 12.12 2.22 2.97 7.04 1.96 Cee | 7.40 6.44 (Esse) 181 1.00 
C-5-4 Re (aSt3) 3.00 7.00 95 och) I s'10) 6.42 2.47 1.81 1.00 
C-3-4 12.12 2.24 2.99 7508 2.00 Bells Tiwi 6.38 2.47 1.8] 0.97 
C-5-3 Ws We aie Gadi 7.00 age! (Aree) 7.40 6.42 OO) 1.88 6.91 
C-4-3 2 2.28 3)..03 7.00 1.94 2309 Techi/ 6.40 (Asr\0) 1.81 1.00 
C-3-3 eae 2.24 SS) 7.00 2.0) 2.76 7.40 6.37 2.47 1.84 0.98 
C-5-2 W252 2.20 Ae) 7.04 {hassle 2.70 7.40 6.47 2.50 1.8] 0.94 
C-4-2 Wan 2.24 2.99 7.00 2.01 2.76 avi 6.37 2.50 1.84 0.94 
C-3-2 Weole Bohl) 2295 7.00 1.96 AMA 7.40 6.46 2.47 1.94 0.88 
C-5-1 et facta} 3.00 7.05 ok 2.68 7.45 6.44 258 ll arAss 1.00 
C-4-] ale Aras 3.01 7.00 ifoeks 2 AG 7.40 6.41 2.44 1281 0.94 
C-3-] igen Bos 3.00 7.04 Wels Ba The 7.19 6.42 2.47 1.81 1.00 
D-3-4 [202 2.16 2.91 7.06 2.08 2.83 We 6.38 200 eS IOS /i2556 
D-3-4A 12.12 Cee 2.97 7.03 2.04 2.79 1525 6.36 2.56, 1.69 1.09/2.56 
D-3-3 Aerie Dale 2.97 7.00 2.06 2.81 Heli 6.34 (2s TS OS 72556 
D-3-2 W252 Cnee 2.97 7.03 2.03 2.78 7.16 Go 2.50 [Se RUQ7 27156 
D-3-1 WAN Oreo 3.00 dene 2.00 Als eee 6.37 2.56 Hts ool 
D-3-0 UiZrie 2.26 3.01 7.03 2.00 (sks) 1325 6.37 2.44 W661 9103/2559 


*A11 dimensions in inches 
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PROPERTIES OF TEST SPECIMENS 


Concrete 
Strength 
Beam No, 
lie sp 
(psi) (psi) 
A-0-4 4425 380 
A-5-4 5480 464 
A-4-4 4708 433 
A-3-4 4254 362 
A-5-3 4372 389 
A-4-3 5140 486 
A-3-3 4460 354 
A-5-2 4549 409 
A-4-2 4737 376 
A-3-2 4390 464 
A-0-1 5383 433 
A-5-1 4896 398 
A-4-1 4407 425 
A-3-] 5686 495 
A-3-0 4708 449 
B-3-4 4932 415 
B-5-3 440] 404 
B-4-3 5135 407 
B-3-3 473] 438 
B-5-2 5109 451 
B-4-2 443) 350 
B-3-2 4670 387 
B-5-] 4930 470 
B-4-] 4684 433 
B-3-] 5279 473 
B-3-0 4384 444 
c-0-4 4237 367 
C-5-4 4516 329 
C-3-4 Sse 484 
C-5-3 4440 336 
C-4-3 5330 423 
C-3-3 5539 493 
(-5-2 A930 489 
(eA 2 N6/b Mn 
Get? AYW?4 TM) 
C-h 1 4/34 Ais 
(-4-1 5005 409 
C-4-] 4519 409 
0-3-4 4879 393 
D-3-4A 4407 382 
D-3-3 5238 475 
D-3-2 4767 407 
0-3-1 4955 484 
0-3-0 4820 449 
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Strand 
(dia) 
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3/8 
3/8 
3/8 
3/8 
3/8 
3/8 
3/8 
3/8 
3/8 
3/8 
3/8 
3/8 
3/8 
3/8 
3/8 
3/8 
3/8 
3/8 
3/8 
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Effective 
Prestress 
(kips) 


Bo 
45. 
43. 
37h 
40. 
43. 
34. 
42. 
Bo) 
48. 
39. 
4). 
42. 
45, 
42. 
64, 
62. 
65. 
63. 
67. 
62. 
63. 
65. 
66. 
65. 
65. 
36. 
40. 
4). 
Soe 
37s 
38%. 

Sy. 

1/7, 

Nir 

WTB 

Wee 
38, 
42.8 


41. 
43. 
as 
40. 


40. 


79 
59 


(in) 


]- —- YO NNN NHYH HFT DY HF NY KH KH YD LY 


feof VaR RT ee a ee er a 


.039 
067 
035 
.920 
94] 
523) 
.872 
.118 
917 
. 266 
-036 
- 169 
063 
Sys 
-964 
.318 
«219 
568 
.290 


617 


044 
301 
- 360 
23) 
. 366 
pole 
.689 
.829 
agen 
79) 
. 708 
77) 
930 
96) 
65/7 
42 1 
-BSe 
870 
478 


oy Op Age at 
fpstlt ost) wipl’s (Gacy 
200 1209 0.168 64 
214 1378 0.154 10 
213 1296 0.169 im 
209 1099 0.163 73 
229 1204 0.173 113 
184 1322 0.156 35 
210 1004 0.144 107 
187 1295 OL 173 36 
224 1158 0.155 63 
167 1508 0.203 12 
194 1176 Onls5 106 
169 1278 0.157 37 
206 1283 0.179 59 
201 1398 0.149 9 
233 1262 0.168 10 
603 1653 0.242 10 
601 1576 0.262 13 
506 1765 0.235 im 
624 1639 0.255 12 
507 1830 0.242 10 
694 1489 0.261 20 
614 1636 0.256 32 
586 1669 0.243 34 
655 1653 0.261 2] 
59] 1687 0.230 29 
615 1662 0.275 12 
496 809 0.162 110 
512 925 0.168 82 
500 976 0.152 8] 
454 800 0.149 99 
503 83) 0.132 108 
498 863 0.130 87 
462 903 0.146 84 
443 891 Os 104 
514 ul/ 0.155 92 
49| 44 0.153 106 
482 8/0 0.142 103 
478 89] O59 92 
361 116] 0.162 8 
369 1116 0.175 9 
360 1203 0.155 7 
367 1185 0.169 8 
329 W1 0.150 10 
370 1109 0.159 9 
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horizontal pieces were dispensed with and the vertical leg was bent 
and extended as a horizontal portion of the stirrup. The two legged 
stirrups for the beams were formed by two separate bars bent as shown 
in Fig. 5.2. The horizontal portion of the bent bars extended to the 
end of the flange width allowing for cover. The stress-strain curve 
and general properties of the #2 smooth bars used for web reinforce- 
ment are presented in Fig. 5.4. 

The concrete mix was designed to yield a nominal 28 day 
compressive strength of 5000 psi; Type III high early strength cement 
was used to facilitate early transfer of prestress. The coarse ag- 
gregate was 1/2-in. maximum size crushed gravel. For each beam and 
its control cylinders the concrete was mixed in one batch of five 


cubic feet. Each batch consisted of the following mix proportions: 


Cement - Type III om DS. 
Sand 279 DS. 
Coarse aggregate 450 lbs. 
Water 1O-1DS. 


As the moisture content of the aggregates varied slightly the amount 
of water was adjusted in each batch to obtain a slump of the order 


of three inches. 


5.3 Fabrication of Specimens 

All the I beams were nominally 12 in. deep, 7 in. wide and 
10 ft.9 in. long. The depth of both top and bottom flanges was 2 in. 
at the outside edge and increased to 2.75 in. at the junction with 


the web. The thickness of the web was 2.5 inches. These beams were 
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FIG. 5.4 STRESS-STRAIN CURVE OF NO. 2 BAR 
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dimensioned to have approximately half the size of full scale girders 
tested by Wyss et al (1969). Both end sections were made rectangular 
to facilitate easy gripping in the torsional heads. To force the 
failure in the test zone, the portions of the beam beyond test zone 
were reinforced with additional mild steel bars, consisting of #3 
stirrups in the rectangular section portions and #2 stirrups in the 
transition and I section parts as shown in Fig. 5.2. Additionally 
the non-test zones were provided with inclined bars on both vertical 
sides of each beam in the probable direction of the principal tensile 
stresses, and four #2 longitudinal bars. The longitudinal bars ex- 
tended into the test zone for about 9 in. and terminated in a staggered 
manner to avoid a plane of weakness. 

The 26 ft. long prestressing bed in the laboratory, permitted 
fabrication of beams in pairs. The strands were individually tensioned 
by a 30 ton electrically-operated hydraulic jack; the force in each 
strand being measured by a load cel?. Each strand was stressed in 
excess by an estimated amount to account for losses and finally to 
yield the required tension at the time of test. After all the strands 
were tensioned, the mild steel cage was assembled. Stirrups instru- 
mented with KYOWA electrical resistance strain gages of 5 mm length 
(Type KFC-5-C1-11, Resistance = 120.0, Gage Factor = 2.12) were 
positioned at estimated locations of failure in the test zone. Budd 
GA-1 cement and a catalyst were used to mount each strain gage on the 
prepared surface of the stirrup. The gages were waterproofed with 
three coats of Budd GW-2 compound followed by two layers of insulating 


tape. Finally a coating of epoxilite compound was applied. For in- 
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strumenting strands, small cavities were created in the bottom and 
top flanges at two or three locations of probable failure, by insert- 
ing styrofoam plugs between the strand and the form before casting. 
After casting the plug holes were cleaned and foil gages of 7 mm and 
Type KFC-1-C1-11 (Resistance = 119.6, and Gage factor = 2.13) were 
installed on the strands using Budd GA-1 cement and catalyst. 

The mould consisted of two steel channels with wood inserts 
to form the I section and the transition portion to the rectangular 
section on either end of test zone. The steel channels were bolted 
to a heavy steel base. At each end of the form a mild steel plate 
was used as an end stop. 

Prior to casting the side forms were oiled and positioned. 
Concrete was mixed in the laboratory and each beam together with its 
control cylinders were cast from the same batch of concrete. One day 
after casting, the beams and the test cylinders were demoulded and 
Subsequently cured for six days by wrapping with wet burlap and plastic 
sheets. 

Prestress was transferred one week after casting. The strands 
exposed between the two beams were cut by heating them gradually over 
a length using an oxy-acetylene torch. The beams were then stored in 
the laboratory until time of test. To measure the prestress losses, 

a series of Demec gage points were glued to the surface of the two 
vertical faces and top face at the middle of the beam. Demec gage 
readings were taken immediately before and after the release of pre- 
stress in order to obtain the prestress loss due to elastic shortening. 


Before the test another set of readings were taken to determine loss of 
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prestress due to creep, shrinkage and relaxation. The difference in 
the readings taken before the application of prestress and before the 


test gave the total losses and hence the effective prestress. 


5.4 Instrumentation of Specimens 

The test specimens were instrumented to obtain the measure- 
ment of angle of twist, vertical deflections, and reinforcement strains. 

(a) Angle of Twist 

Two twist meters were mounted at each end of the gage length 
over which the total angle of twist was to be measured. Each twist 
meter was fabricated by two pieces of 1-1/2 in. x 1 in. aluminum channel 
sections. On the top channel a level bubble was mounted. The base 
and top channels were pin jointed at one end, and at the other end, 
a micrometer screw was mounted between them, for measurement of their 
relative movement. A spring forced contact between micrometer screw 
needle and top channel. The whole assembly was mounted on the beam. 
The pin joint and the micrometer screw needle were 15 inches apart and 
the micrometer least count was 0.001-in. The spirit bubbles were 
initially levelled by the micrometer screws, and as the beam twisted 
one bubble was more displaced than the other. Micrometer screws were 
adjusted until they were level again. The difference between the 
initial and current readings of a twist meter furnished the angle of 
twist of the section. The difference between the angles of twist 
of each twist meter was the total twist over the gage length. Ro- 
tations were measured over a gage length of 48 inches. Locations of 


twist meters are shown in Fig. 5.5. 
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(b) Deflections 
The vertical deflections for each of the beams were measured 
at three locations in the gage length as shown in Fig. 5.5. At each 
of the three locations two small curved nails were glued on the two 
vertical sides of bottom flange and from these nails two scales were 
freely suspended equidistant from the centre of twist of the uncracked 
beam. The least count of the scales was 0.01 inch. The six scales 
were read using two precise levels located on either side of the beam. 
The average of readings of the two scales at each location was taken 
as the deflection of the centre of the cross section. 
(c) Reinforcement Strain Measurement 
SR-4 electrical resistance strain gages were used to measure 
the strain in the vertical legs and occasionally in the horizontal 
portions of the stirrups. Three stirrups having gages were provided 
in each beam at probable locations of failure. Both on bottom and 
top strands on the south side of the beam, foil gages were used at 
two locations. The output from all these gaces was monitored by 


SS-24R switching and balancing box and SM-60 AT strain indicator. 


5.5 Test Setup and Procedure 

A general view of the test setup is shown in Plate 5.1 and 
a schematic diagram is shown in Fig. 5.5. The torsion test rig allowed 
the twisting and shear-flexural loads to be applied separately by means 
of two independent loading systems. 

The transverse load was applied by an Amsler Jack of 44 kip 


capacity, hung from the top girder of the loading frame. An electrically 
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operated Amsler pump fitted with a load indicator, was connected to 
the jack. The pump with its load maintaining device, maintained 
flexural load steadily at any desired level during the test. The 
Amsler jack transferred the load to the beam through a roller assembly 
and a pipe collar. The roller assembly was supported on the pipe 
collar which in turn was clamped to the test beam coaxially. The pipe 
collar provided free rotation of the test beam maintaining the roller 
assembly without any transverse displacement. Further, four hori- 
zontal tie bars, two on each side of the beam were used to keep the 
roller assembly in position. 

The east end of the specimen was supported by a twisting 
head, and the west end by a restraining head. The effective span of 
the beam was 9 feet. The twisting head permitted the specimen to 
rotate freely about its three principal axes while the restraining 
head allowed flexural rotation and translation, but prevented any 
torsional rotation of the specimen. For the twisting head the twist- 
ing couple was produced by two equal and opposite forces in two cables 
spaced 72 inches apart. The cables were connected between the torsion 
arms of the twisting head and the ends of a cross beam located below 
the load bed. The cross beam was designed to move freely within guides 
as it was loaded at its center by two 20 kip hydraulic jacks coupled 
in series. The jacks, operated by a manual pump, were coupled in 
order to increase the travel to 16 inches, which enabled the twisting 
head to develop a rotation of approximately 26°. The maximum rotation 
the twisting head is capable of is slightly greater than 26°. As the 


forces in the cables were essentially the same, only one of these was 
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measured using a load cell and an indicator during the tests. Reports 
by McMullen and Warwaruk (1967) and Mukherjee and Warwaruk (1970) 
present detailed illustrations of the torsion testing rig. 

The Amsler load indicator was preset at a reading appropriate 
to the self weight of the roller assembly and the pipe collar, as 
these were supported by the test specimen. Similarly, the weight of 
the cross beam had to be taken into account in reading the torsional 
load. Initially the cables hung freely as the cross beam rested on 
a support. The zero reading on the torque load cell indicator was 
adjusted, to correspond to a no tension situation in the cables, be- 
fore the cross beam support was removed. The first torsional load 
applied to a specimen was the weight of cross beam. 

Each beam was tested to failure by applyina loads ina 
series of increment which consisted of increasing to a predetermined 
level both the transverse load and the torsional load by a definite 
proportion, simultaneously. As the beam was subjected to a moment 
gradient, in the presence of transverse shear, the predetermined pro- 
portion was chosen in terms of torque-shear ratio which was constant 
in the test zone throughout the test. It took about one minute, for 
increasing the loads, after which, they were held constant, while 
steel strains, deflections, and rotations were measured and cracks 
marked. The holding period after each increment varied from two to 
four minutes. Smaller increments of loads were used as cracking and 
failure loads were approached. In certain beams, concrete surface 
strains were measured by Demec gage readings to relate to the strains 


recorded by strain gages mounted on the strands. In all instances, 
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load application was continued well beyond the stage at which peak values 
were observed. The location of the first crack and the ultimate failure 
plane were noted. Usually 20 to 25 increments were used prior to 

failure and the entire test took about 1-1/2 hours. Photographs of 


the crack patterns were taken at the end of each test. 


CHAPTER VI 
PRESENTATION OF TEST RESULTS 


6.1 Introduction 

The principal test results and the observed behavior of the 
beams are presented in this chapter. Additional experimental details 
consisting of torque-twist, and load-deformation curves for the whole 
range of loading, photographs of two side views of each beam and 
sketches of typical developed fa*lure surfaces are furnished in 
Appendix A. 

In Table 6.1 principal test results such as transverse load, 
torsional moment, bending moment, flexural shear both at cracking and 
ultimate are recorded. The location of the initial crack and the 
centre of the failure surface with reference to the position of the 
transverse load are also listed. The flexural deflections and the 
torsional rotations at cracking and ultimate conditions are tabulated 
in Table 6.2. Figure 5.5 shows the locations where these deformations 
were measured. 

Beams of each series were tested at one of the five nominal 
load ratios of Tbs =.6 Of ~, 10,4, 5.2; 2.229 and, 0.0, One beam .in 
series D was tested at a ratio of 6 equal to 20.8. In the beam designa- 
tions, the last numeral is the load combination corresponding to a 
particular load ratio. These nominal ratios did not include the dead 
load shear. However in presenting the load ratios w and 6 in Table 


6.1 dead load moment and shear were included. For any one of the 
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Beam ' 
No. fc 


(psi) 


A-0-4 
A-5-4 
A-4-4 
A-3-4 
A-5-3 
A-4-3 
A-3=3 
A-5-2 
A-4-2 
A-3-2 
A-0-1 
A-5-1 
A-4-1 
A-3-1 
A-3-0 
B-3-4 
B-5-3 
B-4-3 
B-3-3 
B-5-2 
B-4-2 
B-3-2 
B-5-] 
B-4-] 
B-3-] 
B-3-0 
C-0-4 
c-5-4 
292 
G2529 
C-4-3 
G02 
C522 
Gede2 
G=322 
C-5-1 
C-4-1 
C-3-1 4519 
D-3-4 4879 
D-3-4A 4407 
D-3-3 5238 
D-3-2 4767 
D-3-1 4955 
D-3-0 4820 


4425 
5480 
4708 
4254 
4372 
5140 
4460 
4549 
4737 
4390 
5383 
4896 
4407 
5686 
4708 
4932 
440] 
5135 
473) 
5109 
443) 
4670 
4930 
4684 
5279 
4384 
4237 
4516 
5132 
4440 
5330 
5539 
4930 
4676 
4525 
4734 
5003 
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PRINCIPAL TEST RESULTS 


Cracking Loads 


Ultimate Loads 


@Distance of initial crack from transverse load 
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i ne ict ae te ees ete tt re 
(6) (in) (in) (kips) (in-kips) (in-kips) (in-kips) (in-kips) (in-kips) (kips) crack failure 
oo 34 15 0 36.02 6.0 0.07 3738 6.6 0.01 TF Mode 2 
co) 8 28 0 47.94 6.3 0.05 54.95 6.3 0.05 TF Mode 2 
« 14 33 0 46.59 6.5 0.02 56.03 (i 0.07 he Mode 2 
2 29 29 0 38.64 6.3 0.05 50.64 6.3 0.05 TF Mode 2 
9.811 32 44 5.8 44.2] 74.) Wor'6 50.28 59.2 2.05 Web Mode 2 

10.153 32 28 6.5 50.28 83.5 2.0 56 nao. 101.7 nee Web Mode 2 
9.933 28 4) 655 43.34 79.7 Iara 56.37 Tne Dien Web Mode 2 
5.036 10 25: 1120 39.01 192.4 3.04 46.58 178.8 3.70 Web Mode 2 
5.626 12 22 9.8 38.33 166.5 2.69 41.27 169.8 3.29 Web Mode 2 
5.189 15 205 WekS 47.68 25.5 3.67 58.51 240.6 4.51 Web Mode 2 
ees 1A 7 19.5 25.04 267.5 4.52 36.23 428.8 6.55 BF Mode 1 
anes 2 10 924.0 25.08 276.5 Geos) 2 44.58 502.4 8.04 BF Mode 1 
2.437 6 Soe 27.84 292.2 4.39 39.9 422.3 (Jatsicl BF Mode 1 
2220 ORS 9 23.6, 25.99 293.4 4.70 43.9 502 7.91 BF Mode 1 

0 555 5.5 26.5 0 Dhee 4.06 (0) 593.6 8.89 BF No torsion 
© 30 35 0 Bsa, 6.2 0.05 61.9 5.9 0.08 TF/Web Mode 2 

T0540) 3025 3005 Heal 58.08 98.9 2.29 lp 55 104.4 2.42 Web Mode 2 

WO. 136) -33 33 7.6 60.68 97.1 2.40 65.88 104.9 2.60 Web Mode 2 

10.29 7 30 735 57.21 127.6 2.20 65.6 TAN 2 Zoo) Web Mode 2 
5.153 28 26 14.5 57.43 200.7 4.46 62.74 228.7 4.87 Web/BF Mode 2 
5.05 30 26 14.8 56). 35 188.0 4.39 62.75 233.3 4.97 Web Mode 2 
5.161 al 27 15.0 54.18 219.0 4.18 65.03 231.4 5.04 Web Mode 2 
2.204 10 9° 20/35 27.86 316.4 5.04 37.85 436.9 6.87 BF Mode | 
2.219 8 Qe 2555) 33.44 390.3 6.05 43.65 499.9 7.87 BF Mode | 
Deke 7 Giezien 29.72 S530) 5.38 43.89 486.3 7555: BF Mode 1 

0 6h. 6/5 32.9 0 358.2 5.39 0 724.6 11.02 BF No torsion 
” 17.5 28 0 41.26 6.6 0.0 44.7] 6.3 0.05 TF Mode 2 
oo 8 13 0 45.24 6.5 0.03 49.56 6.15 0.02 TF Mode 2 
20 32 18 0 45.25 6.1 0.06 58.73 8.8 0.00 TF Mode 2 
9.960 34 36 5.4 41.61 66.8 1.67 46.8) 70.6 1.88 Web Mode 2 
9.909 36 2 6.4 48.55 W320) 1.95 55.49 69.3 2.24 Web Mode 2 

10.360 40 21 6.5 49.4) 66nd (ar) 56.46 117.0 2.18 Web Mode 2 
5.328 S055 50 9.0 39.2 123.8 2.89 41.88 70.3 3.144 Web Mode 2 
§.113 17/31 31 11.0 44.75 147.0 3.49 47.68 156.5 34730 TF to BF Mode 2 
5.124 6 30 NO7: 36.84 193.24 3. Ne! 46.38 156.0 3.620 BF Mode 2 
Bisco 6 ) fats) 20.07 226.1 3.63 40.5 478.6 De BF Mode 
Lede VS05 i yatliet| 20.43 220). 3.69 39.58 463.5 7.08 BF Mode | 
Gee lazer: 8 al lips} 19.43 214.7 3.50 39.57 460.7 Helis BF Mode | 

” - 38 0 45.24 6.5 = 56.57 Bi 0.09 TF Mode 2 

19.597. 40 37 3.4 50.29 36.4 1.07 5977 45.4 lee TE Mode 2 
9.79 29 33 7.5 57.27 100.9 2.25 62.9 103.6 2.57 TF/Web Mode 2 
5.120 43.5/6 27 14.0 49.85 114.4 3.95 60.24 216.4 4.7) Web/ BF Mode 2 
2.218 5.5 6 22.5 18.58 22739 Josh) 41.78 501.2 755) BF Mode 1 

0 if 6 24.0 0 25.1 Boe 0 534.2 8.05 BF No torsion 
TF = Top Face 
BF = Bottom Face 


Distance of center of failure surface from transverse load 


“Calculated by statics at the location of initial crack 


dealculated by statics at the center of failure surface 
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EXPERIMENTAL DEFORMATIONS 


At cracking 
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TABLE 6.2 


At ultimate 
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0.145 
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*Rotation measured over a gage length of 48 inches 
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Deflection measured at a distance of 7 inches from the transverse load 
*Deformations at one load stage prior to the ultimate stage 
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following two reasons the actual magnitude of 6 deviated slightly 
from the nominal value. Firstly in each beam test torque and shear 
were applied simultaneously maintaining the same predetermined load 
ratio at the end of each increment throughout the whole loading range. 
However in some beams failure occurred while the torque was still 
being applied. Secondly, for some beams in applying the increment 
leading to failure, a very large rotation was necessary to reach the 
predetermined torque. In this process of twisting the failure torque 
sometimes slightly exceeded the predetermined torque corresponding to 
the load increment. The ultimate loads presented in Table 6.1 relate 


to the centre of the failure surface. 


6.2 Behavior Under Test 

In all the nine beams tested under pure torsion, the initial 
cracks formed in the middle of the top face at torques ranging from 
76 to 96 percent of the ultimate torques. In the case of B-3-4 initial 
cracks occurred simultaneously in the top face and in the web at its 
junct#on with the top flange. The inclination of the initial cracks 
with respect to the axis of the beam was 37, 36, 31 and 36 degrees 
in series A, B, C and D respectively. As the torsional load was in- 
creased, these cracks extended across the width of the top flange and 
down the vertical faces. At about the same load, cracks occurred in- 
dependently in the web. Web cracks extended down to the bottom of the 
web. Generally the web cracks joined the wider cracks on the flange 
sides, but in some cases they were slightly apart. In other cases 


they were bridged by a horizontal crack along the web-flange junction. 
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When the web cracks spread downwards, beams A-0-4 and C-0-4 
without web reinforcement failed suddenly. As the beam approached 
failure load one major web crack opened with the appearance of web 
crushing on the opposite face of the web in a skew bending mode. These 
beams without stirrups developed only one crack. 

In beams with stirrups further increases in torque caused 
more cracks to appear on the top face and the web. As the stirrup 
spacing reduced the cracks were more numerous. In the later stages 
of loading several diagonal tension cracks formed in the bottom flange 
Spreading continuously up the vertical sides of the bottom flange. 

As the failure was approached the diagonal tensile cracks on the top 
face, and the side faces of both flanges were quite wide, and the 
torsional displacements between concrete on either side of these 

cracks were clearly visible. Failure was triggered by yielding of 

the stirrups accompanied by compression distress in the top flange, 
together with web crushing and shear compression failure of the bottom 
flange. The shear compressive distress was more extensive in the 
bottom flange if the eccentricity was high as in the case of 

A-4-4, A-5-4, and D-3-4 or if prestress was high as in the case of 
B.3.4. If the level and eccentricity of prestress was less this dis- 
tress was also less as evidenced in the pure torsion beams of series C. 
Yielding of stirrups occurred only one or two load increments prior 

to failure. Beams in pure torsion failed in skew bending Mode 2; in 
some beams the line of web crushing was somewhat arbitrary and it was either 


on the south or on the north side of the web. Beams with stirrups ex- 
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hibited increased strength and ductility as compared to those without. 

The type of web reinforcement used in the present test beams 
has adequately resisted the torsional load without any anchorage 
failure; the vertical legs developing yield strain. This can be found 
by comparing two beams with and without stirrups. A-0-4 without stirrups 
cracked initially on the top face and the failure followed immediately. 
In the case of all beams with stirrups redistribution of stresses occurred 
after cracking; stirrups carried the tensile forces formerly held by 
the concrete. Further increases of load caused the stirrups to yield 
and the gain in strength beyond cracking was as much as 25 percent for 
pure torsion beams, depending on the amount of stirrups. 

Beams tested at load combination 3 had high torsion, low 
bending moment and shear force. The initial diagonal tension cracks 
in these beams originated on the south face of the web. On this face 
the torsional and flexural shear stresses were additive while they 
were subtractive on the north face of the web. Cracks on the north 
face were comparatively fewer in number than on the opposite face. The 
first cracks were sloped to the axis of the beam at 33, 38 and 36 de- 
grees respectively for series A, B and C. A further one or two incre- 
ments of load developed cracks on the top face which extended over 
the top flange sides. On the south face of the web numerous diagonal 
cracks appeared fairly evenly propagating throughout the web; the 
closer the spacing of the stirrups the closer the spacing of the 
cracks. The bottom flange face developed cracks independently which 


linked up to the cracks on its sides occurring at about the same time. 
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Further increase of load caused opening up of one major crack and sub- 
sequently failure resulted from yielding of the stirrups followed by 
a shear compression failure of both the flanges and compression dis- 
tress on the north face of the web in Mode 2; these were readily noticeable 
in the case of highly prestressed beams of Series B, but not clear in the 
case of other beams. These failures were highly ductile. In the range 
of large rotation, beams with stirrup spacing at 3 inches exhibited 
increasing resistance to load as the rotation increased. For beams 
with larger stirrup spacing increasing rotation resulted in decreasing 
strength as in Series A and B, and increasing strength in Series C. 
This rate of increase or decrease of strength as rotation increased 
depended on the level and eccentricity of prestress, as can be found 
inoFigse 6.) sthrougn Figs 6.0. 

Beams loaded at load combination 2 developed initial cracks 
mainly in the web and »ccasionally in the bottom flange at 83 to 88 
percent of the ultimate load. The inclination of these cracks to the 
axis of the beam was 27, 33 and 26 degrees for series A, B and C beams. 
Crack propagation was similar to that described for load combination 3. 
All beams appeared to fail in a torsion mode accompanied by a shear 
compression failure of the flanges for most beams; however the line of 
web compression was not clear in several beams. 

The beams tested at load combination 1 had a high proportion 
of flexural moment and shear. Except in series D, the T/M ratio is 
approximately equal to the ratio of the ultimate torsional and bend- 


ing capacities. The first cracks to form were flexural either on the 
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FIG. 6.1 TORQUE-MAXIMUM WEB REINFORCEMENT STRAIN 
CURVES OF BEAMS UNDER PURE TORSION 
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FIG. 6.2 TORQUE-MAXIMUM WEB REINFORCEMENT STRAIN 
CURVES OF BEAMS AT LOAD COMBINATION 3 
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STRAIN CURVES OF BEAMS AT LOAD COMBINATION 2 
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bottom or on the south vertical side of the bottom flange. Depending 
on the magnitude and level of prestress the cracking torque was 
between 50 and 73 percent of the ultimate load. The initial cracks 
sloped at 50 to 75 degrees to the axis of the beam. The cracks on 
the bottom flange sides were nearly vertical and linked up those on 
the bottom face. As the load increased cracks formed in the lower 
parts of the web and extended upwards with further increase of load, 
travelling towards the load point. Cracks on the north side of the 
web were fewer than on the south side and occasionally there were 
none. If they were present they could not propagate as far as top 
flange except in the case of Series C beams. The cracks on this side 
were inclined as on the south side, towards the flexural load because 
of the greater influence of flexural shear stresses. On the south side, 
the diagonal cracks reached the top flange-web junction and after 
travelling along this junction extended into flange side at a smaller 
inclination to the horizontal than in the web. At higher loads, 
cracks formed independently in the middle of the top face and these 
joined the cracks on the side faces. At this point the outstanding 
parts of the top flange started to split away from the web, reducing 
the available area resisting the compressive force. Consequently 
failure occurred immediately in a skew bending mode since the middle 
of top face crushed and the compression distress extended into the 
web. This was the general behavior observed for all these beams with 
slight variations among the different series. Beam A-4-1, tested at a 


Slightly higher T/M ratio,exhibited more torsional cracking. In beam 
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D-3-1 the bottom strands fractured at failure, developing nearly the 
flexural capacity and the torque associated with its T/M ratio. Irre- 
spective wai magnitude and eccentricity of prestress and variation 
of amounts cf longitudinal and stirrup einforcements, these beams 
failed approximately at the same loads. The likely reason for this 

is the reduced compressed area of the top flange due to torsional 
cracking. When the torsional cracks reached the top flange-web 
junction, these tended to split off the outstanding flanges. The 
effective compression zone then decreased significantly and the compressive 
force required to balance the yield forces of the tensile reinforce- 
ment could not be developed; hence the failure was sudden and violent. 
Such a behavior will probably be nonexistent in rectangular beams, 

and is applicable to I beams. 

Beams A-3-0, B-3-0 and D-3-0 were subjected to bending and 
shear only. D-3-0 having least longitudinal steel of the above three 
beams, failed in bending in an under-reinforced manner. The initial 
tensile cracks appeared at about 50 percent of the ultimate load 
in the bottom flange and spread to the side faces of the flange al- 
most vertically. As the load increased, these cracks extended into 
the web and propagated towards the transverse load point. They reached 
the web-top flange junction at about 90 percent of the failure load. 
At the same time tensile cracks in the bottom flange were quite wide. 
When inelastic strand strains were reached large deflections occurred 
with little increase in load, finally leading to crushing of the 


concrete. This was followed by fracture of the bottom strands. The 
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beam thus developed its full flexural capacity. 

Beam B-3-0 had a higher prestress than D-3-0 and A-3-0 and 
was slightly over-reinforced. The initial cracks originated in the 
bottom flange at approximately 50 percent of the ultimate load. The 
cracks curved towards the load point in the web as in the case of 
D-3-0, but the number of cracks were fewer in this case and extended 
to only about 75 percent of web depth at 90 percent of ultimate load. 
They did not seem to penetrate further. Failure occurred suddenly and 
violently in that the crushing extended into the web fully and it was 
completely destroyed. The strands in the compression zone buckled. 

Beam A-3-0 was designed to exhibit approximately a balanced 
failure, but failed in a slightly over-reinforced way. First cracks 
appeared in the bottom flange at about 45 percent of the ultimate 
load and the behavior was nearly same as that of B-3-0 except that 
the cracks which curved towards the load point extended higher in the 
web than in the case of B-3-0, almost reaching web-top flange junction. 
In this beam also at failure the compression zone was completely 
destroyed, the compression distress was evident in most of the web 
and the top strands also buckled. 

Generally for all beams failure occurred in skew bending 
either in a bending mode or a torsion mode. However in a few cases, 
particularly in some beams tested using load combination 2, the failure 
mode was not clear. They might have failed in a mode transitional 
between bending mode and torsion mode. In some cases it was diffi- 


cult to establish the failure section, because of discontinuous 
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cracks, and absence of a line of compression distress. 


6.3 Deformations 
6.3.1 Torsional Stiffness 

The torque-twist curves for all beams are presented in 
Appendix A. In Table 6.2 the test secant torsional stiffness at crack- 
ing is tabulated for each beam. The torque-twist curves are essentially 
linear in the initial portion up to about 70 percent of the cracking 
load. A study of the initial torsional stiffness reveals that it is 
apparently little affected by the amount of transverse reinforcement, 
loading combination, level and eccentricity of prestress and the amount 
of flexural shear. When the cracking torque is approached the torque- 
twist relationship becomes nonlinear, more so in the case of beams with 
high cracking torques. The secant torsional stiffness at cracking as 
Shown in Table 6.2 is approximately the same for all the beams except 
for those tested with loading combination 1, that is, with predominant 
bending moment where it was slightly larger. This is probably due to 
the following reason: Beams exhibiting higher torsional stiffness 
cracked in the bottom face in bending mode at considerably smaller 
torques than other beams which failed in torsion mode. At higher 
torque the torsional stiffness was smaller due to microcracking and 
smaller shear modulus. It can be seen from Table 6.2 that the beams in 
C series tested at load combination 1 and cracking at the lowest torques, 
have the highest torsional stiffnesses. 

Post cracking torsional stiffness depended on several factors, 


primarily on the level of prestress, eccentricity of prestress, T/M 
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ratio, and amount of steel. 

The torque-twist diagrams indicate that as the magnitude of 
the flexural shear increased the torsional ductility decreased. This 
can be clearly observed when beams tested at loading combination 1] and 
others are compared. A decrease in T/M ratio caused a decrease in the 
rotation capacity. When T/M ratio was low the loss in torsional stiff- 
ness after cracking occurred very slowly and gradually. This is be- 
cause the progress of flexural cracking is less extensive than that 
of torsional cracking. For a beam under pure torsion, addition of 
small amounts of flexural moment and shear delayed cracking, changed 
the mode of cracking and enabled the beam to retain the precracking 
stiffness up to torques greater than pure torque capacity. However 
subsequent decline in post cracking stiffness was faster. This be- 
havior is readily noted in B-3-4, B-3-3, and B-3-3 and applies equally 
well to other series. Beams with less eccentricity exhibited a less 
gradual reduction in torsional stiffness than beams with higher eccen- 
tricity in the vicinity of the cracking load since cracking is more 
gradual in the latter case. It can be seen in Table 6.2 and also in 
the torque-twist curves for series B beams that the angle of twist at 
failure was generally less than for corresponding beams in other series. 
Stated otherwise, prestress reduced torsional ductility. It is ob- 
served from the torque-twist curves that the stiffness of all beams 
was reduced nearly to zero at ultimate load. 

6.3.2 Flexural Stiffness 
Load deflection curves for all beams are presented in Appendix 


A. The initial portions of the load-deflection curves are generally 
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straight. The flexural stiffnesses of the uncracked beams corresponding 
to the earlier stages of the load-deflection diagrams are about the 
same. These curves are similar in shape to torque-twist curves in the 
uncracked range. After cracking, the flexural stiffness reduced, 
gradually, and also flexural ductility increased as T/M ratio decreased 


after cracking. 


6.4 Reinforcement Strains 

Data obtained from the strain gages is generally qualitative 
because of several uncertainities. However from the measured strains 
several behavioral trends can be observed. The strain gages in the 
present test program were not monitored continuously. Measurements 
were made at the end of each load increment. In the case of abrupt 
failures the readings may correspond to a load slightly less than the 
ultimate load. In a few beams the gages were intersected only by minor 
cracks but not by failure cracks. 

Prior to cracking the stirrups were virtually unstressed. 
It can be seen in Fig. 6.1 through Fig. 6.3 that the maximum stirrup 
strain at cracking was about 0.0001 for beams under pure torsion and 
about 0.0003 to 0.00035 for beams cracking in the web. For most beams 
these strains did not increase significantly immediately upon cracking. 
The strains increased rapidly as the loadings approached the ultimate 
Stage. Stirrups generally yielded when crossed by the failure crack. 
In several beams yielding of stirrups occurred only at or immediately 
after ultimate load. In other beams stirrups yielded at increased 
rotation of the beam at a load lower than the ultimate load as in the 


case of A-3-2 and D-3-2. Stirrups did not yield in B-3-2 and B-5-2 
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though in B-3-4 stirrups exhibited yield strains. In no case did 

strands yield in the beams failing in the torsion mode. Vertical 

stirrup legs in the web close to the compression zone recorded tensile 
strains. Other investigators made a similar observation in prestressed 
rectangular beams. As T/M ratio reduced, the strand strains increased 
and the top strands recorded compressive strains. In the pure torsion 
beams of Series A, the strand strains were of the order of 0.0005. 

When 5 decreased, strain in prestressing strand near the bottom in- 
creased for each additional load increments; the increase being dependent 
on the load ratio 6. 

For typical beams failing in a bending mode, the relationship 
between the transverse load and the bottom strand strain is shown in 
Fig. 6.4. Except D-3-1 which failed developing its full flexural 
capacity and torque corresponding to the T/M ratio, where the bottom 
strands fractured, no other beam developed yield strain in the strands 
because of the premature failure of the top flange. Though the cracks 
were very wide the measured induced strains over and above the pre- 
strain, were of the order of only 0.003 to 0.004 in various beams 
failing in the bending mode. This is probably due to the larger contri- 


bution of torsion to the crack width as compared to that from bending. 


6.5 RatioM /V id 
Oneal 
The shear span of 72 in. adopted in the test program shown in 
Fig. 5.5. was always the same. However, the resulting M Vue ratios 
were different for the test beams since the location of the failure 


section depended on the ratio TeV EDS or TM From Table 6.1 the 
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ratio M Vo for each beam can be computed readily. It ranges from 
2.63 to 3.96 for beams at load combination 3, 3.81 to 4.85 for all 
except one beam of load combination 2, 5.68 to 6.03 for beams at load 


combination 1, and 5.98 to 6.07 for beams tested in flexure-shear. 


6.6 Failure Section 

In tests under combined loadings giving rise to a moment 
gradient, the failure surface assumes importance, as the ultimate 
flexural moment can be established correctly only if the location of 
the failure surface is known precisely. Generally a bending mode 
failure occurs close to the transverse load, in a torsion mode the 
location of the failure depends on the T/M ratio while in the case of 
pure torsion it is arbitrary. Because it is impossible to predict 
the failure surface theoretically and no clearly defined location of 
the failure surface was observed in the tests, the test results at 
best are useful to indicate regions of failure surface location. In 
the bending mode, the failure occurred at a distance of 7 to 10 inches 
from the centre of the transverse load which is applied through a six 
inch wide bearing plate. As the flexural shear and moment were re- 
duced relative to torque the failure section moved away from the 
transverse load. The relationship between load ratio 6 and the lo- 
cation of the failure section relative to the transverse load x is 
shown in Fig. 6.5. Beams tested under load combinations 2 and 3, 
failing in torsion mode, had average A differing by about 10 inches. 
However the difference of Xi for bending and torsion mode failures 
was about 20 inches. The critical section for beams in pure torsion, 


not shown in Fig. 6.5, but listed in Table 6.1 varied from about 30 
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inches to 38 inches. if the strength prediction equations do not have 

a term involving flexural moment, the value of x is not required in 
estimating the ultimate capacity in torsion mode failures. However, 

it will be difficult to predict the flexural moment with equal precision. 
For predicting bending mode failures this location must be known before- 
hand. However a conservative result can be obtained assuming the critical 


section at the transverse load. 
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CHAPTER VII 
APPLICATION OF THEORETICAL ANALYSES 
AND DISCUSSION OF TEST RESULTS 


7.1 Introduction 

In this Chapter, the test results are evaluated; the cracking 
and ultimate strength theories developed in Chapters III and IV are 
applied to the test specimens of the present investigation and several 
others. Comparisons of observed and analytical results are made in 
Tables 7.1 through 7.7. The interaction of torsion, bending and shear 


is discussed. 


7.2 Cracking Strength 
7.2.1 Method 2 (Finite Element Method) 

First, the test results of beams of the present investiga- 
tion are compared with the theoretical cracking strengths obtained by 
the finite element method of analysis. These strengths together with 
the maximum principal tensile stresses and deformations are furnished 
in Table 7.1. For the convenience of discussion beams in which initial 
cracks occurred on the top face, in the web, and on the bottom face 
are defined as cracking in Mode 3, Mode 2 and Mode 1 respectively, 
analogous to the designations of modes of failure of beams. The test 
beams exhibited all the three modes of initial cracking. Except for 


One beam, the theory did not over predict the strength of any of the 
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CAB eical 


COMPARISON OF FINITE ELEMENT ANALYSIS WITH TEST RESULTS 


Beam fe 

No. (psi) 
A-0-4 4425 
A-5-4 5480 
A-4-4 4708 
A-3-4 4254 
A-S-3 4372 
A-4-3 5140 
A-3-3 4460 
A-5-2 4549 
A-4-2 4737 
A-3-2 4390 
A-0-1 5383 
A-5-1 4896 
A-4-] 4407 
A-3-1 5686 
A-3-0 -- 

B-3-4 4932 
B-5-3 4401 
B-4-3 5135 
8-3-3 473) 
B-5-2 5109 
B-4-2 443] 
B-3-2 4670 
B-5-1 4930 
B-4-1 4684 
B-3-1 5279 
B-3-0 -- 

C-0-4 4237 
C-5-4 4516 
C-3-4 5132 
C-5-3 4440 
C-4-3 5330 
C-3-3 5539 
C-5-2 4930 
C-4-2 4676 
C-3-2 4525 
C-5-1 4734 
C-4-1 5003 
C-3-1 4519 
D-3-4 4879 
D-3-4A 4407 
D-3-3 5238 
D-3-2 4767 
D-3-1 4955 
0-3-0 -- 


10.40 
10.40 
10.40 

5.20 
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w 


hoxinim ne ee A a ale 

Pp T.. Test tensile stress FEM T of T 

TRA i dC eae ML Oka Seals cna GUNiy ay Crum Gi) 
FEM (psi) CYCFEM) (Test) (FEM) t 1002 0 (Test) (FEM) At 100% of 
Test T.,. (FEM) Test T., (FEM) 

0 36.02 523 34.63 1.04 140 130 1.08 - 6.9 = 

0 47.94 657 41.33 1.16 165 153 1.08 : 8.3 : 

0 46.59 687 36.12 1.29 175 168 1.04 e 6.2 = 

0 38.64 559 34.50 ee 143 114 1.25 2 5.9 : 
5.10 44.21 689 35.37 1.25 100 160 .63 25 32 .78 
5.80 50.28 800 38.38 1.31 175 182 .96 20 34 .59 
5.00 43.34 723 33.60 1.29 130 157 .83 30 30 1.00 
9.00 39.01 685 31.98 1.22 152 141 1.08 43 52 .83 
8.00 38.33 669 32,21 1.19 175 138 ier, 50 45 ei 
11.00 47.68 861 33.58 1.42 190 S72 1.10 65 64 1.02 
13.50 25.04 772 21.59 1.16 65 1 71 83 7 Test 
13.50 25.08 737 21.81 1.15 74 90 .82 70 73 .96 
13.00 27.84 743 23.39 1.19 75 100 .75 83 72 1.15 
14.00 25.99 675 24.29 1.07 87 83 1.05 85 7 1.20 
=e = =e a= == oD _— ae = ~ Re = 

0 53.33 73 43.0) 1,24 224 193 1.16 = 4.3 : 
6.70 58.08 873 40.33 1.44 210 209 1.00 33 39 .85 
7.00 60.68 893 43.65 1.39 209 220 95 34 39 .87 
6.60 57.21 858 41.16 1.39 245 214 1.14 43 37 1.16 
13.25 57.43 994 39.07 1.47 216 314 69 78 72 1.08 
13.00 56.35 1082 37.07 1.52 197 211 .93 70 68 1.03 
12.50 54.18 952 39.26 1.38 175 203 .86 63 67 .94 
15.00 27.86 602 26.52 1.05 86 105 .81 80 78 1.03 
18.00 33.44 962 25.72 1.30 98 122 80 100 93 1.08 
16.00 29.72 707 26.77 Telit 101 107 94 86 82 1.05 

0 41.26 594 35.57 1.16 149 149 1.00 = 2 : 

0 45.24 681 36.19 1.25 153 169 91 = 2 : 

0 45.25 677 38.03 1.19 150 169 .89 - 2 2 
4.80 41.61 701 32.76 1.27 159 150 1.06 32 26 1.23 
5.60 48.55 869 34.68 1.40 157 182 86 35 29 eal 
5.70 49.41 880 35.55 1.39 175 185 95 25 29 .86 
8.50 39.27 793 30.92 127, 166 147 1.13 56 44 1.27 
10.00 44.75 1000 29.83 1.50 161 168 .96 60 52 1.15 
8.50 36.84 787 28.12 1.31 135 138 .98 50 43 1.16 
10.80 20.07 653 17.92 Wie 45 72 .63 63 56 Tas 
11.00 20.43 697 17.92 1.14 44 74 59 71 56 1.27 
10.50 19.48 629 17.55 wv 50 7) .70 65 55 1.18 

0 45.24 640 38.34 1.18 180 162 eal = 4 : 
2.90 50.29 757 37.53 1.34 191 182 1.05 n 20 .55 
6.60 57.27 988 37.93 1.51 210 208 1.01 40 37 1.08 
11.50 49.85 1088 32.80 1.52 179 18] 99 6] 63 .97 
10.00 18.58 450 19.85 94 57 67 85 60 54 1. 


Average 1.26 -- -- 94 - - 1.03 
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TABLE 7.2 


COMPARISON OF THEORETICAL AND TEST CRACKING STRENGTHS 


een Plastic? Plastic? 

Ream ' T ee —_—_—__—___— [es ay peo theory Top (Test) Top (Test) To, (Test) 
No. Vie (kip) Distance Location cr ere cr pee? TO (FEM) (Plastic? T (Plastic? 

a A rue (in-kip) (in-kip) (in-kip) (in-kip) Ct’ theory) Er theory) 
A-0-4 © 0 34 WS 36.02 34.63 34.84 34.84 1.04 1.03 1.03 
A-5-4 © 0 8 le 47.94 41.33 39.12 39.12 1.16 1.23 1.23 
A-4-4 © 0 14 ile 46.59 36.12 37.30 37.30 1.29 Wars} We25 
A-3-4 ©0 0 29 we 38.64 34.50 35.62 35:62 Wedi2 1.08 1.08 
A-5-3 10.40 Dail 32 W 44.21 Sy by 41.22 40.16 125 1.07 lienlet 
A-4-3 10.40 5.8 32 W 50.28 38 . 38 45.01 44.58 ibssill ol Tels 
A-3-3 10.40 5.0 28 W 43.34 33.60 39.01 Vio 1.29 Weil Uails 
A-5-2 5.20 9.0 10 W 39.01 31.98 34.93 34.20 eee Ua 1.14 
A-4-2 5.75 8.0 12 W 38.33 S221 35.65 34.32 ak) 1.08 1.12 
A= 3=Ze oc Oe lil aO 15 W 47.68 S}slents 37.98 36.77 1.42 1.26 1.30 
A-0-1 sefoe, (Sa 14 BF 25.04 21.59 21.41 Ss? Woe Trea 1.29 
A-5-1 (de (ae ACTER) 12 BF 25.08 21.81 21.91 20835 115 1.14 dees 
A-4-1 Aoi —— [Sul0) 6 BF 27.84 23.39 22.14 22.14 1.19 Te26 1.26 
A-3-1 2.23) 140 10.5 BF 25.99 24.29 22.99 21.76 1.07 des 1.19 
A-3-0 0 -- -- -- -- -- -- -- -- -- -- 
B-3-4 00 0 30 TF/W secs 43.01 47.15 Bel) 1.24 Woes} lens 
B-5-3 10.40 6.7 30/55 W 58.08 40.33 47.82 46.24 1.44 Mere 1.26 
B-4-3 10.40 7.0 33 W 60.68 43.65 52.45 51.83 1.39 1.16 Lentz 
B-3-3 10.40 6.6 17 W Bee 41.16 47.68 46.75 1.39 1.20 ez 
B=5=26 65220) 3425 28 W/BF 57.43 39.07 46.57 43.17 1.47 lees la33 
B-4-2 5.20 13.0 30 W 56.35 37.07 40.93 36.78 Thalgy4 1.38 less 
Bas=ce onc Oliaeo 21 W 54.18 39.24 41.52 39.22 1.38 1.30 1.38 
B-5-1 ZOCS melon) 10 BF 27.86 26.53 26.24 25.04 1205 1.06 1.1 
B-4-] nike. Utah) 8 BF 33.44 PAS A 24.81 24.25 1.30 les 1.38 
B-3-] ansten Hal) 7 BF 29.72 26.77 26.80 26.49 We le l6 lee 
B-3-0 0 -- -- -- -- -- -- -- -- -- -- 
C-0-4 00 0 gets TF 41.26 SN guY/ 40.95 40.95 eG: 1.01 1.01 
C-5-4 00 0 12 We 45.24 36.19 42.07 42.07 ee 1.08 1.08 
C-3-4 oo 0 a2 TF 45.25 38.03 43.82 43.82 sus 1.03 1.03 
C-5-3 10.4 4.8 34 W 41.6) San hs 38.20 34.59 leew 1.09 1.20 
C-4-3 10.4 5.6 36 W 48.55 34.68 41.13 36.83 1.40 1.18 Tes2 
C-3-3 10.4 aif 40 W 49.4) 3On05 42.46 37.95 1.39 1.16 1.30 
C-5-2 5.54 8.5 S0k5 W 39.27 30.92 34.74 29.17 Weey Vas 35 
C-4-2 5.20 10.0 17&31 TF to BF 44.75 29.83 33.84 27.96 1.50 ese 1.60 
C-3-2 5.20 oe 6 BF 36.84 28.12 26.08 26.08 Test 1.4] 1.41 
C-5-1 Aasy | |to)stel lst) BF 20.07 Wipes 16.78 UteTAS Holl? 20 era 
C-4-] Adee Niel) 13.75 BF 20.43 17.92 17.24 15.61 1.14 Ueuks 1.31 
C-3-1 Aavasy — (alets 250: BF 19.48 UWS SE: 16.62 15.30 (lai AUT Vee? 
6-3-9 $96, 0 ay TF 45.24 38.34 41.08 41.08 1.18 1.10 1.10 
D-3-4A 20.8 2.9 40 TF 50.29 3/553 43.21 44.68 1.34 Wo lla lets 
0-3-3 10.4 6.6 29 TF/W YL AC4ll 37.93 44.89 42.60 Qheoit 1.28 1.34 
D-3-2 5.2 Titi) 43.5/6 W/BF 49.85 32.80 38.69 32.00 oe Wer4s) 1.56 
D-3-1 22230 10R0 yak) BF 18.58 19.85 17.76 igets¥/ 0.94 1.05 1.04 
3-0 =- oe ae == aa = as e- ie ~ os 
TF - Top Face BF - Bottom Face W - Web Average ieee) etalts es: 
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ie 
beams. The beams under pure torsion cracked in Mode 3. In this case, 
the maximum theoretical principal tensile stresses varied from about 
550 psi to 700 psi in all the series. When the flexural moment and 
shear were predominant, initial cracking occurred in Mode 1. The maxi- 
mum tensile stresses corresponding to this mode ranged from about 600 
psi to 770 psi. The correlation between the predicted and the test 
strengths is excellent for these beams exhibiting initial cracking in 
Mode 3 and Mode 1 as is evidenced by the mean of the ratios of test to 
theoretical strengths of 1.18 and 1.12 in the former and latter modes 
of cracking. 

For Mode 2 cracking which occurred in the thin web, the theory 
proved to be more conservative than in the other cases. The theoretical 
maximum principal tensile stresses were in the range of 700 psi to 
1100 psi for load combinations 2 and 3. Beams tested at load com- 
bination 2 exhibited higher stresses than those at load combination 3. 
The test strengths are in excess of predicted strengths by 36 percent 
for load combination 3 and 38 percent for the load combination 2, on 
the average. However the maximum excess of strength in a few cases 
was of the order of 40 to 50 percent. This correlation can be con- 
Sidered only fair. Similar trends were also clearly evident in the 
case of web cracking when the strengths were predicted by an alternate 
method, that is by Method 1, which will be discussed later. While 
studying the interaction of torque and shear at cracking, Wyss and 
Mattock (1971) observed test strengths to the extent of about 50 


percent larger than the calculated ones, when the girders cracked 
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initially in the web. This is not in a way unexpected due to the follow- 
ing explanation. There is a high stress gradient across the thin web 
of the I section subjected to either pure torsion or combined loadings. 
It is known that concrete cracks at higher apparent tensile strengths 
when subjected to higher stress gradients. Further if a plastic dis- 
tribution of shear stress is assumed, the entire web will be stressed 
to the same level under pure torsion. Superposition of flexural shear 
stresses on this torsional shear stress distribution will reduce the 
tensile stresses on one half of the web to zero if shear stresses 

due to both forces are equal. Cracking would then occur at higher 
apparent stresses in a web stressed to the same maximum level as 

in pure torsion cases or pure shear cases than in a web one half 

of which is less stressed than the other, as in the case in combined 


loadings. 


7.2.2 Method 2 - Prestressed Rectangular Beams 

Several rectangular prestressed beams in the literature were 
analyzed for cracking strength by the finite element method and the 
results are compared with the experimental values in Tables 7.3 and 
7.4. The analysis was applied in each investigation only to re- 
presentative beams of various loading ratios, eccentricities and levels 
of prestress for reasons of computer time. 

Specimens of the test program of Mukherjee and Warwaruk (1970) 
exhibited the three modes of cracking. The maximum principal tensile 
stresses varied from about 500 psi to 1200 psi leaving one beam V205 


which behaved differently from its group. In the case of Mode 1 cracking 
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COMPARISON OF FINITE ELEMENT ANALYSIS 
WITH TEST RESUDMS 


Max imum 
Effective principal T T 
Beam f! ae tensile stress cr cr cr Location 
Investigator Nes c  Prestress Eccentricity T/M T/Vb ae ee (Test) ( FEM) (Test) of crack 
(psi)  (kips) (in) ae eS 
ing loads by (in-kip)  (in-kip) cr rem) (FEM) 
FEM (psi) 
122 5260 Bones 2 1/3 -- --- 70.0 60.3 1.16 Bot tom 
126 5462 36.37 2 © -- 689 104.0 87.4 1.19 Long side 
127 +5014 38.72 2 1 7/ -- 735 36.0 31..:6 1.14 Bottom 
202 4469 86.54 0 1/3 -- 1196 100.0 65.4 Tae Bottom 
pike 206 5686 84.48 0 0 -- 681 13320 zie 1.13 Long side 
= 207 5005 91.60 0 We -- 1179 50.0 36.8 12365 Bottom 
dees ZEL OSSS 84.53 2 1/3 -- 835 119.0 101.7 ibn by Bottom 
2 226 5633 84 .37 Z 9 -- 821 130.0 89.7 1.45 Top 
z 227 +4778 90.75 2 1/7 -- 866 65.5 570 1,15 Bottom 
= 
' V102 3988 37.78 0 HAS! W267, 643 30.3 Nays gilts) Bottom 
cy V105 4118 eSicil 0 3 19.40 499 78.0 Tord: 1208 Long side 
S. 
& V107 3685 36.04 0 Wie 2927. 703 13.9 AeO 1.26 Bottom 
at 
= V122 4268 32.34 2 1S) Zor 74) 39.0 34.8 isi Bottom 
v125 4870 35,65 2 3 19.50 616 93.6 82.) 1.14 Long side 
V127 4619 34.12 2 Wie WSS 67] 18.6 16.8 1.11 Bottom 
v202 5000 7857 0 Sh osallay/ 480 39.0 40.6 -96 Bottom 
V205 5618 84.8] 0 3 19.50 35 70.2 100.3 .70 Long side 
V207 5041 81.53 0 1/7 .928 926 24.1 19.4 1.24 Bottom 
V222 5766 80.60 2 Sian nO 64] 69.4 66.7 1.04 Bottom 
V225 5389 80.04 2 3 19.59 515 85.8 9i/..3 .94 Top 
Average W158! 
Oe gee 2-0 5620 60.98 0 » -- 660 ib bile 7A} 100.63 Neat Long side 
gus 2-1 5620 60.12 0 1.0 9 == 768 112.98 92.61 1.22 Bottom 
EBD 2-8 5730 59.83 0 We Se 996 32.20 24.58 1.31 Bottom 
WV-1-3 7230 58.0 0 56 1.47 988 70.0 54.7 1.28 Bottom 
IV-2-3 7230 58.0 0 42 eel 785 50.0 45.0 ela Bottom 
IV-3-3 7230 58.0 0 ve Ve965 1085 50.0 37.6 1.33 Bottom 
IV-4-3 7230 58.0 0 12 63 908 70.0 25.0 1.20 Bottom 
V-1-3 4900 75.52 74 als oaks 749 110.0 86.6 Wor Long side 
V-2-3 4900 py .74 walls 50 768 30.0 26.3 1.14 Bottom 
V-3-3 4900 Usyoey? .74 225 5/45 1048 50.0 38.2 leo Bottom 
J V-4-3 4900 7552 .74 Se haze! 786 70.0 59.3 1.18 Bottom 
fon) 
= IV-1-4 5900 56.6 0 wien peed) 784 70.0 58.8 sts Bottom 
= I1V-2-4 5900 56.6 0 29° Wes 1020 50.0 373 1.34 Bottom 
$e IV-3-4 5900 56.6 0 ee thesis) 754 40.0 3551 1.14 Bottom 
eS IV-4-4 5900 56.6 0 wor 05 62] 30.0 29.1 103 Bottom 
Ss 
s V-1-4 4480 elsif) .74 2.44 4.785 598 100.0 88.5 Veils} Long side 
V-2-4 4480 73575 .74 A wrerele 964 95.0 70.4 esd Bottom 
V-3-4 4480 T3075 74 sth hee 973 60.0 47.2 lisey Bottom 
V-4-4 4480 USMS .74 16 2.745 1147 40.0 29.6 Teas Bottom 
Average 7.22 
: eee 
Bee A3 6000 226 4.9 aes 677 300 260 eis Top 
foal 
s38 B6 6000 104 5.0 -- --- 536 200 184.4 .92 Web 
Ss eo c4 6180 0 0 -- --- 549 164 Sses) 5} Web 


Average 1.00 
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Mukherjee - Warwaruk 


Beam 


22 
126 
127 
202 
206 
207 
222 
226 
227 
v102 
V105 
V107 
v122 
V125 
V127 
v202 
V205 
V207 
V222 
V225 
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COMPARISON OF FINITE ELEMENT ANALYSIS 


WITH EXPERIMENTAL DEFORMATIONS 


(Rad per in x 10-6) 
At 70% of At 100% of At 70% of Test 
test T. test T.. Top (Test) 
(Test) (FEM) At 100% of Test 
rep (FEM) 

5] 59 0.93 
88 82 lO, 

24 28 0.86 

79 79 1.00 

92 105 0.88 

27 39 0.77 

78 94 0.83 

97 103 0.94 

38 51 0.76 

25 26 0.96 

69 66 1.05 

6 1] On55 

26 33 0.79 
85 90 0.94 

1] 15 0.73 

15 33 0.45 

58 59 0.98 

V7 20 0.85 
47 59 0.80 
69 33 0.95 


Angle of Twist 


Flexural Deflection 


(in x 10-3) 


Ws 


At 70% of At 100% of At 70% of Test 


test ee 
(Test) 


70.0 
100.0 
110.0 
120.0 
105.0 
165.0 

45.0 

5.0 
40.0 
65.0 
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50.0 

5210 
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a0 


test LAA 


(FEM) 


74, 


16. 


{8 GOe COINS) (eso Ge Co 


55 


Tp (Test) 
At 100% of Test 
(FEM 


(lf 


' 
' (ee) 
= 


cS) (ay =) 


[SY Ke) (Ss te) SS erie (Sy Ks! 


~™ 6 6G) OO OO UNS Sm COO Go to * 
oO ODN A TN MO W DO Nw 


Ww oO 
(Sy 


} 


Le 
ws Waa bes | ; 


7 


% 


aie ot 
; an 
_ 


F Red 
ae os oe rae - 


28.0 
8.0 
2.0 
- 
- a 7 
7) ee oo 


. 


2538 
cattery 


2 
7 


SBR" RERBAG 
ebbizuee 


> 


RES RRS SRSER TST AEsa ws | 


e235. 
92023 
_ 


LAs) 
the values of Test TU ,/ Theoretical Lan varied from 0.96 to 1.36 except 
for one beam for which it was 1.53 and the mean value was 1.18. The 
higher test to theoretical ratios occurred in beams having a uniform 
stress due to prestress. This is probably due to this: two beams one 
having uniform and the other triangularly eccentric prestress if acted 
on by a particular bending moment will yield finally triangular and 
uniform stress over the depth;in the former there is a high stress 
gradient, while in the latter the stress if uniform. If torsional 
stresses are superimposed on such triangular and uniform stress dis- 
tributions, cracking would occur at apparently higher stresses in the 
beam having a uniform stress due to prestress only and hence higher 
test strengths. 

For beams cracking in Mode 2, the ratio Test TU ,/Theoretical ie 
was 1.12 disregarding beam V205 as mentioned earlier. 

Similarly the strength ratio was 1.45 for beam 226 tested 
under pure torsion and the theoretical cracking occurred on the top 
face; but it was reported that the cracking commenced near the upper part 
of the side face. If the theoretical stress is used according to the 
test location of the crack on the side, the ratio 1.45 comes closer 
to unity. The same explanation offered for Mode 1 cracking explains 
the higher observed strength in Mode 3 cracking of eccentrically 
prestressed beam 226 subjected to pure torsion. Of the 20 beams of 
Mukherjee and Warwaruk analyzed herein, 3 beams were subjected to pure 
torsion, 6 beams to combined bending and torsion, while the remaining 


had additionally shear force. On the average the strength was con- 
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servatively predicted within 15 percent of the test strength. 

Three rectangular prestressed beams of Gangarao and Zia 
(1970) subjected to combined bending and torsion, and 16 rectangular 
prestressed beams of Henry and Zia (1971) were also analysed and the 
results furnished in Table 7.3. Beams cracking in Mode 1 developed 
principal tensile stresses ranging from 750 psi to about 1100 psi. 
These stresses are the maximum stresses given by the finite element 
analysis and correspond to the location of initial cracking. However 
the torque-moment ratio used in the cracking analysis by Henry and Zia 
related to the failure section rather than the first crack section. 
This as pointed out earlier in Chapter II, may result in computing a 
cracking strength at a section different from where the crack had actually 
occurred. If their method is adopted slightly lesser theoretical 
stresses would occur resulting in better correlation than that reported 
here. The ratio of observed to theovetical strengths varied from 1.03 
to 1.35 and the mean is 1.23. Of these 19 beams only in 3 beams 
initial cracks occurred on the vertical side at theoretical stresses 
ranging from 600 psi to 750 psi. Theory predicted the cracking 
torque 17 percent in excess of the observed torque showing good agree- 
ment for Mode 2. 
7.2.3 Method 2 - Prestressed I Beams 

Three I girders of Wyss et al (1969) under pure torsion were 
analyzed by this method. The maximum principal tensile stresses in 
these cases ranged from 536 psi to 677 psi and the ratios of observed 


to analytical strengths varied from 0.92 to 1.15 indicating very good 
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correlation. 
7.2.4 Method 1 

Experimental cracking strengths of the I beams tested in the 
present program are compared with the analytical results obtained by 
Method 1 in Table 7.2. In most beams, the centre of the ultimate failure 
surface and the point of initial cracking were not the same. In some 
cases they were close to each other. It is not possible to predict 
precisely either of these locations. In computing the cracking strengths 
by Method 1, two procedures were adopted. The only difference between 
these procedures is in fixing the location of the cracking section. 
If the moment in the test zone is uniform, this difference vanishes. 
¥-c.4.1 ‘Procedure’ | 

The w and 6 ratios corresponding to the test cracking 
section were used. The 6 values are listed in Table 7.2 and do not 
include the dead load shear. At this section, on bottom face for Mode 
1, on top face for Mode 3, and at both the web-flange junctions and the 
center of web for Mode 2 the cracking strengths were calculated. The 
smallest of these values is taken as the cracking load. The average 
Tacios Of test to theory are Vil, i221 and 1218 for Mode"3,"Mode* 1, 
and Mode 2 respectively and the combined average is lit eee his hase 
very good correlation. The results are slightly conservative for 
Mode 2 and Mode 1. Mode 1 cracking is usually associated with modulus 
of rupture, which is normally more than 6 Yfas the tensile strength 
used in this method. This is probably the reason for lower predictions 


in Mode 1. The cracking modes predicted by the theory agreed with the 
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observed cracking modes in virtually all of the test beams. 
7.2.4.2 Procedure 2 

In this procedure the critical crack is considered to occur 
on a cross section at a distance equal to half the effective depth of 
the section from the transverse load measured in the direction of de- 
creasing moment without any regard for the experimental values. Both 
procedures then would result only in negligible differences in strengths 
for Mode 1; similar is the case for Mode 2 if the cracking is assumed 
to commence at the centroid of the web. As initial cracking is assumed 
to occur anywhere in the web Procedure 2 may yield lower strengths than 
Procedure 1 in many cases. In the case of web cracking for load com- 
bination 2, Procedure 2 more underpredicted the strengths than Procedure 
leet cihis»is-evident, fromscolunns 8. 9.5.11 and 12 of. Tabley/ 22. »The 
reason is the variation of bending moment in the test span. If the 
possibility of web cracking is considered only at the neutral axis as 
was done by Wyss and Mattock (1971), the theoretical results for Mode 2 
presented in Table 7.2 came closer to the observed strengths. The average 
Otyy7at10S.0f Trey test/TO. calculated 71s. l.o22 shor.all pthe beams... 101s 
ratio is higher as explained previously than that for Procedure 1. 
7.2.5 Deformations 

No useful method is available for estimating the post crack- 
ing stiffness of prestressed concrete beams. Several investigators ob- 
served the initial portion of the torque-twist relationship to be 


essentially linear up to the commencement of cracking implying elastic 
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behavior. In this range torsional stiffness was computed using elastic 
methods and compared to the test stiffness. However differences exist 
in selecting the experimental stiffness. Henry and Zia (1971) used 
the torque on the torque-twist curve where the curve deviated from a 
straight line and the corresponding rotation, Chander et al (1970) corre- 
lated the secant modulus of torsional stiffness at 50 percent of the 
cracking torque to the stiffness derived by the elastic method for 
prestressed beams subjected to pure torsion. Woodhead and McMullen 
(1972) used calculated cracking torque but not the test torque which 
was usually higher and the corresponding experimental angle of twist 
for prestressed rectangular beams under combined loadings. Johnston 
and Zia (1971) chose the stiffness of the initial portion of the torque- 
twist curve, having observed slight nonlinearity beyond 60 to 70 
percent of test ee in the case of hollow prestressed beams. Wyss, 
Garland and Mattock (1969) took the average slope of the initial linear 
part of the torque-twist relationship for prestressed I beams under 
pure torsion. Wyss and Mattock (1971) made comparison of test and 
theoretical precracking stiffnesses for prestressed I beams under com- 
bined loadings with excellent correlation. But the way they defined 
the test stiffness was not given, presumably it is the average slope 
of the initial portion of the curve. However it is generally agreed 
that as the cracking load is approached inelastic deformations do occur, 
often to a significant extent. 

In the present study, the finite element method is used to 


predict the torsional and flexural deformations. This is an elastic 
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method but yields as a rule lesser deformations than the exact theory. 
The theoretical and experimental angles of twist at test cracking loads 
and 70 percent of test cracking loads, respectively, are compared in 
Table 7.1 for the I section beams of the present test program. The 
ratio of test to theoretical angles of twist varied from 0.6 to 1.25; 
only 5, out of 41 beams exhibited ratios close to these extreme values. 
On the average this ratio is 0.94 which indicates a very good correla- 
tion. 

Flexural deflections are also compared at the same torques as 
the torsional deformations and are listed in Table 7.1. Measured de- 
flections included also those due to prestress. However the deflections 
obtained from the finite element analysis do not include them. Hence 
the deflections due to prestress were calculated using the elastic 
method and added to those obtained by the finite element analysis. 

On the average the predicted values are within 3 percent of the observed 
values. 

Since for the prestressed concrete beams of Mukherjee and 
Warwaruk (1970) deformations in the whole range of loading are furnished, 
they are computed and listed in Table 7.4 at about 70 percent of test 
cracking strengths under combined loads. When compared with the finite 
element results, the rotations are conservatively predicted, varia- 
tions are less than in the case of I beams and except for two beams 


the correlation in each case is very good. 
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7.3 Ultimate Strength 


At any given load level the test beams were subjected to a con- 
Stant torque and shear while the bending moment varied in the test zone. 
As the torque was related to flexural shear and moment by the ratios 6 
and respectively, once torque was known other forces were obtained. 
The ultimate strength analysis developed in Chapter IV based on the ob- 
served two modes of skew bending failure, is applied to the beams of pre- 
sent investigation and several I beams of other investigations; the re- 
Sults of the former are listed in Table 7.5 together with the governing 


mode of failure. 


7.3.1 Beams Under Pure Torsion 

For all the beams tested under pure torsion, though initial 
cracking occurred in Mode 3,failure took place in Mode 2. Beams A-0-4 
and C-0-4 not provided with stirrups failed after one or two incre- 
ments of load beyond cracking. The cracking strength was 92 to 96 per- 
cent of the failure torque. Little ductility was exhibited by these 
beams. The theory cannot be applied to beams without stirrups. For 
such beams the cracking strength is taken as the ultimate strength. 

The beams were not provided with any longitudinal mild steel. 
For these beams to develop the torsional capacity corresponding to the 
volume of the web reinforcement used, the prestressing strands supplied 
the necessary longitudinal force that would otherwise have been pro- 
vided by the longitudinal reinforcement. The minimum amount of web 
reinforcement in the beams was slightly more than the minimum defined 
by the following equation suggested by Hsu (1968d) for reinforced con- 


crete beams 
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TABLES /25 
COMPARISON OF TEST RESULTS WITH THEORY - 
ULTIMATE TORQUE 


Beam Tu/Mu TusYubv Tu test  "u theory ‘utest Observed 
No (y) (4) (in-kips) (in-kips) Ty Theory mode of 
failure 
A-0-4 »& co 37533 34.84 1.07 Mode 2 
A-5-4 & © 54.95 41.84 1634 Mode 2 
A-4-4 & oo 56.03 43.29 1.29 Mode 2 
A-3-4. & oo 50.64 45.20 hel Mode 2 
A-5-3 0.849 9.811 50.28 36.93 1.36 Mode 2 
A-4-3 0.554 10.153 56.35 42.2 1.34 Mode 2 
A-3-3 0.778 9.933 56.37 40.5 1.39 Mode 2 
A-5-2 0.261 5.036 46.58 36.6 We27 Mode 2 
A-4-2 0.272 5.626 46.27 39.3 1.18 Mode 2 
A-3-2 0.243 5.189 58.5] 46.4 1.26 Mode 2 
A-0-1 0.085 22218 36.23 41.80 0.87 Mode 1 
A-5-1 0.089 2.218 44.58 44.92 0.99 Mode ] 
A-4-] 0.094 2.437 39.90 41.52 0.96 Mode 1 
A-3-] 0.087 2.220 43.90 45.83 0.96 Mode 1 
A-3-0 0.0 0.0 0.0 0.0 - No torsion 
B-3-4 »& © 61.90 57.39 1.08 Mode 2 
B-5-3 0.590 10.400 6155 44.9 Wee vs Mode 2 
B-4-3 0.628 10.135 65.88 5250 N27 Mode 2 
B-3-3 0.590 10.290 65.60 5235 1.25 Mode 2 
B-5-2 0.274 Solas 62.74 46.3 1-25 Mode 2 
B-4-2 0.269 5.050 62.75 46.3 1.35 Mode 2 
B-3-2 0.281 5.161 65.03 51.6 1.26 - Mode 2 
B-5-] 0.087 2.204 37-89 * 44.64 0.85 Mode 1 
B-4-] 0.087 2.219 43.65 43.53 1.00 Mode 1 
B-3-] 0.092 ZI329 43.89 47.04 0.93 Mode 1 
B-3-0 0.0 0.0 0.0 0.0 - No torsion 
C-0-4 @& © 44.7] 40.95 1.09 Mode 2 
C-5-4 & oo 49.56 37.250 ese Mode 2 
C-3-4 = © 58.73 45.52 1.29 Mode 2 
C-5-3 0.663 9.960 46.81 35.16 isch Mode 2 
C-4-3 0.801 9.909 55.49 Sor 1.42 Mode 2 
C-3-3 0.482 210.360 56.46 42.6 TESZ Mode 2 
C-5-2 0.596 5.328 41.88 shiei/ (Ll Mode 2 
C-4-2 0.305 Jellies) 47 .68 Sieh 1.26 Mode 2 
C-3-2 0.297 5.124 46.38 40.0 1.16 Mode 2 
C-5-1] 0.085 aces 40.50 40.19 1.01 Mode 1 
C-4-] 0.085 2.236 39.58 = (AOR) 0.98 Mode 1 
C-3-] 0.086 2.214 39.57 39.74 1.00 Mode 1 
0-3-4 & ce) 56.57 45.47 1.24 Mode 2 
D-3-4A 1.317 19.597 59.77 42.9 1.39 Mode 2 
D-3-3 0.607 9.790 62.90 42.9 1.47 Mode 2 
D-3-2 0.278 5.120 60.24 39.8 Sei Mode 2 
D-3-) 0.083 2.214 41.78 33.82 1.24 Mode 1 
0-3-0 0.0 0.0 0.0 0.0 - No torsion 
Averages 


Mode 2, Pure torsion = 1.20 
Mode 2, Combined loading = 1.317 
Mode 1] = 0.98 
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The test specimens had varying amounts of transverse rein- 
forcement. The reinforcement factor, XV Agtey/s ranged from 4.88 to 
8.13. This is a rather narrow range. The best way to use a given 
amount of reinforcement is to have a maximum lever arm. In the thin 
web of I specimens, the available lever arm between the two legs of the 
stirrups is only 1.25 inches, about one quarter of that in a 6 inch wide 
rectangular beam. 

The addition of stirrups increased the strength and ductility 
Significantly; the strength increase corresponding to the volume of stir- 
rups used in the beam. Of the three amounts of stirrups used even the 
smallest resulted in vary large ductility. An increase in the level 
of prestress increased the cracking strength but the post cracking 
strength was linearly proportional to the web reinforcement factor as 
suggested by Hsu and Kemp (1969) for prestressed concrete members. 

The higher observed strengths for beams A-5-4 and A-4-4 than that of 

A-3-4 which had more transverse reinforcement are due to higher unintended 
prestress. However the difference between the ultimate and cracking 
strengths was maximum for beam A-3-4 having the largest amount of web 
reinforcement of all the three beams. Increase in the amount of web 
reinforcement resulted in increases of strength but this was realized 

only at very large rotations where stiffness was almost negligible. 

It then becomes doubtful whether the increase in the ultimate strength 


obtained by using more transverse reinforcement is useful. Comparing 
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D-3-4, C-3-4 and A-3-4, it can be observed that beams of different ec- 
centricities developed the same post cracking strength. In other words 
these beams differing nominally only in the eccentricity exhibited the 
Same pure torsional strengths. 

As the prestress increased, the contribution from the web 
reinforcement reduced slightly as can be seen from the comparison of 
B-3-4 with A-3-4, C-3-4 and D-3-4. This means the web steel did not 
yield and an amount of stirrups smaller than that used in B-3-4 might 
have rendered the beam under reinforced transversely. The strain gage 
readings showed that the stirrups yielded. It is however difficult to 
draw firm conclusions based on the strain gage data, as much depends 
on the location of the gage relative to the crack. 

The theory predicted conservatively the test strengths of beams 
under pure torsion, failing in Mode 2. The ratios of test to theoretical 
strengths varied from 1.08 to 1.32 and the average is 1.20. The reasons 
for the conservatism are discussed in Section 7.3.3. 

Slightly better correlation would result by using the modi- 
fication (equation 2.36) suggested by Wyss et al (1969) to equation (2.35) 
proposed by Hsu (1967) applicable to concentrically prestressed beams 
subjected to pure torsion. However it predicted an unconservative re- 
sult for B-3-4. 

7.3.2 Influence of Transverse Shear Force 

In the equations derived for the torsion mode of failure, 

the shear force is present and its effect is to reduce the ultimate 


strength in combined loadings. This was also noted by Henry and Zia 
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(1971) in the case of prestressed rectangular beams. As the effect of 
bending moment is not taken into account for Mode 2, if the magnitude 
of the transverse shear were to be constant, the torsional strength also 
would remain constant. If the shear force increases the torsional 
strength decreases. However it will be clear from the interaction 
diagrams of test results that small amounts of shear force helped in 
increasing the test strength in combined loadings. This increase can- 
not be predicted by the theory. 

For Mode 1 failure, the developed equations do not involve 
shear force and can be applied only to beams which do not develop 
premature shear failure. None of the specimens tested in the present 
program failed prematurely in shear. 

At cracking, the effect of shear is to increase the tensile 
stresses on one vertical side and thus reduce the cracking loads. 
Moderate amounts of bending moment and shear have increased the cracking 
Strength, but further increases of both decreased the strength. This 
can be observed from the interaction diagrams at cracking presented in 
ek hepa aall 
7.3.3 Beams Failing in Mode 2 Under Combined Loading 

Beams tested at load combination 3 and 2 generally failed in 
a torsion mode. The theory assumes that Mode 2 failures are independent 
of flexural moment. The theoretical predictions for these beams were 
very conservative under predicting the strengths by 32 percent on the 
average. Several reasons may be responsible for such low theoretical 


strengths. 
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(a) Small amounts of bending moment and shear increased the test 
strengths in Mode 2. This increase was not reflected in the theory. 
Similar increases were observed in prestressed rectangular beams by 
other investigators. 

(b) As part of the shear compression forces were resisted by the 
outstanding flanges on the compression side, assumption of a compression 
hinge along the web might have resulted in a smaller lever arm of shear 
compression forces. 

(c) More importantly, the strain compatability factor used to 
obtain the induced strain in the prestressing strands from the yield 
Strain of stirrups most likely underestimated the induced longitudinal 
Strains. 

According to the strain compatability factor used the strain 
in the strands can assume a maximum value of the yield strain of 
stirrups, i.e. 0.0013, if the crack angle is 45 degrees. As this angle 
is far less than 45 degrees, the induced strain in the strands is less 
than 0.0013. However strain measurements on the strands have indicated 
induced strains above 0.002. Hence the relationship between the yield- 
ing and non yielding steels needs to be established more precisely in 
partially over-reinforced beams such as tested in the present program. 
The assumption used in the theory regarding the compatability condition 
iS appropriate at cracking only. A more reasonable assumption is 
warranted but cannot be made at the present time. Until then signifi- 
cant progress can not be made in the skew bending method of analysis 


for partially over-reinforced beams failing in Mode 2. The theoretical 
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predictions obtained by the present theory constitute a lower bound 
solution to strength of beams developing a Mode 2 failure. 

Similar conservatism was observed in prestressed rectangular 
beams using the same compatability of strains. Good correlation was 
obtained for the prestressed beams tested by Henry and Zia (1971) 
failing in Mode 2, only when they over estimated the forces in the 
longitudinal mild steel violating their assumption of compatability of 
strains between stirrups and the longitudinal steel. They assumed yield 
force in the longitudinal steel. But according to the compatability 
conditions the force is far less than that at yield. In analyzing 
several prestressed rectangular beams in the literature failing in 
Mode 2, and partially over-reinforced, Woodhead and McMullen (1972) 
reported test to theoretical strength ratios of 1.30 for majority of 
beams, more than 1.30 for some, and an average of 1.174. 

7.3.4 Beams Failing in Mode 1 

Beams of load combination 1 were tested under high bending 
moment; torque to bending moment ratio equalled approximately the ratio 
of pure torsional strength to flexure-shear strength in Series A, B 
and C. These beams failed in Mode 1; but as described in Chapter VI, 
the failure was premature due to the splitting off of the oustanding 
flanges. This occurs when the magnitude of torsion is sufficient 
enough to cause cracking in the top flange which is under compression. 
In such a case only the web portion of the top flange is effective in 
resisting the shear-compression forces. However the full width of the 


top flange can be completely effective in the following two cases: 
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(a) If the magnitude of torsion is small enough torsional 
cracking cannot destroy the top flange. 

(b) If by the time the initial torsional cracking commences on the 
top flange, the beam attains its flexural capacity as that in D-3-1. 

For simplicity and conservative prediction it is proposed 

to neglect the outstanding flanges in isolated girders having reinforce- 
ment detail similar to that used in the present test program. Instead of 
full flange width, as only the web portion is considered, the strength 
of D-3-1 is underestimated giving T ,test/T, theory = 1.24. If full flange 
width is considered this ratio is 1.13. For Mode 1 failures the theory 
predicted the test strengths closely. The ratios of test to theoretical 
loads varied from 0.87 to 1.24 with the exception of one beam and the 
average is 0.98. Beams in Series A, C and B at load combination 1 were 
expected to fail in an under-reinforced manner in Mode 1 as they were 
balanced or slightly over-reinforced when transverse load alone was 
acting. Beams having varying percentages of transverse reinforcement 
failed nearly at the same load because of reduced effectiveness of the 


top flange at ultimate. 


7.4 Interaction Behavior 
7.4.1 General 

The net effect of a simultaneous application of torsion, 
bending and shear on the strength of a member can be expressed by an 
interaction surface in a three dimensional rectangular coordinate system. 
Each axis then represents a type of force. Each point on the surface 
represents the strength of a beam subjected to a certain amount of 


torsion, bending and shear. It is more convenient to represent the 
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interaction of forces on two dimensional diagrams. The torsion-bending 
interaction curve can be obtained by projecting onto the torsion bend- 
ing plane the intersection of the interaction surface and a plane de- 
fined by the loading ratio of torsion and bending. Similarly the 
interaction curves on the other two planes can be obtained. To con- 
struct a three dimensional interaction surface, the interaction curves 
on three planes must be determined. It is possible to conduct tests 
and obtain these curves on torsion-bending and bending-shear planes, but 
it is difficult to get such a curve on the torque-shear plane as 
shear force cannot be easily simulated without the presence of bending 
moment in a finite length of a beam. 

In defining the pure shear strength Wie Mukherjee and Warwaruk 
(1970) and Henry and Zia (1971) adopted the equations of ACI 318-63. 
These equations are usually conservative. Hence a theoretical curve 
based on these equations may seem too conservative near the shear axis 
as beams may have test strengths in excess of theoretical predictions 
of Nene Other investigators, Ersoy and Ferguson (1968), Mattock (1968), 
Woodhead and McMullen (1972) defined Ate as the shear strength corresponding 
to the ultimate flexural capacity. This definition suffers from the 
disadvantage that Life depends on the bending moment in that geometry 
and loading patterns control the strength. Further, it may be uncon- 
servative near the shear axis if a premature shear failure occurs. 
In the present investigation, the latter definition is used in present- 


ing the two dimensional interaction diagrams. 
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7.4.2 Interaction of Torsion, Bending and Shear at Cracking 

The results of interaction between torsion and bending at 
cracking are presented in Fig. 7.1. It can be observed that the inter- 
action behavior is little affected by the amount of web reinforcement 
and that except for beams with stirrup spacing at 5 inches in Series 
A and C, small amounts of bending moment helped in increasing the tor- 
sional cracking strength. The reason for the absence of such an in- 
crease for A-5-3 and C-5-3 is the smaller prestress for them than that 
in corresponding beams A-5-4 and C-5-4. The maximum increase was about 
12 percent. The cracking torque was greatest when the applied bend- 
ing moment was equal to the internal moment created by the eccentric 
prestress. In other words a beam exhibited highest cracking strength 
when the normal stress due to prestress and bending was uniform in all 
the series. The torsional strength at cracking did not drop below the 
pure torsional cracking strength until the coexisting moment was about 
60 to 70 percent of bending shear capacity at cracking. As expected 
increases in level of prestress increased the cracking loads for all 
load combinations. The beams did not have tension due to prestress alone. 

A theoretical cracking curve is plotted for each series in 
accordance with Procedure 2 of Method 1 cracking analysis. If cracking 
initiates in the web it is assumed to occur at the centroid of the web. 
For Mode 3, initial crack is assumed to occur on the face of the top flange 
ata distance of half the flange width from the edge of the supporting plate 
of the restraining head, i.e., 8 in. from the center of support. 
7.4.3 Interaction of Torsion, Bending and Shear at Ultimate 

The non-dimensionalized interaction curves at ultimate between 


torsion and shear, and torsion and bending are shown in Fig. 7.2 through 
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(ALL BEAMS WITH STTIRRUPS AT 3" CENTERS ) 
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FIG. 7.2 NON-DIMENSTIONAL INTERACTION CURVES - TEST RESULTS 
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FIG. 7.3. THEORETICAL NON-DIMENSIONAL INTERACTION CURVES AND TEST RESULTS 
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O BEAMS WITH STIRRUPS AT 3” CENTERS 
4 BEAMS WITH STIRRUPS AT 4" CENTERS 
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FIG. 7.4 NON-DIMENSIONAL TORSION-SHEAR INTERACTION CURVES - TEST RESULTS 
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FIG. 7.5 NON-DIMENSTONAL TORSION-BENDING INTERACTION CURVES - TEST RESULTS 
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7.5. In non-dimensionalizing the forces the shear-flexural capacity of 


beams in Series C is taken the same as that of Series A. The only 
difference between these two series is eccentricity of prestress. 
Flexure-shear tests are made on beams with stirrup spacing at 3 inch 
centers and it is assumed that the variation of stirrup spacing does 
not affect the flexure shear strength. Pure torsional tests in 
Series C beams with stirrup spacing at 4 inches, and in Series B with 
stirrup spacing at 4 inches and 5 inches were not made. For interaction 
study, these strengths are obtained from pure torsion tests strengths 
of the same series with stirrups at different spacing, by adding or 
subtracting as the case may be, the difference in the contribution of 
stirrups towards the torsional strength as obtained by theory. 

It is evident from Fig. 7.2 that the torsion-shear interaction 
curves are remarkably similar to the torsion bending curves. This is 
partly because the control points on the bending and sheur axes are ob- 
tained by tests, and in flexure-shear tests shear failure did not occur 
prior to the attainment of flexural capacity whether the beams were 
under reinforced or over-reinforced flexurally. Examining Fig. 7.2, 
it becomes evident that with the exception of Series C having the lowest 
eccentricity, the torque capacity in combined loading increased to a 
peak value and did not fall below the pure torsional strength until the 
bending moment was about 40, 55 and 65 percent of the pure flexural 
strength in Series B, D and A respectively. This increase was about 
6 to 12 percent; beams in Series A with the highest eccentricity showed 


the highest increase. Regarding the magnitude of increase in the tor- 


sional strength these test results differ significantly with the ob- 
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servations of Bishara (1969) who reported the increases in torsional 
strength as much as 100 percent in the case of I beams with stirrups 
under combined loadings. In non dimensionalizing the moments Bishara 
used theoretical pure torsional strength. The bending strength was 
reduced only 6 and 15 percent of the pure bending capacity when the 
applied torque was 75 and 90 percent of the pure torsional strength for 
Series D and A. For Series B though designed to achieve a square inter- 
action curve this was not obtained due to premature failure of the 
compression zone. This premature failure was evident in varying degrees 
in other series. While the gain in torsional strength above pure 
torsional strength with small amounts of bending and shear was more 
gradual; the decline in the torsional strength beyond the peak with 


further increases in other forces was very sharp. 

For Series C, all the beams at load combination 3 had 6 to 
10 percent less prestress than the corresponding pure torsion beams. 
This might have been responsible for the absence of gain in the tor- 
sional strength above pure torque capacity for the beams at load combin- 
avron. se 

In Fig. 7.3, the theoretical interaction curves based on the 
theory are shown. The test strengths represented by HOA 08S are non- 
dimensionalized by the theoretical] Ue and Ble, The theoretical shear- 
flexural capacities of the beams A-3-0, B-3-0 and D-3-0 are 684.7, 691.2 
and 464.5 in. kips, respectively. These diagrams indicate the close 
predictions for Mode 1 and conservative predictions for Mode 2. 

Beams that failed in Mode 1 have magnitudes of bending moment 


and torque associated with the corner region away from the coordinate 
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axes of the nearly rectangular torque-moment interaction diagram. These 
test results are closely predicted by the theory of Mode 1 and the 
failure torques are of the same order as those associated with Mode 2 
failures. The theoretical torque in Mode 2 in some cases is slightly 
less than that in Mode 1 as the Mode 2 failures are conservatively 
predicted by about 20 to 30 percent on the average. Therefore it is 
necessary to determine at which T/Vb ratio Mode 1 changes to Mode 2 
failure. This transition occurs at a particular Ai ratio which must 
be established from the tests. Sufficient test results are not accumu- 
lated here to find the transitional T/V DS ratio. However the transition 
can be reasonably assumed to occur ata T/Vb,, ratio where the predicted 


bending moment in Mode 1 becomes greater than that in Mode 2. 


7.5 Application of Theory to Beams in Literature 

The theory presented in Chapter IV was applied to several I 
section beams reported in the literature. The test and theoretical 
Strengths of specimens failing in Mode 1 and Mode 2 are furnished in 
Table 7.6 and 7.7 respectively. If data reported in some cases con- 
cerning the stress-strain characteristics of strands was inadequate, 
missing information was obtained from standard stress-strain curves 
available from the publication by British Ropes Limited, 1967. 

The I beams of Gausel (1970), and Wyss and Mattock (1971) 
exhibiting Mode 1 failure were analyzed and the theory has predicted 
the test strengths within three percent on the average. i test/T theory 
ranged from 0.90 to 1.14. One half of Gausel's beams did not have web 


reinforcement in the test zone and the other half had. The correlation 
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Investigator 


Gausel 
(1970) 


Wyss and 
Mattock 
(1971) 


Beam 
No. 


TABLE 7 +6 


COMPARISON OF TEST AND THEORETICAL ULTIMATE 


ihe 
Cc 
(psi) 


8236 
7483 
7810 
7653 
8747 
7923 
7455 
7085 
7384 
6887 
7597 
7795 
8335 
8094 
6943 
6659 


7100 
5950 
7350 
7300 
7050 
6650 
6740 
7050 
7300 


Prestress 


force 


(kips) 


70. 
69. 
64. 
62. 
66. 
66. 
67. 
736 
70. 
69. 
64. 
62. 
66. 
66. 
67. 
67. 
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“All! 
220. 
224. 
230. 
224. 
223) 
AVE 
All 
218. 


84 
70 
54 
60 
08 
56 
92 


84 
70 
54 
60 
08 
56 
92 
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Eccentricity 


Size and 
spacing 
of stirrups 
in inch 
#2 @ 6.3 
#2 @ 6.3 
#2 @ 4.49 
#2 @ 4.49 
#2 @ 2.88 
#2 @ 2.88 
#2 @ 2.97 
#2 @ 2.97 
#4 @9 
#409 
#409 
#409 
#409 

#4 @ 6 

#4 @ 6 

#4 @ 6 


Test Results 
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42. 
55. 
byAc 
59. 
52% 
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(in-kips) 


2 
0 


My 


586. 
586. 
759). 
607. 
595. 
649. 
75 
Sih 
-20 
586. 
759. 
607. 
595). 
514. 
80 


541 
54] 
651 


54] 


476. 


6500. 
4108. 
4257. 
4520. 
4345. 
4206. 
Nese 
4232. 
3819. 


(in-kips) 


08 
08 
00 
20 
32 
44 
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00 
00 
30 
00 
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MOMENTS FAILING IN BENDING MODE (MODE 1) 


Theory 
uf M 
(in-kips) (in-kips) 

46.77 650.07 
51.23 548.20 

-- 767.82 
48.98 563.26 
54.49 544.89 
47.24 581.06 
64.90 538.69 
65.85 546.58 
43.50 669.93 
60.66 535.59 
46.43 668.63 
64.13 547.06 
[Eas 533.76 
78.27 529.87 
82.04 529.16 
84.44 Bie Gil 

Average 

0 6000.4 
353 3921.1 
386.6 4175.6 
385.7 4436.0 
400.1 4521.3 
402.8 3959.5 
434.8 4091.6 
446.4 4392.7 
347.8 3861.1 


Average 
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with the test strengths has been equally good in both cases. 

Wyss and Mattock reported the outstanding flanges split away 
from the girders at ultimate. Because the horizontal legs of the stir- 
rups were separate bars and did not have the full development length 
their contribution to the strength is neglected. However it can be 
noted that in Mode 1, the number of horizontal legs of the stirrups 
intersected by a failure crack is usually one or two and in the girders 
of Wyss et al none because of the higher inclination of the crack to 
the axis of the beam and the narrow web. Even if the intersected 
horizontal legs are assumed to yield and their strength contribution 
taken into account, the difference in strength is only of the order of 
] to 3 percent. In their case there is no contribution from the hori- 
zontal legs. For the reasons stated in Section 7.3.4, the outstanding 
flange parts were neglected in analyzing Gausel's beams. 

Several full scale I girders, prestressed and non prestressed, 
tested under pure torsion by Wyss, Garland and Mattock (1969) were 
also analysed. These beams failed in Mode 2. ie test/T theory varied 
from 0.81 to 1.37 for prestressed and reinforced girders. On the 


average the theory over predicted the strength by only 2 percent. 
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CHAPTER VIII 
SUMMARY AND CONCLUSIONS 


8.1 Summary 

This investigation was devoted to the study of the behavior 
of prestressed concrete I section beams loaded in combined torsion, 
bending and shear. The experimental phase of this study consisted of 
tests on 44 eccentrically prestressed I beams to evaluate the effect 
of eccentricity and level of prestress, varying combinations of torsion, 
bending and shear and various amounts of transverse reinforcement on 
the behavior and strength of these beams. The analytical phase of the 
study included strength predictions at cracking and ultimate of pre- 
stressed concrete I beams. Test results were discussed, and theory was 
applied to the beams of the present investigation and that of others. 

The following conclusions are drawn, based on the findings 


of the tests reported here. 


8.2 Conclusions 

|. Prestressed concrete I beams subjected to torsion, bending, 
and shear can develop initial cracks either on the top face, on the web, 
or on the bottom face. The location of the initial crack depends on the 
relative dimensions of the cross-section, on the magnitude and eccentricity 
of prestress, and on the relative amounts of torsion, bending and shear. 

2. An increase in prestress causes a corresponding increase in the 
cracking torque of prestressed I beams; initiation of cracking appears 
to be governed more by the total prestress rather than stress at a point. 
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3. Moderate amounts of bending moment increase the torsional 
cracking strength. This is greatest when the apptied bending moment is 
equal to the internal moment created by the eccentric prestress. Bending 
moment along with the corresponding shear in excess of an optimum amount 
results in rapid reduction of the torsional cracking strength. 

4. Prior to cracking the initial torsional stifness is nearly in- 
dependent of the amount of transverse reinforcement, and the magni tudes 
of torsion, bending and shear. However the secant torsional stiffness 
at cracking is higher for the beams loaded at a very low T/M ratio than 
for beams loaded at a high T/M ratio. Flexural stiffness of the uncracked 
section is unaffected by the presence of torsion. 

5. The three dimensional finite element analysis developed here 
using hexahedrons predicts conservatively the cracking strength of pres- 
tressed and reinforced concrete beams subjected to combined loadings. 
This method of analysis is more conservative for the case of web cracking 
in I beams than that for cracking at other locations on the cross-section. 
The theory estimates the test strength within 26 percent, on the average. 
The torsional and flexural stiffnesses at 70 percent of the test cracking 
strenath can be satisfactorily estimated. This analysis has reinforced 
the possibility of neglecting the intermediate principal tensile stress 
in formulating the failure criteria of concrete for slender beams. 

6. Assuming full plastic distribution of torsional shear stress 
and elastic distribution of stresses due to other forces, and using the 
principal stress failure criteria coupled with the tensile strength of 


concrete which is estimated equal to the splitting tensile strength, the 
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cracking strengths are conservatively predicted. This method also, 
like the finite element method, proves to be more conservative for the 
case of web cracking than that for cracking at other locations on the 
cross-section. The case of the web cracking is more closely predicted 
if the initial crack is assumed to occur at the centroid of the cross- 
section. The theory conservatively estimates the test strength within 
17 percent, on the average. 

7. Similar to reinforced and prestressed concrete rectangular 
beams, prestressed I beams can fail in skew bending either in a bending 
mode (Mode 1) or a torsion mode (Mode 2). The bending mode of failure 
is characterized by the formation of a compression hinge in the top 
flange joined by an approximately continuous tensile crack on the other 
three sides. The torsion mode (Mode 2) is recognized by a nearly con- 
tinuous tensile crack on the two flange faces and on one side of the web. 
Based on the type of crack patterns observed, it is concluded that the 
location of the compression hinge is close to one face of the web rather 
than along the sides of the flanges. This location is in line with the 
findings of Wyss, Garland and Mattock (1969) and in contradiction with 
that suggested by Gvozdev et al (1968) and by Woodhead and McMullen (1972). 

8. The failure mode need not be the same as the initial cracking 
mode; for example as in the case for beams of this investigation subjected 
to pure torsion. 

9. Isolated I beams with web reinforcement such as those tested 
in the present investigation in which the compression zone is not provided 


with any anchored binding reinforcement may fail in an over-reinforced 
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manner, even though such I beams are under-reinforced in pure torsion 

and in pure bending. The reason for such behavior resulting in pre- 

mature failure and lesser strength is the reduced effectiveness of the 
compressed flange peeier tna the shear-compression forces because the 
outstanding flange splits off due to torsionally induced cracking. However, 
it can be expected that such a premature failure can be avoided by pro- 
viding confining Peintercenent in the compression flange. 

10. The configuration of the web reinforcement used in the present 
investigation can adequately resist torsion. The vertical legs should 
continue across the full width of the flange. It is not necessary to 
use closed rectangular stirrups as a web reinforcement. 

11. The addition of web reinforcement increases the strength and 
ductility of beams under pure torsion and under combined loadings. 

12. The minimum amount of web reinforcement used in this investi- 
gation corresponds approximately to that suggested by Hsu (1968d) for 
reinforced concrete rectangular beams with similar width and depth 
dimension of the stirrup. This minimum is sufficient to enable the beam 
under pure torsion or combined loadings to fail at loads in excess of 
cracking loads. 

13. The ultimate strength of prestressed I beams subjected to 
combined loadings can now be predicted by the theory presented in 
Chapter IV which is based on a skew bending failure mechanism taking into 
account strain compatibility. The bending mode analysis presented here 
takes into account the splitting off of the outstanding compression flange 


due to torsional cracking and the resulting premature failure by neglecting 
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the flange tips. If the amount of torsion applied is small such a 
premature failure may not occur and the full flange width can be con- 
Sidered effective. For such cases if the outstanding flanges are 
neglected, the analysis yields conservative results. The agreement 
between the test and theoretical results is good (Tu test/Tu theory = 0.98) 
for bending mode failures. 

The torsion mode analysis can be applied only to beams rein- 
forced longitudinally and transversely and it predicts the test strength 
conservatively; the prediction of strength is reasonably good (Tu test/ 
Tu theory = 1.20) for beams under pure torsion, but rather conservative 
for combined loadings (Tu test/Tu theory = 1.32). 

14. Until a more accurate assumption concerning strain compatibility 
conditions of reinforcement is established, significant progress cannot 
be made in the skew bending method of analysis for partially over- 
reinforced beams failing in Mode 2. 

15. The post cracking torsional and flexural stiffnesses depend 
primarily on the loading ratio and also on the amount of reinforcement 
and eccentricity of prestress. 

16. In web reinforced beams as long as the web reinforcement is 
sufficient to prevent failure at cracking, the torsional and the flexural 
ductilities do not vary significantly as the web reinforcement is varied. 

17. Torsional and flexural ductilities depend mainly on the loading 
ratios. Rotation capacity at ultimate reduces as the T/M ratio decreases 


or level of prestress increases. 
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APPENDIX A 


A.1  Torque-Twist and Load-Deflection Curves 

Torque versus angle of twist and load versus deflection 
curves for beams feted in the experimental phase of the program are 
presented in Fig. Al through A20. Each figure contains curves for 
several beams and for clarity the origin of each curve is offset by 
a fixed magnitude. Torque-twist and load-deflection curves for each 
beam represent the entire range of loading. In some cases the 
curves during the last increment of loading are shown as broken lines 
since measurements could not be available. In the load-deflection 
curves of a few beams, the last point on the curve represents the 
twisting of the beam without any increment in bending load. 
A.2 Typical Failure Crack Patterns 

Sketches of typical failure crack patterns for the beams 
in the form of developed surfaces are shown in Fig. A21 through A32. 
The location of the transverse load is shown by cross hatching on the 
top face of beam; this serves to identify the location of the failure 
surface relative to the transverse load. 
A.3 Photographs of Beams 

The photographs of the beams indicating the development 
of the crack patterns during the test are shown in Plate A.1 
through Plate A.7 . For each beam views of the sourth and north 
faces on which shear stresses due to torsion and shear are 


respectively additive and subtractive are presented. 
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APPENDIX B 
THREE DIMENSIONAL FINITE ELEMENT ANALYSIS 


B.1 General 

This computer program relates to a three dimensional finite 
elemert analysis for predicting the cracking strength of concrete beams 
subject to combined loadings. A general description of the program, 


formulation of element stiffness, and a listing are given below. 


B.2 Usage and Limitations 

The program was written in the Fortran IV Language for the 
IBM 360-67 installation operating under MTS at the University of Alberta, 
Edmonton. It is hardware dependent and will provide a solution for 


problems that do not exceed the following size. 


Number of elements 560 
Number of nodes 950 
Number of materials ] 


The limits can be easily extended by changing the Dimension statements 
provided the material used is the same. Even small three dimensional 
problems require large storage and computation time. In such problems, 
the system stiffness matrix requires most of the share of the total 
storage. Economy in the storage of the stiffness matrix can be achieved 
recognizing its banded nature and symmetry. Further economy can be 
effected by generating the stiffness elements in blocks and transferring 


the data to auxiliary storage keeping only two blocks in core for current 
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computations. In the present program, only two fixed size square blocks 
are in core, each block having a size of one-half band width. This 
method hence requires only (2 half bandwidth {half band width + 1}) 
words of in-core storage. Thus there is only a limit on the half band 
width, and the total number of equations that can be solved is unlimited. 
Using this program a problem having a one-half band width up to 300 can 
be solved by the IBM 360-67 computer with a usable memory of 1024 k. 
However with large systems of equations round-off starts to play an 
important role. A clear remedy is to do the computations in double 
precision. But this roughly doubles the core storage requirement. 

Above all there is a limit on the time of computation, one would not 
like to exceed. Hence it is recommended that the band width be limited 


to 200, with a maximum number of equations. 


B.3 Salient Features 

A two dimensional finite element program with a block solver 
for solving equations, developed by Wilson (1966), and an equation solver 
coded by Gopalakrishnaiah (1973), formed the basis in the development 
of this program. An eight node hexahedron is used in this program for 
the three dimensional cracking analysis of reinforced and prestressed 
concrete beams. The same elelent is specialized for a triangular prism. 
B.3.1 Improvement of Bending Performance 

Before use is made of any element in an analysis it is desirable 
to test its performance in structural idealizations made of beam type 


structures. Clough (1969) has shown that eight noded hexahedra cannot 


adequately represent plate bending mechanism though its performance is 
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good in short deep beams. However, based on comparison of several ele- 
ments he recommended an eight node element for analysis of general 
elastic solids. Gopalakrishnayya (1973) used eight node hexahedra in 
the analysis of the cracking of earth dams. His formulation, although 
not meant for bending problems, has poor bending characteristics for 
beams of usual proportions even at an aspect ratio of about 2.5. 

Hence as suggested by Wilson (1966), the general hexahedron 
element is modified here to improve its bending behavior. The modifi- 
cation consists of the addition of extra degrees of freedom, within 
the element. This causes a violation of inter-element displacement 
compatibility. The generalized displacements associated with the extra 
degrees of freedom within the element are eliminated at the element 
level after the element stiffness is formed by static condensation. The 
resulting hexahedron has significantly improved bending characteristics. 
B.3.2 Data Transfer in Solver 

The structure stiffness matrix is generated in blocks and 
solved in blocks, keeping only two consecutive blocks at any time in 
the core; economizing the requirement of core storage. This process 
requires transfer of data from sequential disc files to core and vice- 
versa. An efficient method of data transfer between core and disc 
files using certain system oriented subroutines and appropriate equiva- 
lent statements resulting in considerable savings in the cost of com- 
putation has been available (Soil Mechanics report number 20, University 


of Alberta, 1973) and is used in the program. 
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B.3.3 Element Stiffness Formulation 


The stiffness is the same for elements of same size and orien- 
tation provided the maerial is same, wherever they occur in the model. 
Further economy in computation time is achieved by formulating the 
element stiffness only once for a group of identical elements and storing 
it in a sequential disc file for subsequent use by an identical element. 
B.3.4 Solution of Different Load Vectors 

If more than one loading condition is to be solved economically 
it is not necessary to invert the stiffness matrix for every load vector. 
Hence, the reduced matrix is stored in a temporary sequential disc file 
and used every time for reducing the new load vector. This facilitates 


solution of several problems with inversion made only once. 


B.4 Description of Subroutines 

The computer program uses eleven subroutines with the main 
program. Further, it requires the following system-oriented subroutines 
in connection with data transfer between core and auxiliary storage. 

RCALL, SETDSN, ADROF, WRITE, READ, NOTE, POINT 

Additionally, another system subroutine TIME computes and prints the 
CPU time. The details of these system subroutines can be obtained from 
the MTS system manuals. A brief description and purpose of each of the 
eleven subroutines used, follows. Figure B.1 shows the order of calling 
the subroutines. 
MAIN PROGRAM 

The main program calls subroutine THREED. The dimensions of 


certain arrays depending on the one-half band width and size of the 
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THREED 
READIN ASTIF BANDI STRESS 
| 2 s) 6 
ELSTIF MODIFY ERSTE PRINCI 
3 4 if 8 


MAIN 


THREED : 


READIN 
ASTIF 
EST tk 


BANDI 
STRESS 
ESA 
PRINCI 
POLRT 
EIGEN 


Note 


> Changes 


> Reads and outputs nodal, element, 
> Forms 


ELGs Bel 


dimension 
statements as per size 
of program 


POLRTM EIGEN 
9 10 


material data 


Calls all subroutines 


system. stiffness 


> Assembles element stiffness 
MODIFY : 


Modifies displacement boundary conditions 
> Solves equations 
: Prints stresses 
: Computes stresses 
| Calls (ROLRI sand -eEIGEN 
> Computes principal stresses 
> Computes principal stress directions 


Numbers represent the order of execution of 
subroutines 
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problem are supplied to other subroutines. 
THREED 

This subroutine calls all the other subroutines, computes 
modulus matrix and the coefficients necessary for numerical integration. 
It further reads load vectors for each problem, if more than one loading 
condition is studied. 
READIN 

This subroutine reads the nodal data, element data, material 
properties and the first loading condition. It computes the one-half 
band width. READIN incorporates a data generating facility not only 
across the cross section but also along the length if needed; only a 
minimum amount of information, is then needed as input to specify the 
problem topology and geometrics. 
ASTIF 

This subroutine assembles the load vector and the system 
stiffness matrix after obtaining element stiffness from subroutine ELSTIF. 
The assembly is through a direct stiffness approach. Element stiffness 
contributions are directly added to the structure stiffness by using a 
simple labelling technique in blocks of fixed size (half band width) x 
(half band width + 1). Subroutine MODIFY, called for each node, imple- 
ments the displacement boundary conditions in blocks. As the blocks are 
formed, they are stored in a temporary sequential disc file (file number 2). 
At the time of formation and subsequent solution of the stiffness equations, 
only two consecutive blocks are needed in the core at any time. Hence 
storing of other blocks in auxiliary storage, on needed for the current 


computation saves core space, allowing the solution of larger problems. 
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EESTI 

This subroutine formulates the element stiffness matrix and 
stores it in a temporary sequential disc file (file Number 3). The ASTIF 
subroutine reads element stiffness from this disc file. The subroutine 
is skipped if it has already formed once the stiffness of an element 
having the same size and orientation. It further calculates the element 
and nodal stresses and strains when called by subroutine STRESS. 

MODIFY 

This subroutine modifies the stiffness matrix and the load 
vector in accordance with the displacement boundary conditions. 
BANDI 

This subroutine solves in single precision, the modified 
stiffness matrix equations in blocks, by the method of direct Gauss 
elimination. For reduction, or back substitution, only two consecutive 
blocks of the stiffness matrix are required in the core. Blocks are 
transferred from disc files to core and vice-versa when required. For 
reduction, the structure stiffness matrix is obtained from file Number 
2 as formed in ASTIF. The reduced equations used for back substitution 
are stored in file Number 1. 

If more than one load condition is to be solved, it is not 
necessary to invert the matrix each time. The same reduced matrix 
available in file Number 1 is stored once for all the load conditions 
in file Number 2 where, previously unreduced stiffness elements were 
available. For solving each new loading condition the reduced matrix 
is read from file Number 2, computations using the new load vector are 


completed and the equations are then stored in file Number 1. Back 
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substitution is carried on from file Number 1 as before. 
STRESS | 

The subroutine computes stresses and strains by calling sub- 
routine ELSTIF. It also calls subroutine PRINCI for calculating principal 
stresses and their directions. 
PRINCI 

This subroutine calls subroutines POLRT and EIGEN for the com- 
putation of principal stresses and their directions. 
POLRT 

This is a subroutine of the IBM SS package and calculates the 
principal stresses, by reading element and nodal stresses computed in 
ELSIE: 
EIGEN 

This is another subroutine from the IBM SS package and computes 


the principal stress directions. 


B.5 Changes Necessary in the Main Program Before Execution 

The dimension declaration of certain arrays in the program is 
to be in accordance with the one-half band width and number of equations 
of the problem to be solved. This can be accomplished in the main program 
without recourse to the rest of the program once these dimension values 
are known exactly. Because of a flexibility provided, the program can 
be terminated after the generation and output of the geometry of the 
problem, number of equations, and one-half band width. This allows to 
check the input data and gives the values of the one-half band width and 
number of equations. To execute up to this stage, some arbitrary values 


can be assigned to the one-half band width, number of equations, and 
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program run. Once the exact values are known, dimension and equivalent 
statements are set accordingly in the main program and execution is 
allowed to continue until the end. The statements requiring a change 
in the main program are detailed below: 


MB1 


Half band width (MBAND) 

NLD = MBAND X Number of blocks 
DIMENSION B(NLD), A(MBAND, 2 X MBAND), BL(MBAND), BR(MBAND), AL(MBAND X 
MBAND), AR(MBAND X MBAND) 
EQUIVALENCE (B(1), BL(1)), (B(MBAND+1), BR(1)), (A(1,1), AL(1)), 
(A(1,MBAND+1), AR(1)) 


B.6 Gauss Quadrature, Load Conditions, and Output 

The computer program permits the use of both a hexahedron or 
a triangular prism which is specialized from the hexahedran. Two point 
Gauss quadrature is used in the integration of the element stiffness 
matrix. In the case of a triangular prism, as it is formed by combining 
two pairs of nodes of hexahedron, to avoid quantities allowing to tend 
to infinity, stresses and strains are computed very close to the nodes 
but not exactly at the nodes. 

The program provides for the solution of any number of load 
conditions with the inversion of stiffness matrix accomplished only 
once. Nodal code value - 1 is to be specified for the nodes where 
forces are specified for different load conditions. By supplying unit 
weight of the material, the self weight of the beam can be taken into 
account. 

The output from the program consists of a reprint of all the 


input data, nodal displacements, and element and average nodal stresses 
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and strains. It also includes element maximum shear stress, principal 
stresses and strains and their directions at the centre, on the re- 


quired faces, and at each node. 


B.7 Estimation of Computation Time 

It is advantageous to estimate in advance the time required 
for the solution of a problem. The computation time essentially con- 
sists of amounts of time required for the calculation of element stiff- 
ness, assembly of structure stiffness, solution of equilibrium equations 
and computation of stresses and strains. Once the number of elements, 
nodes and one-half band width of a model are known, the computation 
time can be reasonably evaluated. 

Time required for one element stiffness formulation is approxi- 
mately 2.5 seconds. As stiffness is not formed for identical elements, 
the total time for stiffness formulation may be obtained by multiplying 
the number of groups of dissimilar elements by 2.5 seconds. Time required 
for stress calculation can be arrived at by multiplying the number of 
nodes by about 0.5 seconds. Based on the solution of several problems 
of different sizes (different one-half band widths), it has been possible 
to draw a relationship between the average computation time required for 
the solution of equations in one block and the one-half band width. This 
relationship is furnished in Fig. B.4. As the solution time for a 
given number of equations and band width depends on the number of dis- 
placement boundary conditions implemented, the time predicted by means of 
Fig.B.4is approximate in a situation different from the ones used to 
plot the curve in this figure. Total solution time of equations may be 


approximated as the product of the number of blocks and the time required 
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FIG. B.4 APPROXIMATE CPU TIME FOR THE 
SOLUTION OF EACH BLOCK OF EQUATIONS 
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for one block as obtained from PiGas By 4s 


B.8 


to represent its boundary displacement condition. 


Code for Nodes 


Zoe 


Each node is assigned a code consisting of a two digit number 


The codes used to 


identify various boundary displacement conditions are detailed below: 


Code 
00 
01 
02 
03 
04 


X force 


or displacement 


X-force 
X-force 
X-force 
X-displacement 
X-displacement 
X-displacement 
X-force 


X-displacement 


Y force 


or displacement 


Y-force 
Y-force 
Y-displacement 
Y-force 
Y-force 
Y-displacement 
Y-displacement 


Y-displacement 


B.9 Code for Face of an Element 


Z force 
or displacement 


Z-force 
Z-displacement 
Z-force 
Z-force 
Z-displacement 
Z-force 
Z-displacement 


Z-displacement 


Each element is provided with a code consisting of a two digit 


number to identify the faces on whose centres stresses are to be computed. 


A maximum of two faces can be specified. 


sents a face. 


Each digit of the code repre- 
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Code Faces 
0 Centre of element 
] X-coordinate positive face 
rs Y-coordinate positive face 
S X-coordinate negative face 
4 Y-coordinate negative face 
3) Z-coordinate negative face 
6 Z-coordinate positive face 


B.10 Nodal and Element Data Generation 

Automatic data generation within the program reduces the 
amount of data to be input. Nodal data is to be fed in ascending order. 
If a certain series of nodes are equally spaced in a straight line and 
have their numbers in increasing order, intermediate nodes can be auto- 
matically generated with their nodal code and nodal displacements or 
loads equal to zero, by supplying only the two nodes at the beginning 
and end of that part. 

For example as shown in Fig. B.5 9 to 13 are equally spaced 
in the Z-direction, and nodes 13 to 16 in both the x and z-directions. 
Each extent is a straight line. All nodes from 9 to 16 can be generated 
by giving information only for nodes 9, 13 and 16. If any intermediate 
nodes 10, 11, 12, 14, or 15 have nodal codes or displacements or forces 
different from zero, they can be taken care of individually at a later 
stage. Nodes 24 to 30 cannot be generated in one sequence and have to 
be read individually. 


The elements can be read in at random one by one, however, 
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CROSS SSEGIION.. 2] CROSS SEC UON 2 


Note : Element numbers are circled 


CROSS CROSS CROSS CROSS 
SECTION SECTION SECTION SECTION 
1 2 3 4 


FIG. B.5 A GENERAL LAYOUT OF ELEMENTS AND NODES 
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it is preferable to do it in order. If a group of consecutive elements 
have their numbers and their corresponding nodal numbers in increasing 
order, the element data of all the elements in the group can be generated 
by the computer by feeding in only the final element. All the elements 
created in one extent assume the same material number, element code, 
element size code, and element face code, as the first element of the 
group. If any elements, in the process of automatic generation assume 
any of these code values different from their actual codes, the proper 
values can be fed in individually. The variable "NUMTIE" represents 
the total number of groups of elements in the model each group being 
fed in one extent. If all the elements are input individually NUMTIE 
will equal the total number of elements in the model. Again, referring 
to Fig. B.5, elements 1 to 7 can be generated by giving only the data 
of element Number 1. All the elements 2 to 7 have the same material 
number, element code, element size code, and element face code as 
element Number 1. However, elements 13 to 19 have to be read in in- 
dividually as the nodes from 54 to 60 are not in increasing order. 

If the nodes are spaced in the same way across all the cross 
sections, and if the element properties are the same throughout the 
length as of the elements enclosed by the first and second cross sections, 
then by providing the nodal and element data of first cross section alone, 
the properties of the whole beam can be generated. This facility reduces 
the input data to approximately five percent of the total required if 


the beam is modelled using 20 tiers of elements along its length. 
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B.11 Size and Creation of Disc Files 


A total of three temporary sequential disc files of sufficient 
size are to be created for the execution of the program. Two files, 
one for storing the structure stiffness matrix (file Number 2), the 
other for storing reduced stiffness matrix (file Number 1), are needed. 
The third file (file Number 3), stores the element stiffness matrix, 
and HMUL5, a matrix generated in connection with static condensation. 
The size of a file can be expressed in terms of the number of tracks 
required. Each track cf the file corresponds to about 7000 bytes. 
The number of tracks needed for the first and second files are given 

(Number of blocks) x 


Number of tracks = 


(half band width)(half band width + 1) (ha) 


The number of tracks for file Number 3 can be calculated as follows: 


Number of 

groups of 4 
Number of tracks = (24 x 24 +9 x 24) |dissimilar (=00) 

elements 


The following control cards are necessary for the creation of the files: 
$CREATE $ - T & TYPE = SEQ $ SIZE = nT 
$CREATE J - TEMP26 TYPE = SEQ # SIZE = nT 


$CREATE & - T1 ~ TYPE = SEQ $ SIZE = nT 


$RUN 6 - LOAD # 6 1 = - T B 2 = - TEMP26 3 = - Tl 


where "n" represents the number of tracks. 


B.12 Variables Associated with Size of Problem 
The arrays whose dimensions depend on the size of problem, 


in terms of number of elements, number of nodes, and number of element 
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size code, are given below. These dimensions as declared in the program 


need not be changed if the following limits are not exceeded: 


Number nodes = 950 
Number of elements = 560 
Number of element size code = 100 


In case they are exceeded, the common or dimension statements of the 
arrays are to be accordingly specified. Further, statement 24 in sub- 


routine ASTIF is to be changed in accordance with the maximum element 


size code. 
Variable Description Minimum Size 
KODEL code assigned for each element (NUMEL ) 
to identify as hexahedron or 
triangular prism 
KELLS code for specifying the faces (NUMEL ) 
of each element over which 
stresses are required 
KOST code for size of each element (NUMEL ) 
for formulating the element 
stiffness matrix 
KODE code for each node to identify (NUMNP ) 
its boundary displacement con- 
dition 
KOUNT counter used for averaging nodal (NUMNP ) 
stresses and strains 
LCHECK variable used to determine whether (as many as the 
element stiffness is to be formed number of groups 
or already has been done of dissimilar 
elements) 
MB1 one-half band width 
MAT material number alloted for each (NUMEL ) 


element 
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NLD 


NP 


SIGPC 


STGP] 


SIGP2 


SIGA 


STNODE 
U 


WEIGHT 


number of blocks x one-half 
band width 


vector of eight nodes stored in 
order of each element 


vector of element principal 
stresses, strains, and maximum 
shear at the centre of an 
element 


vector of element principal 
stresses, strains, and maximum 
shear at the centre of face 1 
vector of element principal 
stresses, strains, and maximum 
shear at the centre of face 2 
nodal stress vector 

nodal strain vector 

force or displacement at each 
node in the direction of x- 
coordinate 

force or displacement at each 
node in the direction of y- 
coordinate 

force or displacement at each 
node in the direction of z- 
coordinate 


weight at each of the 8 nodes 
of an element 


x-coordinate of a node 
y-coordinate of a node 


z-coordinate of a node 
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(8, NUMEL) 


(NUMEL ,7) 


(NUMEL ,7 ) 


(NUMEL ,7) 


(NUMNP ,6) 
(NUMNP ,3) 
(number of blocks 


X MBAND/3) 


(number of blocks 
X MBAND/3) 


(number of blocks 
X MBAND/3) 


(as many as number 
of groups of 
dissimilar 
elements ) 
(NUMNP ) 
(NUMNP ) 


(NUMNP ) 
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B.13 Input Data Format 
General Data 


1. Read: HED, NUMNP, NUMEL, NUMAT, NUMTIE, NUTRPR, NLODC (18A4/616) - 


two cards 

HED = Heading 

NUMNP = Number of nodal points 

NUMEL = Number of elements 

NUMAT = Number of materials 

NUMTIE = Number of series of hexahedra in the model 
If each element is read in separately, this 
will be equal. to the total number of elements 
(NUMEL ) 

NUTRPR = Number of triangular prisms in the model 

NLODC = Number of total load conditions 


2. Réad: E, PR, RO (Fl2e0sdF6233 7.3) one card 


E = Modulus of elasticity 
PR = Poissons ratio 
RO = Unit weight (If RO is non zero, the self 


weight of the material is included in the 
negative z-direction 


3. Read: NRSEPA (14) - one card If equal to zero, only the nodal 
and element data on the first cross section will be read 
and the information for the rest of the model will automati- 
cally be generated. If equal to non zero, nodal and element 
data has to be read at every cross section along the length 
of the model. 


NRSEPA = This variable is used to determine how data is 
going to be read. 
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4. Read: NREP, NSUMN, NS1 (314) - one card If NRSEPA # 0, this 
read statement is skipped 


NREP = one less than the number of nodal sections along 
the length of model 


NSUMN = Number of nodes in the cross sections 
NS] = Number of elements in one tier along the length 
Nodal Data 


5. Read: N, KODE(N), X(N). Y(N),.Z(N), UCN); VIN), WIN), (215, 6F5.0) 
Number of cards depends on the problem and one of the two ways 


nodal data is fed. One card for each node. 


N = Nodal number 

KODE (N) = Nodal kode (Refer Sec. B.8) 

X(N) = x-coordinate of node 

Y(N) = y-coordinate of node 

Z(N) = z-coordinate of node 

U(N) = force or displacement at node in the x-coordinate 
direction 

V(N) = force or displacement at node in the y-coordinate 
direction 

W(N) = force or displacement at node in the z-coordinate 
direction 


If NRSEPA = 0, nodes on one cross section are generated by 
the execution of Read 5. If NRSEPA # 0, by this all the 
nodes of the model are generated. 

6. Read: YIN (F5.0), As many cards as NREP 


YIN = y-coordinate of the nodal cross section along the 
length, starting from the second cross section 


If NRSEPA = 0, execution of Read 6, generates all the nodes 


of the beam model. 
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Read: NKOCH] (14) - one card 
If NKOCH] = 0, earlier nodal code generation does not 
require a change. And Read 8 and 9 are skipped. 
Read: NKOCH2 (14) - As many cards as NKOCHI 


NKOCH2 = one less than the number of nodes where 
nodal code is to be changed in one extent 


Read: N, KODE(N) (214) - As many cards as NKOCHI. 
N = node number 


KODE(N) = nodal code 


Element Data 


10. 


Tle 


Read: NUMEE (14) - As many cards as NUMTIE 


NUMEE = one less than the number of hexahedra in each 
series 


If NUMEE = 0, only one element will be read 
Read: M, NP(1,M), NP(2,M), NP(3,M), NP(4,M), NP(5,M), NP(6,M) 
NP(7,M), NP(8,M), MAT(M), KODEL(M). KOST(M), KELLS(M), 
(1315) - As many cards as NUMTIE. 


M = element number 

NP(I,M) = nodal numbers at the eight corners of an 
element in counter clockwise direction 

MAT (M) = material number 

KODEL (M) = type of element. If equal to zero, hexahedron. 


If equal to one, for triangular prism 
KOST(M) = code for the size and orientation of an element 
Elements having same size and orientation have same code 
for size, for purposes of calculation of element stiffness. 
KELLS(M) = Element face code. Two digit number, each 
digit identifying a particular face on the 


centre of which stresses are to be computed 
(refer Sec. B.9). 
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M, NP(I,M), MAT(M), KODEL(M), KOST(M), KELLS(M), (1315), 


As many cards as NUTRPR 


If NUTRPR = 0 as given in Read 1, then Read 12 is skipped. 


Read: 


NELFC (14) - one card 


If NELFC = 0, Read 14 is skipped. 


NELFC = number of elements if any having face code 
different from what is generated earlier 

Read: M, KELLS(M), (214) - As many cards as NELFC. 

M = element number 

KELLS(M) = face code of element 

Read: NELF (14) - one card 

NELF = number of elements if any having element 
code size different from what is generated 
earlier 

Read: M, KOST(M), (214) - As many cards as NELF 

M = element number 

KOST (M) = element size code 


If NRSEPA # 0, data of elements is generated by the time Read 16 


is executed, if NRSEPA = 0, element data on the first cross 


section only is read. 


Read: 
KADD 


KADD (14) - as many as (NREP-1) 


= A number to be added to each of the tiers of 
elements along the length, to specify the 
element size code. 


If all the tiers along the length of beam have 
the same length in the y-direction, then KADD 
for all the tiers is equal to zero. 
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First Load Condition 


ites 


D2 


20. 


Read: NOCHNO (I4) - one card 


| NOCHNO = number of nodes where nodal forces are to be 


assigned 
Read: N, KODE(N), U(N), V(N), W(N) (214, 3F6.0) - as many cards 
as NOCHNO 


N = node number 
KODE(N) = nodal code 
U(N) = x-force at node 
V(N) = y-force at node 
W(N) = z-force at node 


Read: NPARA (14) - one card 

NPARA = If > 0 execution terminates after the generation 
of nodal and element data, semi band width and 
number of equations. 


If NPARA = 0 solution is continued. 


Other Load Conditions 


Alls 


aap 


Read: NCHNO (14) - one card. As many times as number of load 


conditions - 1]. 


NCHNO = number of nodes where nodal forces are to be 
assigned for the particular load condition 


Read: N,KODE(N),U(N), V(N), W(N), (214, 3F6.0) - As many cards as 


NCHNO. As many times as number of (load conditions - 1). 


N = node number 

KODE (N) = nodal code (-1) 
U(N) = x-force at node 
V(N) = y-force at node 


z-force at node 
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B.14 Element Stiffness of Isoparametric Hexahedron 

An eight node hexahedron is used for the three dimensional 
cracking analysis of reinforced and prestressed concrete beams, in 
this investigation. It is based on the element developed by Irons 
and Zienkiewich. The element performance in bending is greatly improved 
by the addition of incompatible deformations as suggested by E.L. 
Wilson (1970). 

An arbitrary hexahedron is shown in Fig. B.6. The nodes 
of the element can be described by the three axes (x,y,z) constituting 
an orthogonal right handed system. This is referred to as the Global 
Coordinate System. A local coordinate system (also designated as 
natural coordinates) (&,n,7) is adopted, in which, for this element, 
each pair of opposite faces have values of +] and -1] on an axis. Thus 
in local coordinates the solid is a cube while it may be a distorted 
and warped six faced solid in space. The local coordinate system 
generalizes and simplifies the formulation. It further facilitates 
the integration necessary to obtain the element stiffness. 

The local coordinates are related to global coordinates 


by a set of linear interpolation functions. 


X P0450 X 
Kee oeeren i (B.1) 
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where P = <P, Po P. Pa Ns Pe Po Po> 


Pas Pos P. etc. are interpolation functions 
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(a) GLOBAL COORDINATES 


(b) LOCAL COORDINATES 


FIG. B.6 THREE DIMENSIONAL COORDINATE SYSTEMS 
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X12 Vz> 24> etc. are coordinates of nodes in Global System. 


The interpolation functions are in terms of local coordinates of the 


element. 
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In an isoparametric element the same interpolation functions 


are used to describe the geometry and the displacements of the element. 
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Accordingly the relationship between the nodal displacement vectors 
and the element displacements may be written, assuming displacements 


in the element vary linearly between the nodes. 


u P 0 0 u 
Vv [= Prat) Vv (8.2) 
Ww Ut0e PB W 
where mL =< Usei eee eae ee 
eine aoe ee Oe ey wcO 
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| ih uo Me Va Ve 6 Va Nog 


We = <Wy Wo W3 Wa We We Wy Nog 


u, V, W are nodal displacement vectors. 


To improve the bending characteristics of the element, the displace- 


ments within the element are assumed to be of the following form: 
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Pig (1-n ) 
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fo =e eel cee 
m mi m2 “m3 
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Wis = <Wen WOW 
m ml Wn2 m3 
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ee are displacement amplitudes associated with the extra 
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degrees of freedom. 


For a three dimensional body the strain components and the 


derivatives of displacements are related as follows: 
To derive stiffness matrix define strains in terms of nodal displacements. 
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{e} = [B]{r,} (B. 4b) 


ee is the nodal displacement vector including displace- 
ment amplitudes associated with extra degrees of 
freedom 


[B] is the strain-displacement matrix 


As the displacements U, V, W are not expressed directly in terms of 
global coordinates x, y, z the chain rule of differentiation has 


to be adopted. The chain rule in matrix form is expressed as: 
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where [J] is called the Jacobian matrix. 
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By equation (B.6) the derivatives of displacements U, V, W with re- 
spect to global coordinates can be expressed in terms of natural co- 
ordinate derivatives of the displacements. 

From the principle of minimum potential energy the stiff- 


ness matrix may be expressed as: 
K = | B'CB dV (BS?) 


where the matrix [C] is the stress-strain relationship for an elastic 


solid of isotropic material. 
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Thus once the matrices [B] and [C] are known, the element stiffness 


can be formulated. As the matrix (B) is expressed in natural co- 
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ordinates the integration is to be carried out in the same coordinates. 


K = fi [ [ B'cB dV (B.9) 


using the relationship dV = dxdydz = |J| d& dn dz. 

As it is impossible to perform the integration explicitly 
in equation (B.9), numerical integration is employed. In carrying 
out numerical integration by Gauss quadrature the matrices [J], [B], 
and the argument in equation (B.9) are evaluated at specific set of 
values of coordinates €, n, 3. The number of integration points re- 
quired in Gauss quadrature to achieve satisfactory results depends 
on shape of the element and finess of the mesh. In this investiga- 
tion two point quadrature formula is employed. It is programmed to 
obtain the matrix product B'cB detJ for each of the integration points 
and they are appropriately combined to form 33 X 33 element stiffness 
matrix. The number of integration on points required depends on the 


Shape of the element and on the finess of the mesh. 


B.15 Static Condensation 
The force-displacement relationship of the element may be 


expressed as 


(FJ = (kIr] (B. 10a) 
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ym? Eoin are the external forces associated with the 


generalized displacements Uns V 


where Fn? F 


<0 Wa Since there are no external 


forces associated with the displacements Uns V.,W., they may be ex- 


moe 
pressed in terms of the eight corner nodal displacements. From equation 


(B.10b), it can be shown 
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Now the nodal forces and the displacements are related as follows: 
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The final element stiffness matrix is 24X24 in size. 
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B.16 Evaluation of Stress 


The stresses at any point are given by: 
[o] = (Clle] = (cltetr,] (B. 13a) 
If [B] is partitioned 


[o] = [C][B, 8.) 


(B.13b) 


The stress may be expressed in terms of Conner displacements 


[o] = (C](8*)Lr] 


A: -] 
where [B*] = [B.-B Kay 


Kpad 
The generalized displacements {U, Va? W calculated at the time of 
element stiffness formulation, are saved and used for the evaluation 


stresses. 


{. | Pore eat _ 
abs ee 
at | pen % ory 


a er 
emaraosigeth vonnoa Yo eames at bazastqxs od am a2ende 


oe 


- 
_ 


Cea alah = [a] eat 

le= | il vee | 

to anf? ot tp perefuales 1,0 .,V ..0) etnemeosigeth bestievsnep oT 
“nob?outsve odd wo? beey bne bavee ave ,notsefumey canqtite snomets | 


[x}0*8909) = fo) 


274 
B.17 Theoretical Cracking Load 
The stress components for the three coordinate planes and 


the principal stresses are related by the following equation. 
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where Oys Ov Os Tyy? Ty? Tyy are the six stress components, in the 
present case these yield the stress distribution at cracking load with 
the effect of prestress eliminated making use of principle of super- 
position after the analysis is run. 

S = Principal tensile stress 

If the load is increased in a proportion f, the stresses also in- 


crease in the same proportion assuming the linear elasticity. 


Equation (B.14 ) may be written as: 
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where f = load factor and also stress factor. If prestress P is to be 
Superimposed on this stress distribution it adds to the component Oy 
which is the stress in the longitudinal direction. 


3 2 2 
S~ - Flo toto, \eSe = PS (o,0,t00,+ Oo, 


Cisare waa 2 2 
“Tyo tye Txy) f°S + (o,+0,) f°SP 


Z 
- (9, G22 ToT xy OxT yz 


-0.T _-0 x ) rome (ao see ) pre = 0 (B.16) 


From this equation for the given values of the stress components, 
prestress and a maximum principal tensile stress equal to the tensile 
strength of concrete, the factor f can be obtained. This factor gives 
the proportion in which the applied load associated with the stress 
components ys Oy etc is to be increased or decreased to yield a 


principal tensile stress equivalent to the tensile strength of con- 


crete at a given prestress, and hence the theoretical load. 
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